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Design data for free-air chambers measuring cobalt-60 and cesium-137 gamma ray 
It has been shown that the 
the saturation current is adequate for air pressures of about 4 to 12 atmospheres. 


roentgens are presented. 


s in 
Jaffé-Zanstra method of obtaining 


Also, 


radiation measurements of the gamma rays from cobalt-60 and cesium-137 made by a cavity 
chamber and a free-air chamber agree to within the experimental errors. 


1. Introduction 


The standard instrument for measurement of 5 to 
500 kv X-rays in roentgens is the free-air chamber 
[1]... Both free-air chambers [2, 3, 4] and cavity 
chambers are used for measurement in roentgens 
of the gamma rays from radium, cobalt-60, and 
cesium-137. Kaye and Binks [2] used both methods 
to measure the same source of radium and obtained 
agreement to the order of 1 percent between the two 
results. However, Taylor and Singer [3] have since 
pointed out that the free-air chamber measurements 
by Kaye and Binks were too low by the order of 
2 percent because of inadequate plate separation. 
Ionization from scattered photons, for which no 
correction was made, may have compensated for 
part of this error. This contribution was not meas- 
ured either by Kaye and Binks or by Taylor and 
Singer. On the other hand Attix and Ritz [5] have 
indicated that the cavity-chamber measurements of 
Kaye and Binks require a 1}4 percent correction for 
chamber wall attenuation. An additional possible 
source of error, inadequate guarding of the collector 
region, has been pointed out [6] for free-air chamber 
measurement of low-energy X-rays. Preliminary 
estimates of this effect on the radium [3], cobalt- 
60 [4], and cesium-137 [4] measurements indicate 
that the dose rates given there may be up to 5 per- 
cent low. Thus, there is not yet an unambiguous 
comparison of the two methods. 

At the present time, neither type of instrument is 
designated [1] as the standard instrument in the 
million-volt region. The ultimate selection of one, 
largely hinges on their individual capability of abso- 
lute accuracy. Therefore, measurement of the ex- 
posure dose rate of the same radiation sources by 
the two methods and an analysis of the factors 
affecting the accuracy seems desirable. Cebalt-60 
and cesium- 137 gamma-ray sources have already 
been measured in this laboratory by cavity-chamber 


—_—____., 


' Figures in brackets indicate the literature references at the end of this paper. 
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methods and an analysis of the accuracy has been 
performed [5]. Similar treatment is needed for the 
free-air chamber method and this paper describes 
such an analysis. 

One of the factors which may limit the accuracy of 
the free-air chamber measurements is the determi- 
nation of the saturation current. In order to limit 
the dimensions of the free-air chamber, the ionized 
air is often at high pressure. Under this condition 
it is impossible to collect all of the ionization; one 
must rely on a theoretical extrapolation to an infinite 
value of the collecting field. 

One method of performing this extrapolation, the 
Jaffé-Zanstra theory [7], has been checked and 
verified experimentally by Taylor et al. [8] for air 
pressures of up to 10 atm.’ According to this 
theory [7], a plot of the inverse of the collected 
ionization current versus a particular function of 2, 


f(x), should give a straight line and the intercept 


at f(x) =0, where the field strength is infinite, should 
give the inverse of the saturation current. Here 
r=c(} ’/DP)’; where V=collecting potential in volts; 
D=plate seperation in em; P=pressure in atmos- 
pheres; and ¢ is empirically determined from the 
requirement that the inverse of the collected ionization 
versus f(z) should be a straight line. Curves of f(x) 
versus x are available in the literature [7]. However, 
at least two very different values of ¢ have been pub- 
lished [8, 9]. The latter reference also recommends 
a field strength per atmosphere of at least 80 v em™! 
atm, contrary to the earlier work [8] where fields 
of less than 10 v cm™ atm give data that fall on 
the straight line plot. 

A second method of extrapolation to an infinite 
field has been suggested by Kara-Michaelova and 
Lea [10]. They plot the inverse of the measured 
ionization current versus the inverse of the collecting 
potential. For sufficiently high potential, a linear 
relation exists which may be extrapolated to an 


2 The authors neglected the correction for the scattered photon contribution. 
When such a correction is made for cesium-137 gamma rays, the extrapolated 
current per atmosphere at 10 atm is at least 2 percent less than that at 1 atm. 
The correction should be larger for the 390-kv X-rays used for this test by Taylor 

et al. 
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infinite collection field. The authors also pointed 
out the equivalence of the two methods for obtaining 
saturation current when the collecting potentials are 


sufficiently high. Therefore, an exploration of the 
feasibility of the suggested extrapolation methods 
seems desirable for the present application. At the 
same time, it is worthwhile to obtain the radial dis- 
tribution of the ionization from the primary electrons 
and from the scattered photons and the distance 
required for electron equilibrium. These are neces- 
sary [11] design data for the construction of a free- 
air chamber for measuring the gamma rays from 
cesium-137 and cobalt-60 in roentgens. 

The main objective of this report is, therefore, to 
carefully measure a cobalt-60 and a_ cesium-137 
gamma-ray source with a free-air chamber, and com- 
pare these measurements with those obtained with 
the cavity chamber. With a careful analysis of the 
accuracy of each type of measurement it may be 
possible to indicate which is capable of greatest 
accuracy at the present time and therefore, can be 
designated as the standard instrument in the million- 
volt region. However, the free-air chamber measure- 
ments require an exploration of the feasibility of the 
extrapolation methods for saturation and the de- 
termination of the design data necessary for the con- 
struction of free-air chambers for measuring cesium- 
137 and cobalt-60 gamma rays in roentgens. 


s 
< 


2. Experimental Arrangement 


A small beam of gamma rays from a radiation 
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diaphragm (fig. 1). This aperture is at the end of a 
cylindrical hole parallel to, but displaced from the 
axis of the cylindrical shutter-diaphragm. Rotation 
of the shutter cuts off the beam and permits measure- 
ment of the leakage radiation without introducing 
a modification of the scattering conditions. The 
beam passes first through a thin aluminum window 
(0.8 mm thick), then midway between and parallel 
to the two vertical plates of a parallel-plate ionization 
chamber, and then through a glass exit window. The 
ionization chamber is sealed inside a large tank filled 
with air under high pressure so that the mass of air 
between the plates is adequate for totally absorbing 
all the energy of the highest energy electrons pro- 
duced by the beam of gamma rays. 

One of the plates of the ionization chamber is 
maintained at a high potential with respect to ground 
(see sec. A—A of fig. 1); the other is divided into 3 
collecting regions (plates marked 2, 4, 5a, 5b, 6a, 6b, 
7a, 7b, and 9) and 4 large and 6 small grounded guard 
plates (see elevation view of fig. 1). The center 
collecting region is divided into 7 separate insulated 
plates (marked 4, 5a, 5b, 6a, 6b, 7a, and 7b). Each 
collector plate is separated from the adjacent plate 
by an air gap of 0.6 mm. A flat grounded plate is 
installed between the collector plates and the tank 
at about 3 em from the collector plate to reduce the 
possibility of ionization collection from behind the 
collector plate. 

The ionization to any of the collecting plates can 
be measured. During an exposure the plates which 




















































































































source is defined by an accurately known aperture | are not used for ionization measurement are 
(approximately 2 cm in diameter) in the shutter- | grounded. After careful alinement of the plate 
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The experimental arrangement. 





system, the displacement of the plane of each collector 
from the adjacent guard plate is measured. Gener- 
ally, this is less than 0.03 mm but increases to as 
much as 0.1 mm very near the top of plate 2 and the 
bottom of plate 9 because of warpage of these plates. 
This causes a distortion of the field. However, 
previous measurements [11] indicate that as the 
ionization collected in this latter region is small, this 
distortion should cause no more than 0.1 percent error 
in the collected ionization. 

This collector plate arrangement provides a means 
for determining the adequacy of the plate separation 
at a given pressure and also provides a means for 
determining the air absorption of the gamma-ray 
beam. With proper guarding of the electric field, 
each plate collects ionization from a volume defined 
by the plate height, the length of 20.15 em (fig. 1), 
and the distance between the collector and high- 
voltage plates. Thus, by observing the relative 
amount of ionization collected by each of the plates 
4 to 7, it is possible to determine if the effective plate 
separation (the product of plate separation and 
pressure) is adequate to collect all the ionization 
from electrons produced in the beam. 

The ionization current collected by plate 9 is 
smaller than that collected by plate 2, due to air 
absorption. Therefore, the ratio of these two meas- 
urements may be used to compute the air absorption 
coefficient. Comparison of these computed values 
at the different pressures for each source provides 
an independent estimate of the precision of the 
current measurements. A comparison with the 
theoretical absorption coefficient indicates whether 
there is a difference in any systematic error between 
the current measurements to the two plates. 

Guard bars are used to assure that the grounded 
tank does not distort the field in the collecting region. 
These are T cross section aluminum bars equally 
spaced between the chamber plates and_ parallel 
to these. In order to avoid extra photon attenuation 
and scattering, vertical bars are not used in either 
the front or back of the guard system where they 
might be struck by the beam (fig. 2). In these 
regions, fine aluminum wires are used instead of bars. 
A resistance network with a tap for each bar divides 
the potential uniformly between the high voltage 
and ground. A determination of the ability of these 
guard bars to eliminate the effect of the pressure 
tank on the electric field between a collector and the 
high-voltage plate requires an auxiliary experiment. 
As pointed out previously [5, 11], this is accomplished 
by noting the variation in collected ionization current 
as the potential of the tank is varied. lonization 
measurements to each of the collectors indicate no 
variation greater than +0.1 percent as the tank 
potential varies from ground potential up to the 
potential of the high-voltage plate. Thus, the tank 
produces no distortion of the collecting field which 
changes the collecting volume by more than the 
precision of measurement (about +0.1 percent). 
All subsequent measurements are performed with 
the tank at ground potential. 
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FiGguRE 2. A view of the back of the tank with the end beil 
removed showing the guard bars for the chamber. 


The air pressure in the tank is determined from the 
reading of a bourdon-type gage and the barometric 
pressure. The gage is calibrated by a piston gage 
and shows a short-time (over about 2 yr) maximum 
drift of calibration of about +0.2 psi and a long- 
time stability (over 8 yr) of about +0.3 psi except 
for a progressive zero shift (0.4 psi over 8 yr). For 
pressures of less than 5 atm, a mercury manometer 
also is used. In this range, the two methods agree 
with each other to within 0.2 percent. Filters and 
an ice-cooled condenser in the pumping system, as 
well as silica gel on the inside of the tank, maintain 
a low water vapor pressure. Mass spectrometric 
tests of the air under pressure indicate that the 
composition is the same as room air. 

Separate thermocouples are used to test for possible 
temperature differences of the high-pressure air in 
the various ion collecting regions of the tank, and to 
calibrate as a function of pressure a mercury ther- 
mometer placed inside the tank near the exit window. 
For this test, thermocouples are placed at the center 
of each of the collecting regions as well as near the 
thermometer. After a short stabilizing period, all 
of the thermocouples give temperatures agreeing 
to better than 0.2° C. Moreover, the calibration 
of the thermometer at a given pressure appears to 
be constant to 0.2° C, although its calibration varies 
by nearly 1° C over the pressure range used. 





In 
scattered photons, the ionization current also is 
obtained with a 10-cm diam 2-m long, Bakelite tube 
surrounding the collimated beam of photons inside 


order to determine the ionization due to 


of the tank. This tube is of sufficient thickness 
(about 3 mm) so that it stops all electrons originating 
in the air within the tube, and of sufficient diameter 
so that it is not struck by the primary beam of gamma 
rays. 

To reduce the possible distortion of the electric 
field by the tube, colloidal graphite is painted in 
strips along the length of the tube. Adjacent strips 
are insulated from each other by uncoated tubing. 
Each of these strips is maintained at a potential 
which minimizes the electric field distortion by con- 
necting it to the proper guard bars. Observation of the 
ionization from background radiation * with and 
without the tube in place gives an estimate of the 
field distortion produced by the tube. With the 
tube in place the current should be the same as those 
without the tube for plates 5a, 5b, 6a, 6b, 7a, and 7b. 
However, with the tube in place plate 4 will not 
collect ionization from the area in the tube shadow, 
that is, between the tube and the high voltage plate 
and within the tube but the difference can be com- 
puted from the geometry. The measured current 
to the different collectors are in the same ratio to 
within 5 percent as those predicted from the geom- 
etry. As the scattered photon contribution is only 
2 or 3 percent of the total ionization, field distortion 
caused by the presence of the tube does not affect 
the final result. The ionization produced by elec- 
trons generated by the primary photons is, within 
the precision of measurement, proportional to the 
difference between the ionization measured with * 
and without this tube. 

The ionization-current measuring system consists 
of a null indicating electrometer with feedback, a 
carefully calibrated condenser, and a_ variable 
bucking potential provided by a _ potentiometer. 
Timing of each exposure is accomplished by means of 
a scaler and a 1-ke oscillator {12].. This chrono- 
graph gives timing accuracies of about 0.01 percent 
for the times used here. 

Measurements are made with the shutter closed 
in order to correct for the background. As the 
shutter-closed measurement is not changed when a 
solid lead plug is inserted into the hole in the lead- 
shot-filled shield, the leakage radiation through the 
shutter is negligible and the shutter-closed measure- 
ment indicates the background. 

In addition to corrections for leakage radiation and 
air-scattered photons, it is also necessary to correct 
for any ionization from photons scattered by the 
source shield or diaphragm. Consider any rays 
which are scattered by the source shield and then 
pass through the diaphragm aperture. It may be 
shown for the geometry of figure 1 that the path of 


3 The background radiation is assumed to produce a uniform density of ioniza- 
tion over the whole collecting region. 

‘ A slignt correction is required for the shadow cast by the tube and a larger 
correction for lack ofsaturation. The reduction of the scattered photon ionization 
by photon absorption in the tubewallis much smaller than 0.1 percent of the tube- 
out ionization. Therefore this reduction can be neglected. 
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such rays will pass through the collector regions at 
distances greater than 5 cm from the axis of the 
primary beam. Thus, the ionization that they 
produce in the collecting region will be measured 
both with and without the Bakelite tubing in place, 
and the difference in readings will not include ioniza- 
tion either from photons scattered by the source 
shield or from the air in the chamber. 

The effect of scattering from the diaphragm is 
explored by comparing the ionization current with 
diaphragms 15 cm and 22.5 em long. 

The source holders are made as light as is cons’st- 
ent with dimensional stability so as to reduce the 
amount of scattered radiation from them. These 
source holders fit into the large hole in the rear end 
of the source shield. The position of the source is 
determined by observing the position of the end of 
the long rod attached to each source holder. The 
radiation is measured perpendicular to the axis of 
the 3.2 mm diam by 3.2 mm long cobalt-60 source 
and along the axis of the 4.6 mm diam by 5 mm long 
cesium-137 source. The cobalt-60 source is en- 
capsulated in aluminum with a 1.6 mm thick wall 
and the cesium-137 in 1 mm thick stainless steel. 

A separate piece of equipment similar to that 
previously used [13] provides a means of determining 
the distance required for electronic equilibrium [11]. 
A grid chamber (fig. 3), made of colloidal-graphite- 
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Figure 3, Schematic plan view of the thread chamber and 


electron diaphragms. 


For simplicity the threads in each plane are drawn as a single line. 





coated nylon threads (~ 0.2 mm diam) in 5 parallel 
planes spaced about 1 cm apart, is placed in the 
pressure tank with the planes perpendicular to the 
gamma-ray beam. The threads in any one plane 
also are spaced about 1 cm apart. The threads in 
the central plane form a collector at near ground 
potential, and the two adjacent planes are connected 
to the collecting potential. The collector plane has 
two collecting regions separated by uncoated thread 
for insulation. A central collecting region, which is 
coaxial with the gamma-ray beam, is about 6.5 cm 
diam (AB in fig. 3). The outer collecting region has 
an inner diameter of 8.5 em (CE) and an outer diam- 
eter of 51 em (DF). The gamma-ray beam diam- 
eter at the position of the collector plane is about 
6.5 em. Thus, the central collector measures the 
ionization in a cylindrical region of about 2 em 
depth and adiameter of 7.5 cm, and the outer collector 
measures the ionization in a ring of about 2 em depth, 
7.5 cm inner diameter and about 51 cm _ outer 
diameter. 

A circular electrically-conducting Bakelite sheet 
(electron diaphragm), which is thicker than the 
range of the electrons produced by the gamma rays 
in air, is positioned perpendicular to the gamma-ray 
beam and in front of the grid chamber. A similar 
sheet is positioned behind the chamber. Each sheet 
has a center hole, which is slightly larger than the 
gamma-ray beam and coaxial with it, and each can 
be moved along the beam by remote control from 
outside the pressure tank. 


3. Results and Discussion 


The ionization current to each of the collector 
plates is measured both with and without the Bake- 
lite tube and for a variety of pressures and collecting 
potentials. Data are obtained for five different 
pressures from 3.96 to 10.75 atm with the gamma 
rays from cesium-137 and for pressures of 10.75 and 
12.27 atm for cobalt-60 gamma rays in order to 
choose a value for the Jaffé-Zanstra constant, c. 
With the Bakelite tube in place part of the ionization 
from the scattered photons is produced in the shadow 
of the tube and is thus not collected on plates 2, 4, 
and 9. However, the ionization produced in a ring 
about the axis of the beam decreases only slowly 
with increase of the ring radius®. Therefore, it is 
relatively simple [14] to make a correction for this 
shadowing, and also to change from the rectangular 
coordinate system of measurement for plates 4 to 7 
to a cylindrical set of coordinates. For a given 
pressure and collecting potential, the difference be- 
tween the ionization measured without the tube in 
place and the shadow-corrected determination with 
the tube in place is proportional to the ionization 
produced in the collecting length (20.15 em) by the 
electrons generated by the primary gamma rays. 





* Scattering in the source shield or collimator might distort this monotonic 
decrease but such effects prove to be small in the present arrangement. Correec- 
tion is still required for lack of saturation. 


If the pressure is too low, so that some of the 
energy of the electrons is dissipated in the plates, a 
linear extrapolation to larger radii may be per- 
formed [14] on a log-current-loss versus radius plot. 
For the data reported here, this extrapolation does 
not exceed % percent. 

Data from the large range of pressures used with 
cesium-137 gamma rays permits a study of the 
techniques for extrapolating to infinite collecting 
fields. Figure 4 shows plots of the inverse of the 
ionization currents,’ J, against the Jaffé-Zanstra 
parameter, f(z). Four values of ¢ are used to obtain 
f(x) in order to determine the value of ¢ which 
gives a straight line. The value of 1.24107 is 
that used by Taylor et al., [3] and by Grove [4]. 
It is seen that the point for the largest value of 
f(x) falls well below the straight line through the 
remaining points. However, for both e=3.2107° 
and c=2.5 1075, all data fall on straight lines with 
slightly different intercepts at f(x)=0 (infinite 
collecting field). For c=1075, the points for large 
values of f(x) fall above the straight line for the 
remaining points. The proper intercept on the 1/J 
axis therefore probably lies between that obtained 
for c=10~ and for c=1.2410~*. 
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1 2 3 4 f(x) 


FicurE 4. Sample plots of the extrapolation procedure with 
cesium-137 gamma rays at a tank pressure of 3.96 atm. 
Four Jaffé-Zanstra type plots are indicated using different valves of c. A 


Kara-Michaelova and Lea type plot (marked J0/kv) does not give a straight 
line relation. 


Also included in figure 4 is a plot (marked 10/kv) 
of the inverse of the ionization current, J, against 
the inverse of the collecting potential (Kara-Michael- 
ova and Lea). The intercept on the 1/J axis for 
the inverse collecting potential plot is indeterminate 
as the points do not fall on a straight line even 
for the three smallest values of 1//, where the 
effective electric field strengths are about 171, 114, 
and 85.5 v em! atm™?. 

While the uncertainty in the intercept by these 
two methods is small at this pressure because the 
extrapolation is only about 1% percent, the uncer- 
tainty becomes larger at the higher pressures where 
the extrapolation is as large as 10 percent. For- 
tunately, there is another criterion for choosing 
the proper extrapolation procedure. 








6 Tube-in minus tube-out determination after correction for tube shadowing. 
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The exposure dose rate at the aluminum window 


is independent of the pressure in the tank. The 
exposure dose rate is proportional to the ionization 
current per atmosphere extrapolated to an infinite 
field and corrected for air attenuation, i.e., to 
T.e’*/P. Here J, is the current at infinite field, 
P is the pressure in atmospheres, d is the distance 
between the center of the collector plate and the 
entrance window, and yw is the linear absorption 
coefficient of the gamma rays inair. Valuesof /,,e°%/P 
at different pressures are listed in table 1 for cesium- 
137 gamma rays. These results verify the earlier 
conclusion that a proper extrapolation is obtained 
for a value of ¢ between 1.24107' and 107°, As 
the extrapolated values for e=2.5*107-> show a 
slightly smaller deviation from the average than 


c=3.2107°, this value is used for all subsequent 
determinations for both the cesium-137 and the 
cobalt-60 source. 
TABLE a 10! P eP du for cesium-137 
Pressure 1.24x 10-4 3.2x10 2.5\10 I 
m 

10.75 1. 956 1. 992 2.005 2 O54 
8.45 1. 977 1. 997 2. 009 2. O45 
6.70 1. 990 2.013 2.018 2.041 
+.79 2. 001 2.015 2. 018 2. 030 
3.06) 1. 99S 2. 02 2. 004 2.01 

Average 2 OO4+0.008 &  2.011-+0.006 @ 

t Mean deviation 


The computations of e? assume that the photons 
which generate the electrons are attenuated for a 
distance, d, in air. Actually, most of the electrons 
which produce the ionization in the collecting region 
are generated before and a few after the photons 
reach the collecting region. The amount by which 
d must be reduced because of this is the distance, 
a, from the center of the collecting region to the 
mean position of origin of the electrons producing 
ionizations in the collecting region, <As the value 
of a is nearly a constant fraction of the distance for 
electronic equilibrium for the photon energies con- 
sidered here, it is convenient to first determine the 
distance’ for electronic equilibrium. These dis- 
tances for the two sources may be obtained from 
the thread chamber data. 

Results obtained with the thread chamber for 
cesium-137 gamma rays at 6.7 atm are shown in 
figure 5 and for cobalt-60 gamma rays at 12 atm 
in figure 6. The solid curves give the ionization 
current to the outer collector ring for different 
positions of each Bakelite sheet with the other sheet 
positioned as close as possible to the collecting 
region; the dotted curves are for the ionization 
current to the central collector. Initially the solid 
curves rise rapidly with the distance between the 
sheet and the collecting region, but eventually reach 
a plateau. The rise indicates that electrons generated 
at larger distances from the collecting plane are 
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Figure 5. Data obtained from the thread chamber and cesium- 


137 gamma rays at 6.7 atm, 
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dotted curve to the 


e the ionization current to the outer collecting ring and the 
center electrode 
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FIGURE 6. Data obtained from the thread chamber and cobalt-60 
gamma rays at 12 atm, 


rhe solid curves give the ionization current to the outer collecting ring and the 
lotted curve to the c 


nter electrode 

contributing to the ionization. The distance at 
which the ionization becomes constant is the distance 
for electronic equilibrium. From the solid curves it 
appears that the distance required for electronic 
equilibrium with cobalt-60 gamma rays is about 350 
em-atm in front of the collector and about 180 
em-atm behind the collector. The corresponding 
distances for cesium-137 are 90 and 50 cem-atm. 
Aglintsev and Ostromukhova [15] also obtained a 
distance of about 350 em-atm for cobalt-60 gamma 
rays. 

The mean distance from the collecting region to 
the mean position of origin of electrons producing 
ionization in the collecting region is, of course, less 
than the equilibrium distances. Calculations by 


Roesch [16] indicate the distance from the collecting 
region to the mean position of origin of the electrons 
is about 0.3 times the distance for electronic equilib- 
Aglintsev and Ostromukhova [15] computed 


rium. 





alues of 0.36 for cobalt-60 and 0.34 for cesium-137 
gamma rays. For these distances in a 
chamber they compute the transmission, K= 
of cesium-137 and cobalt-60 gamma rays to be 0 996 
and 0.992, respectively. In a carbon- walled cavity 
chamber, because of the waka effective absorption 
coefficient there, these transmissions are 0.998 and 
0.995. These factors must be multiplied by e? 
in order to obtain the attenuation of the photons 
whose electrons produce the measured ionization. 

The attenuation cofficients for the gamma-ray 

beams are computed from data obtained for plates 
2 and 9. The difference in the currents obtained 
with the tube out and with the tube in, after cor- 
rection for the shadowing and saturation, gives the 
current produced by the primary beam. The ratio 
of the plate-2-to-plate-9 currents are given in table 2 
for different pressures with the two sources. The 
linear absorption coefficient computed from the ratio 
and that given by Grodstein 117] are also shown. 
The maximum deviation of the computed absorption 
coefficient from the mean is approximately 6 percent 
and the maximum difference between the experi- 
mental and theoretical values is also 6 percent. The 
first number indicates that the precision of the 
current ratio is at least 0.5 percent and the second 
that the systematic errors if they exist for two cur- 
rent determinations are equal within the precision 
of measurement. As there is no reason to assign 
different prec isions to the two current determinations, 
the precision of each is about 0.35 percent. 


Pen 


TABLE 2. Air absorption coefficient 


uecTD 
source Pressure Te 
Experiment Grodstein [17] 
atm 
10. 75 1.094 | 1.04X10-4 1.0K10-4 
Cc S. 457 | 1.068 | 0. 9610-4 
, | 6.70 1. 057 1.04% 10-4 
1 | 4.797 1. 042 1. 06X 10-4 
12. 27 1. O76 7.5X10-5 | ican ’ 
wo 10.75 1. 069 7.7X10-3 1. 3X10 


Possible scattering from the diaphragm is investi- 
gated by obtaining data with both a 15 em and a 
22's em long diaphragm. The radiation scattered 
through the diaphragm and into the chamber from 
the source holder is decreased slightly by lengthening 
the diaphragm. However, this decrease is estim: ated 
to be less than 0.1 percent of the electron contribu- 
tion to the ionization. Such a small decrease is not 
observable from the tube-out measurement as this is 
than the precision. Large angle scattering 
from the shutter-diaphragm is considerably attenu- 
ated by the steel surrounding the entrance window 
of the pressure tank. Probably for this reason no 
difference is obtained in the current to plate 7a and 
7b when the diaphragm length is increased. Cur- 
rents to plates 5a, 5b, 6a, and 6b are decreased 


less 


, (3,690) (2.998) (16°) 7 b2A e? He? Al e! 
PADL 


93 


ip e-air | 


slightly (0.1% of the electron contribution) as the 
diaphragm length increased for the tube-out 
measurement. There is a similar decrease for the 
tube-in measurement. Therefore, the differences in 
the tube-out and tube-in measurements are inde- 
pendent of diaphragm scattering. 

From the measured currents to plates 4, 5, 6, and 
7, it is possible [14] to compute the loss in ionization 
(electron loss) obtained beyond different radii and 
pn the contribution to the ionization provided by 
the scattered photons. In order to compare the 
results obtained at different pressures for each 
source, the radii at each pressure are corrected for 
the diameter of the gamma-ray beam [13] and for 
the density of the air. The abscissa of figures 7 and 
8 is in equivalent centimeters at 1 atm for a zero 
diameter beam. It is seen that the data of figure 7 
agree to within about +0.2 percent for different 
pressures and that the data of figure 8 agree to 
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within about +0.1 percent for different pressures. 
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FiGguRE 7. The loss of ionization from primary photons be- 


yond different radii for cobalt-60 and cesium-137 gamma rays. 


Results have » | een corre 
with cesium-137 at 3.96 atm 
it 10.75 atm= 


‘ted to 1 atm and for the 
, at 4.797 atm=Y; 


beam diameter. Results 
with cobalt-60 at 12.27 atm=O, 


The gamma ray output, G, in roentgens per hour at 
1 m (rhm) for each source is obtained by substituting 
numerical values in eq (1). 


+0 40) pu / 0.001293 


rhm 
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FIGURE 8. 


Contribution of secondary photon 
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ionization in 


percent of total primary electron ionization for different radii 
from the cesium-137 and cobalt-60 beam. 


Results have heen corrected to 1 atm and for the heam diameter. 
with cesium-137 at 10.75 atm : 


Results 


>), at 8.457 atm =O, at 6.70 atm =A; with cohalt-00 


at 12.27 atm=0, at 10.75 atm=Y. 
where 
3,600=see hr, 
(2.998) (10°)=esu of charge coulomb, 

/ ionization current in amperes pro- 
duced by the primary beam in 
the collecting region, 

b=0.871 m (see fig. 1), 

K=0.992 for cobalt-60 and 0.996 for 
cesium-137, 

P=air pressure in atmospheres, 

d=distance from collecting region to 
entrance window, 133.9 em, 

u=air absorption coefficient in em! 
atm~!, 

u4;—aluminum absorption coefficient in 


¢ 100 (6 +0. 40) pu/0.00129 


p 
0.001293 


— 


em", 

air absorption correction between 
source and window, 

density of laboratory air in g em~$, 

density of air at a pressure of 760 
mm mercury and a temperature 
of 0° C in g em, 

measured area of aperture of dia- 
phragm, 3.147, em, 

effective area of aperture 

» and 


collector length (20.15 em). 


The effective area of the aperture is greater than the 
area A because of photon leakage through the edges 


of the aperture [3]. 


for this leakage. 


The computed factor, D, allows 


Table 3 gives the rhm for each of the sources com- 
puted from the equation with the data obtained with 
the pressure chamber. 


TABLE 3. Output determinations 


Source Free-air chamber Cavity chamber Difference 


{ (1.003) 0.758 rhm 0.3 

Cr 0.756 rhm (1.016) ® 0.768 1.6 
| (1.003) ® 0.758 0.3 

| { (1.007) 2.5 | 1.6 

Cs | 2.61 14 (1.026) ® 2.61; 0.3 
| (1.013) 2.586 ___- 1.0 


“The numbers in parentheses are the ratio of the stopping power of carbon to 
that of air for **‘ Bakker-Segre,” *‘ Caldwell,” and “NCRP” mean excitation po- 
tentials, respe tively, for each of the sources. 


The carbon cavity chamber previously used by 
Attix and Ritz [4] is also used to measure the sources. 
The largest single uncertainty in such a determina- 
tion is in the stopping power. Attix and Ritz have 
indicated two sets of values for the mean excitation 
potential to be used in the stopping power equation. 
For one set, labeled “Bakker and Segre,”’ they give 
76.4 and 80.5 ev for carbon and air respectively; 
for the other set, labeled “Caldwell,” they give 78 


and 94.9. An analysis by a task group of the 
National Committee on Radiation Protection and 


Measurements [18] of the information presently 
available provides values of 78.4 and 85 ev for carbon 
and air respectively. The resulting stopping power 
ratios (in parentheses) obtained after taking into 
account the density correction and the outputs de- 
termined from the cavity chamber data for each are 
listed in table 3. 

The emission constant of the cobalt-60 source may 
also be computed. The activity in curies is obtained 
from the measured value of Minton [19] with a radio- 
active decay correction. Corrections are also applied 
to the values in table 2 for photon absorption in the 
aluminum capsule (1.011) and in the cobalt-60 
evlinder (1.030). The resultant emission constant 
is 12.7, rem? me"! h-!.) The theoretical value for 
this constant is 12.9 if the average energy to produce 
an ion pair is 34.0 ev. In view of the uncertainty in 
the activity determination (196) and the uncertainty 
in the decay correction for the four vears between the 
activity and ionization-chamber determinations, the 
agreement seems reasonable, 


4. Accuracy 


In order to assign an accuracy to the free-air and 
cavity-chamber values in table 3, it is necessary to 
review the accuracy of each of the factors used in 
determining the values. These factors are listed in 
table 4. 

The error for item 1 is estimated from the repro- 
ducibility of individual measurements. Part of this 


fluctuation is caused by slight drifts in the collecting 
and bucking potentials and possibly from charging 
of some of the insulators. 





TaBLE 4. Possible errors in rhm determinations 


Percent error 
Item Errors due to » 
Pressure 
chamber 


Cavity 
chamber 


1 Ionization measurements ® 2 a().2 (1 
2 | Time, capacitance, and potential a 3 
3 | Voltage saturation. : b3 i 
4 | Air absorption a l 
5 | Window absorption soa om a 
6 | Plate separation correction__-_- ae ‘ 
7 | Air pressure __ on Soil 
& | Temperature . ae on 
9 | Effective collector length ¢ 2 
10 | Diaphragm area of 
11 | Diaphragm leakage_- . —_ 
12 | Chamber volume > l 
13 | Diaphragm-to-source distance E <4 2 
14 | Scattered radiation and chamber wall alter- 
nate.... - : .3 
15 | Seattered photon contribution od 
16 | Electron contribution 2 
7 1) «. —— .2 
18 | Stopping power ratio-_ -- -5 


Of individual points. 
Including relative measurement of currents for different pressures but not 
correction for error for item 2. 
¢ Including field distortion. 


The timing accuracy depends principally on the 
difference in action time of two relays [12]. An 
investigation indicates that the accuracy is about 
().01 percent for the times used here. The condenser 
and potentiometer are each calibrated to an accuracy 
of 0.05 percent but the capacitance drifted during 
the experiment by almost 0.1 percent. However, 
it does not seem probable that the determination of 
the charge collected on the condenser is in error by 
more than 0.1 percent (item 2). 

The accuracy of the saturation extrapolation 
depends on the adequacy of the Jaffé-Zanstra theory. 
Each extrapolation is obtained by a least squares fit 
so there is no subjective error. The amount of extra- 
polation at 3.96 atm for cesium-137 is only about 
1's percent. However, the constancy of (/.e?%/P 
in table 1 for cesium-137 and similar values for 
cobalt-60 (which agree to within £0.15%) indicates 
that the error in extrapolation is not likely to be over 
0.3 percent (item 3). For the cavity chamber the 
extrapolation is only 0.1 percent so the uncertainty 
must be small compared to 0.1 percent. 

The air attenuation between plates 4 to 7 and the 
source varies from about 7 percent to 17 percent for 
cesium-137 and is about 12 percent for cobalt-60. 
The distance used in the computation is known to 
at least 0.1 percent. There appears to be no accurate 
experimental determination of the absorption co- 
efficient for air although the values for aluminum 
agree with theory to within a few tenths of a percent 
for the sources considered here. Therefore, it seems 
reasonable not to expect an error of more than 
0.1 percent for item 4 because the maximum correc- 
tion for air absorption is only about 17 percent for 
the pressure-chamber measurement. For the cavity- 
chamber measurements, the maximum correction is 
only about 0.7 percent so the uncertainty in this 
item should be negligible. 
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The absorption in the aluminum window of the 
pressure chamber is of the order of 1 percent. It is 
difficult to imagine that variations in thickness, 
absorption coefficient, or composition could cause 
uncertainties of more than 0.1 percent (item 5). 

Several extrapolations for inadequate plate separa- 
tion were tried. For any reasonable assumption on 
this extrapolation, the uncertainty is probably not 
more than 0.1 percent (item 6). The small variation 
with pressure indicated in figures 7 and 8 corroborates 
this. 

The density of air in the pressure chamber is de- 
termined from the gage pressure and temperature 
in the tank and the barometric pressure and room 
temperature. Errors in the room pressure and tem- 
perature could cause less than 0.1 percent error in 
the air density within the tank at the pressure used. 
However, as the gage could be calibrated to only 
+0.1 psi and it was not calibrated at the exact 
pressures used, there is a possible uncertainty of the 
pressure of about 0.1 percent at the higher pressures 
(item 7). The same uncertainty is estimated for 
measurements below 5 atm where the mercury 
manometer is used. Temperature fluctuation in the 
tank as indicated by the thermocouples could not 
have caused an uncertainty in the density determina- 
tion of more than 0.1 percent (item 8). 

For the cavity-chamber measurements, the air 
pressure is given by a barometer calibrated to 0.1 
mm and the pressure was very constant during a 
set of measurements. Thus, the error in item 7 
should be less than 0.1 percent. However, the 
temperature did vary by about +0.1 percent during 
a set of measurements so a possible error of 0.1 
percent is assumed for item 8. 

The effective length of the collector plate is deter- 
mined by the measured length and by any field 
distortion caused by the presence of the grounded 
tank and by warpage of the plates. The measured 
length could not be in error by as much as 0.1 percent 
but, as indicated earlier, the distortion produced 
by the grounded tank is determined by ionization 
means. Therefore, the reproducibility of the ioniza- 
tion measurements (0.1 % for this investigation) de- 
termines this possible error. From previous work 
[11] it appears that the warpage of the plate could 
produce an uncertainty of no more than 0.1 percent. 
Therefore, for these two possible errors it is estimated 
that the effective plate length could not be in error 
by more than 0.2 percent (item 9). 

The diaphragm area is determined by measuring 
four radii at different depths. These measurements 
show the hole to be slightly coned (about 3 parts 
in 7,800) at the two ends, but at any one depth the 
radii agree to 8 parts in 78,000. The uncertainty 
in the physical determination of the hole area is 
therefore not more than 0.1 percent (item 10). 

The effective area of the diaphragm is slightly 
larger than the hole area because of leakage of 
radiation through the edge of the holes [3]. The 
increased area amounts to about 1.8 percent for 
cesium-137 and about 3.5 percent for cobalt-60 
radiation. It is estimated that the uncertainty 
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[3] in this leakage could not cause an error of more 
than 0.2 percent (item 11). 

The uncertainty in the cavity chamber volume 
(item 12) is assumed to be the same as that estimated 
by Attix and Ritz [5] (0.1%). 

The distance between the collimator and the source 
(87.1 em) is determined by a cathetometer and 


vernier calipers reading to 0.1 em. It is estimated 
that the error introduced by the uncertainty of 
this distance is than 0.1 percent (item 13). 
However, the uncertainty of the source-to-cavity- 
chamber distance is somewhat larger—about 6.2 
percent. 

The accuracy of the correction for scattering from 
the walls of the room and from the cavity-chamber 
walls, and for the photon attenuation in the wall is 
estimated to be the same as that given by Attix and 
Ritz [5] (item 14). 

The computed scattered photon contribution to 
the ionization for cesium-137 gamma rays amounts 
to about 4 percent of the electron contribution at 
10.85 atm and about 2 percent at 3.96 atm; for 
cobalt-60 this contribution amounts to about 3!5 
percent for the pressures used. The agreement (see 
fig. 8) in the scattered photon contribution obtained 
at various pressures for a given radiation indicates 
that it is reproducible to about 0.1 percent (item 15). 

The agreement (see fig. 7) in the electron contri- 
bution obtained at different pressures indicates that 
this quantity is reproducible to about 0.2. percent 
(item 16). 

The theoretical values of A determined for cobalt- 
60 gamma rays by Roesch |16] are about 0.3 percent 
larger than those determined by Aglintsev and 
Ostromukhova [15]. Values obtained by Aglintsev 
and Ostromukhova are used here, but the uncer- 
tainty of this correction (item 17) has been assigned 
a value of 0.2 percent. 
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Probably the largest single uncertainty in the out- 
put determination by the cavity-chamber is in the 
stopping power ratio (item 18). Three different 
values for carbon-walled air ionization chambers 
have been suggested recently. The maximum differ- 
ence between them is about 2 percent. However. 
the values labeled “Caldwell” in table 3 do not fit 
the data of Attix and Ritz [5] and of Whyte, who 
conducted experiments with cavity chambers having 
walls of different atomic number. The experimental 
data from these authors fit the “Bakker and Segre” 
and “NCRP” stopping-power ratios slightly better 
than “Caldwell” but neither obtained a 
direct measurement of the stopping power for air. 
The data of table 3 also indicate better agreement 
between the free-air chamber and cavity chamber 
for the “Bakker and Segre” and “NCRP” stopping- 
power ratios with the cobalt-60 source but poorer 
with the cesium-137 source. <Attix and Ritz esti- 
mate the accuracy of the stopping-power ratios as 
0.3 percent, but Whyte gives a value of 0.5 percent. 
The latter hes been adopted for the present report, 


those of 


The maximum possible error is obtained by adding 
the errors of item 2 through 18, and is 1.9 percent 
for the pressure chamber and 1.5 percent for the 





cavity measurements. Item 1 is greatly reduced 
by the multiplicity of observations. If one assumes 
that all of these items give random errors, the square 
root of the sum of their squares gives 0.6 percent 
for the pressure-chamber method and about 0.7 
percent for the cavity -chamber method. 


5. Summary 
1. By a proper choice of the constant in the 
Jaifé-Zanstra method of extrapolating the ionization 
current to an infinite field strength, it has been 
shown that consistent saturation currents can be 
computed from ionization measurements made with 
noninfinite fields, at least for pressures of 4 to 12 
atm in air. However, the best constant seems to 
be somewhat different than that previously used. 
It has also been shown that field strengths per 
atmosphere of greater than 85 and probably greater 
than 114 vem7! atm! are needed for the Kara- 
Michaelova and Lea extrapolation technique. 

2. Design data for the construction of free-air 
chambers used in the measurement of cesium-137 
and cobalt-60 gamma rays in roentgens are presented, 

3. Determination of the gamma-ray output from 
a cesium-137 and a cobalt-60 source by means of a 
free-air chamber and a cavity chamber are indi- 
cated, The experiments are not sufficiently accurate 
to unambiguously decide which value of stopping- 
power ratios should be used. 

4. It would seem from the estimate of errors that 
the uncertainty in the two methods of measurement 
is approximately the same and, therefore, that 
neither can be designated as the standard method, 


The author thanks T. P. Loftus, who made all 
of the cavity chamber measurements; F. H. Attix, 
with whom many of the problems arising during 
the investigation were discussed; and G. H. Demp- 
sev, who was of great assistance in setting up the 
equipment. 
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Emissivity in the Infrared 
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Apparatus and methods are described for measuring in the infrared the normal spectral 


emissivity of metals and coatings or oxides which tightly adhere to metals. 


Examples of 


the use of this apparatus are given in measurements of the emissivity of platinum and 


of oxidized Inconel within the spectral region of 1.5 to 15 microns. 


reproduced to better than 5 percent. 


1. Introduction 


The measurement of the spectral emissivity of 
heated substances may be accomplished by com- 
paring the energy emitted by the substance under 
study with that emitted by a standard substance 
for radiation at various wavelengths. For such a 
standard one usually uses a blackbody, that is, a 
hollow body provided with a hole for viewing the 
radiation inside. The inside surface of the black- 
body must be maintained at a uniform temperature. 
The wavelength separation may be accomplished 
using a double-beam monochromator or a single- 
beam monochromator. 

The specimen under investigation may be formed 
into a rod, bar, ribbon, or other shape. Some 
workers have even formed the specimen into a 
blackbody. For such a specimen both the inside 
surface (the blackbody) and the outside surface 
(the specimen) are observed. 

The double-beam monochromator and associated 
method have certain advantages especially with 
relation to the elimination of the effects of absorp- 
tion by carbon dioxide and water vapor, but the 
optical adjustment is critical as compared with a 
single-beam instrument. Furthermore through the 
use of a single set of optics and equal path lengths 
in the single-beam monochromator the effects of 
both atmospheric absorption and differential reflec- 
tivity in the auxiliary optics are practically elimi- 
nated. The apparatus described herein is of the 
latter type. 

The specimen may vary in temperature from 
point to point depending upon its shape, size, and 
method of heating. If made into the form of a 
blackbody the inside temperature may differ from 
that on the outside. A convenient shape is that of 
a long thin ribbon heated electrically by contacts 
at the two ends. The data presented in this paper 
apply to samples of this type for which the tem- 
perature is measured by means of thermocouples in 
metallic contact with the opposite face from the 
pomt of radiant energy Measurement. 








Measured values were 


2. Apparatus and Method 


This apparatus was set up as shown in figure 1. 
The blackbody and sample under investigation were 
mounted side by side on an optical bench (con- 
structed from a lathe bed and table) so they could 
be alternately placed at the focal point of the optical 
system. The auxiliary optics, consisting of a plane 
and a spherical mirror (65-cm radius of curvature), 
were mounted rigidly on the monochromator base. 
bes 
The placing of the sample and blackbody aperture 
alternately at the same position ensured essentially 
equal light paths and the use of identical opties. 


_— 


FicgurRE 1. Instrumental setup for the measurement of the 
normal spectral emissivity of heated substances 
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monochromator 


A single-beam 
sodium chloride prism and thermocouple was modi- 
fied to vield the arrangement shown in figure 1 for 


incorporating a 


use in this work. It was calibrated for wavelength 
in the usual manner by using atmospheric and poly- 
stvrene absorption bands. The light beam was 
chopped at 10 eps and the thermocouple voltage 
amplified through the use of a ten-cyele thermo- 
couple amplifier giving an approximately linear 
output. 

The linearity of the amplifier was of importance 
as the emissivity was obtained as the ratio of the 
measured outputs for the blackbody and specimen, 
both being maintained at the same temperature. 
For the high emittance measurements the linearity 
of the amplifier was checked using filters of known 
transmission. In the case of the low-emissivity 
measurements an attenuator switch was used to 
bring the blackbody and specimen readings within 
a comparable range. This attenuator was then 
calibrated by comparing the measured values of the 
radiant energy from a blackbody at a high tempera- 
ture with those from the blackbody at a lower 
temperature. Using the Planck blackbody radiation 
equation the ratio of the radiation from these two 
was calculated and used to calibrate the 
amplifier and attenuator. 

In this work the wavelength bandwidth (half in- 
tensity width) averaged about 0.15 » although it 
varied somewhat with the wavelength setting. 

The blackbody cavity was constructed of an alloy 
casting of essentially 80 percent nickel-20_ percent 
chromium (3 in. outside diam. by 6 in. length, 
with a wall thickness of about '4 1in.). The low re- 
flectivity of this (oxidized) metal coupled with the 
fact that the aperture was only 0.3 percent of the 
interior surface area resulted in a blackbody of ex- 
tremely high effective emissivity (probably 0.999 
or above). The high heat capacity of the furnace 
employed in heating the blackbody gave it a high 
thermal stability and made manual temperature reg- 
ulation easy (see fig. 2). The large size of the fur- 
nace, the method of heating, and the small furnace 
aperture used to permit viewing the blackbody 


sources 


vielded a very uniform temperature within the 
blackbody. 
T : | 
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Figure 2. Blackbody and furnace cross section. 





The temperature of the blackbody, as well as that 
of the sample under investigation, was measured 
through the use of platinum versus platinum—10 per- 
cent-rhodium thermocouples. Checks against stand- 
ard thermocouples indicated that the absolute tem- 
peratures of the blackbody recorded herein are cor- 
rect to within about +6° K. 

The specimen consisted of a metal strip about 8% 
in. long by 4 in. wide and about 4 in. thick. The 
thermocouple used to measure the specimen tem- 
perature was welded to the back side of it directly 
behind the surface area from which the measure- 
ments were made. The thermocouple used was of 
rather small diameter wire (5 mils) so that the cool- 
ing of the surface which resulted from conduction 
along the thermocouple wire was quite small. The 
useful re adiating surface area was about * in. by lig 
in. and the radiation used was essentially normal to 
the surface. Observations with an-optical pyrom- 
eter indicate that the temperature gradient over 
the observed surface was less than ~4° K. 

The specimen under examination was heated by 
passing through it an electric current of about 80 
amp at 5v. It was clamped rigidly at the two elec- 
trodes and held centrally located within a water- 
cooled metal tube, about 144 in. inner diameter. 
Although the inner surface of this tube was highly 
oxidized, reflection from its inner surface affected 
the emissivity measurements by about 1 percent. 
However, its use reduced air currents and added 
greatly to the stability of the temperature of the 
sample so that its temperature could be manually 
held constant within about £1.5° K even at 1,400° K. 


3. Examples of the Use of This Apparatus 


As examples of the use of this apparatus, the 
emittances of two very different substances, plati- 
num and oxidized Inconel, were measured. The 


oxide coating on Inconel has a very high spectral 
emittance while the surface of platinum being highly 
reflective has a very low emissivity. 


3.1. Oxidized Inconel 


The Inconel strips were prepared by first cleaning 
the surface by sandblasting, and then heating in air 
for 6 hr at a high temperature in order to form an 
opaque oxide coating. 

The spectral emittance of oxidized Inconel is given 
in table 1 and figure 3. These data are the average 
of the values obtained on three separate specimens 
each of which was run several times at each tem- 
perature. 

The oxide coating has a very high emissivity 
characteristic of many oxides. The most prominent 
features of the emissivity curve are two minimums 
at around 6 and 6.5 uw. In this same wavelength 
region there are two strong absorption maximums 
for atmospheric water vapor. Although care was 
exercised to make the pathlengths for blackbody and 
specimen equal, the atmospheric absorption may be 
the cause of these minimums. 
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TABLE 


Normal spectral emittance of oxidized Inconel at 


1,400°K, 1,050°K, and 780°K from 1.4 to 15 p 


different emissivity from a perfectly smooth surface, 
but infrared radiation, since it has relatively large 
wavelength, is rather insensitive to small surface 














The two platinum strips examined were first 
annealed by heating to about 1,525° K. One speci- 
men was held at that temperature for 1 hr while the 
other was heated for only 10 min. Both specimens, 
although initially smooth surfaced, acquired a crys- 
talline appearance after the annealing process. This 
crystalline surface would probably have a slightly 











beta 1,400°K | 1,050°K 780°K irregularities. 
The emissivity of the platinum strip is given in 
’ table 2 and the results are shown in figure 4. The 
-. “= = = low emissivity is characteristic of most metals. This 
2. 5 91 - 86 78 follows from the fact that the reflectivity of metals 
3.5 90 "85 78 is very high. 
0 - aa om Price [1]! measured the spectral emissivity of 
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asia - - * perfect blackbody which he used his emissivity 
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Ficgure 4. Normal spectral infrared emissivity of commercial 
grade platinum at 1,400° K. 


Observations by D. J. 


Price shown by triangles (A); calculated emissivity 
shown by squares (0). 


Figures in brackets indicate the literature references at the end of this paper. 
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In the long wavelength region the emissivity of 


metals increases with increase in temperature. This 
fact has been verified for platinum by Hagen and 
Rubens [2]. Thus our results should be higher than 
measurements made at lower temperatures. As 
expected the reflectivity values for platinum at 
lower temperatures reported by other workers vield 
emissivity values lower than ours. 

The measured spectral emissivity values for 
platinum may be compared with the theoretical 
emissivity calculated by using the Drude-Zener 
relation. This is valid only for metals for which the 
electrons are perfectly free. According to this 
relationship [3] the long wavelength spectral emis- 
sivity is given by 


1 
ope il P\3 .4?P 
36! “—() . 
€ 0.365 (>) 0.0464 


Here p is the electric resistivity of the metal in ohm 
centimeters and \ is the wavelength in centimeters. 
Table 2 and figure 4 include the theoretical emis- 
sivities calculated using this equation and_ the 
resistivity value 46.52 wohm em measured at 1,400° 
K by Wm. F. Roeser [4]. 

At long wavelengths the measured emissivity is in 
excellent agreement with the theoretical values. 
There is, however, an unexpected rise in the observed 
values around 8, 9, and 10 uw. This may be due to 
very small amounts of impurities on the surface of 
the metal. The deviation at shorter wavelengths 
may be partially due to the fact that the theoretical 
curve is expected to apply only to the long wave- 
length region. 


4. Conclusion 


The apparatus described in this paper may be used 
to measure the emissivity of metallic substances or 
substances capable of forming tightly adhering 





| 
| 
| 
| 
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coatings on metals to within about +3 percent for 
substances with a high emittance and to about +5 
percent (of measured value) for substances with a 
low emittance. For highly emissive substances the 
temperature Measurement is probably the principal 
source of error while for low-emittance specimens 
the amplifier stability becomes a major source of 
error. The reproducibility of the results has been 
found to be very good. It was checked by making 
repeated measurements on the same sbecimen and on 
different specimens of the same material. In addi- 
tion a different operator was used to check some of 
the work. 

One could also measure the emissivity by changing 
the blackbody temperature so as to equate the radia- 
tion of sample and blackbody. The principal source 
of error would then be due to the temperature mea- 
surement for all emissivity values. <A big disad- 
vantage of this procedure is the slowness of making 
measurements due to the large thermal inertia of the 
blackbody. 


In conclusion the authors express their gratitude to 
Wm. F. Roeser of the Building Technology Division 
for helpful suggestions; also to Fred Shorten for 
assistance in taking some of the observations. 
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Electrostatic Deflection Plates for Cathode-Ray Tubes. 
I Design of Single-Bend Deflection Plates With Parallel 
Entrance Sections. 

II. Deflection Defocusing Distortion of Single-Bend Deflec- 
tion Plates With Parallel Entrance Sections. 
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In section I, a plate design system is offered which allows rapid and accurate determina- 


tion of mechanical plate parameters to achieve given electrical plate characteristics. 
design is suitable for single-bend plates with parallel entrance sections. 


The 
The design curves 


were calculated under the conventional assumptions of small deflection theory, but correc- 


tions for the entrance and exit fields are included. 


and corrections is indicated in an appendix. 


The method of calculating the curves 


Section II deals with deflection defocusing of such plates, a convenient formula for 
calculating the defocusing distortions of single-bend deflection plates is derived and compared 


with experiment. 


tion is shown to lead to increased distortions. 


I. Design of Single-Bend Deflection Plates 
With Parallel Entrance Sections 


1. Introduction 


Designing a cathode-ray tube is as much an art as 
a science; experience and the correct choice of com- 
promises between various parameters are both im- 
portant. Experience suggests that the best tubes 
are designed in close cooperation between instru- 
ment designers and tube engineers. In many in- 
stances the tube designer in these discussions has 
nothing to go by except his experience, and any 
information of a quantitative nature is welcome. 
Such information is provided for the gun_ per- 
formance, for instance, by the work of H. Moss [1],? 
and it is hoped that the design system offered here 
will contribute to establishing this kind of informa- 
tion for the design of plates. 

The design of cathode-ray tube deflection plates 
of specific characteristics is complicated by the choice 
of parameters available. Thus the shape of the plate 
may range from parallel plates to continuously flared 
plates, and the initial spacings from 0.030 in. to 
0.1 in. 

Given a definite plate shape, one can solve the 
equation of motion and thus determine the sensi- 
tivity and sean. The converse is not true, however, 
and tube designers have long needed a method which 
enables them to decide on suitable compromises 
without having to check design possibilities by labo- 


' Formerly, Tektronics—Portland, Oregon. 
DuMont Laboratories in 1953 and 1954. 
? Figures in brackets indicate the literature references at the end of section I. 





This type of distortion is proportional to the square of the deflection 
angle and is shown to be sensitive to plate design. 
lower distortions while the most ‘economical’ plates yield larger distortions, 


In general, long deflection plates give 
Post deflec- 


rious experiment. The design method offered here 
has been checked in many experiments over a span of 
five years and has been found accurate within the 
limits imposed by mechanical tube assembly. 


2. Qualitative Considerations 


A successful cathode-ray tube is one which will 
allow the ultimate utilization of the characteristics 
of the vertical amplifier for which it was designed. 
Such amplifiers are generally limited in voltage swing 
and bandwidth by the capacitance of the load which 
they have to drive. The requirements on the 
cathode-ray tube are therefore to provide best de- 
flection sensitivity consistent with the required sean, 
capacitance, and writing speed. The most sensitive 
plates would be flared in such a manner that the 
fully deflected beam smoothly follows the positive 
deflection plate. Such optimum plate shapes have 
been described in the literature [2]. 

In practice continuous curvature is difficult to 
obtain, and approximation by bending is preferred. 
For the present design, single-bend plates with paral- 
lel entrance sections have been chosen. Such plates 
approach the optimum very closely, particularly in 


scan-limited tubes, which are used in the most 
advanced applications. 
3. Theory 


The most important design characteristics of a 
plate pair are its inverse sensitivity (say, in volts per 
inch) and the sean (in inches). The product of these 


This work was done at the A. B. | parameters determines the maximum potential dif- 


ference applied to the plates. 
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When, 
ranged symmetrically 
then one-half of this potential, 
dimensions of the plate. 


the 
about 


as is usually case, the plates are ar- 

the undeflected beam, 
Wie. determines the 
Consider the parallel sec- 
tion first, as in figure 1. We arbitrarily stipulate 
that the width of the beam, w, should be 2/3. s, 
where s is the spacing of the plates, and that \ of the 








beam should be cut off by the plate edge. The 
deflection at the exit end is then given by: 
> la { ly y . /i /} 
Y,=8/H5=—e/mM (J max/S)@E ¢ (1) 
4 
Zi ae) 
je = A 
— —_- paar cg A 
Ww ad ir aller cae ae csi oe 

















Frat RE g The parallel section. 


whiere a is the length of plate. 
to establish a as a function of s 
During transit 
further deflected. 
section alone gives, 


This equation serves 
and Vmax. 

in the flared section the beam is 
The deflection due to the flared 
with reference to figure 2 


Yo m| | (Vinax/ar) dt dts (2) 


where we have assumed cylindrical field symmetry. 
The total deflection Ys at the exit end of the combined 
plate is due to the deflection in each section plus the 
displacement in the flared section due to acceleration 
in the first, 1. 


Y3=¢ m| : (Vv, max/S) dt dt 
+ ¢ mt | (Vaas/s) dt | 
vr, "lp z 
r ¢ m| | Ey wees ar) dt dt; (33) 
where ¢t,, is the transit time from ¢ to b. We now 


stipulate arbitrarily that at the exit end of the 
combined plate one-half of the beam should be cut 
off, i 


Y3 a)/2, (4) 
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Figure 2. The flared section. 


Thus, eqs (3) and (4) serve to give a as a function 
of s, a, Vinax, and any arbitrary length of plate 
(L=a+b—c). In practice, a must be determined 
graphically. The sensitivity or inverse sensitivity 
may now be calculated for known screen distances, 
Additional assumptions, corrections, and the general 
formula are given in appendix 1 and 2. It should be 
noted that the determination of plate shi ape involves 
the simultaneous solution of eqs 1 to 4 and eq (al) 
of appendix 38, after arbitrary decisions on spacing 
or le ngth have been made. 

On the other hand, the designs presented here are 
based on a large number of theoretical calculations 
in which the sensitivity was calculated as above after 
spacings, total lengths, and maximum voltages were 
arbitrarily chosen. As a result of calculating a 
large number of designs and tabulating the results, 
sets of curves are obtained which permit the selection 
of a plate shape for a given sensitivity and sean. 


4. Method of Design 


There are 3 sets of curves concerned with plate 
design parameters. All are based on the character- 
istic peak to peak plate cutoff voltage of 2 Vinay. 
This may be defined as follows: Given a standard 
accelerating potential of 1 kv, 2V. max is the product of 
the desired total scan, times the inverse sensitivity 
(in volts per unit deflection). 

The first set of curves Al to A7 gives the inverse 
sensitivity as a function of There are 
several sheets of curves each for a given total length 
of plate LZ. On each sheet there are several curves 
each for a specific entrance spacing s. The first set 
of curves thus serves to determine immediately what 
length and spacing are necessary to achieve a given 
scan and sensititivity. 

The remaining two sets of curves B and C1 to C7 
serve to establish the design completely by giving 
the length of the straight section and the angle of 
bend as functions of the par ameter 2 Vinx and the 
length and spacing chosen from the first set. Thus 
it is not necessary to consult these sets for the initial 
discussion of a new cathode-ray tube design. 


2 \ max. 





5. Post Acceleration 


Post acceleration does not affect the value of V max 
for a given plate but increases the inverse sensitivity 
by an amount variously called intensifier factor or 
compression factor. The factors vary depending on 
the bulb shape and the shape and the degree of 
continuity of the post accelerator. It differs also 
for the front and back plates of the same tube. To 
facilitate the design of post deflection accelerator 
tubes a set of curves D giv Ing compression ratios for 
various tube types as a function of accelerator po- 
tential to anode potential ratio is included. 

The compression factors for types quoted will 
apply approximately to all tubes of similar design. 
Type | has a straight-sided (T-shaped) bulb with 3 
accelerator bands; type IT has 1 accelerator band on 
a conical bulb; and type IIL uses a continuous accel- 
erator in a bulb of exponential shape. For exact 
information the factors must be established experi- 
mentally for any given bulb shape, accelerator bands, 
and tube geometry. In general the continuous 
accelerator spiral gives the lowest loss in sensitivity. 


6. Correction for Tube Length 


The values of inverse sensitivity must be corrected 
by the ratio of the effective distance to the screen to 
the effective distance applicable for each sheet. 
This may be done as follows: The effective distance, 
with sufficient accuracy, may be taken as the distance 
of the screen from the exit end of the plate plus half 
the plate length. The distance D quoted on each 
sheet is the distance of the exit end of the plates on 
that sheet from the screen of the tube. The cor- 
responding assumed effective distances therefore are 
greater in each case by about one-half of the plate 
length applicable to that sheet. The reason for this 
arrangement is that the design is needed mostly for 
the vertical plates of sensitive instrument tubes, and 
in considering the effect of changing plate lengths it 
is convenient to fit all designs to a standard length of 
tube. 


7. Beam Width 


The pl: ite designs are based on an assumed beam 
width of % of the spacing s of the parallel section of 
the plates. If the limiting aperature size is taken 
in this manner cutoff is sharp and clean. For the 
highest performance types of tubes one might sacrifice 
edge performance, and values of 4 s for the diameter 
of the limiting aperture have been lused in some 
tubes with consequent gain in brightness and defini- 
tion in the center area of the sereen. 


Appendix I. Summary of Steps in Designing a Plate 


lL. Hstablish the desired characteristics. 
potential; (b) Post deflection accelerator 
Desired inverse sensitivity; (d) 
(e) Tube length. 


Anode 
(Cc) 
and 


(a) 
ratio; 
Desired scan; 
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2. Find 2 Viner. Multiply the inverse sensitivity 
— by the desired sean. 


Normalize inverse sensitivity and 2 Ving. This 
is in by dividing both quantities by the anode 
potential in kilovolts. 


Adjust normalized inverse sensitivity for post 


deflection. This is done by finding the compression 
factor from curves D and reducing the inverse 


sensitivity by it. 

5. Adjust normalized inverse sensitivity for tube 
length. It is now necessary to select an approximate 
value of plate length and to find the corresponding 
value of D from the curves. The value for inverse 
sensitivity is then adjusted as outlined previously. 

Steps 1 to 5 are easily accomplished, and the infor- 
mation is now available for a complete design. 

6. Hstablish a spacing s and complete design. By 
consulting the set of curves A, a spacing may be 
selected for each of various possible plate lengths. 
This involves a compromise in terms of beam size, 
etc. When this choice has been made the straight 
length and bend angle may be found directly from 
the remaining two sets of curves. If it is found that 
the design chosen leads to unrealistic spacings or 
length it is necessary to modify the electrical char- 
acteristics and to try again. Since the whole 
procedure takes only a short time, a suitable com- 
promise can be found quickly. 


Appendix II. Corrections 


An accurate small signal theory is given by Hutter 
[3]. More specific formulas [2, 4] have appeared in 
the literature from time to time, and the present 
method of design based on the conventional 
assumptions of constant transit velocity and cylin- 
drical field geometry in the flared part. Three 
corrections are, however, included in the final cal- 
culation which should be stated here. 

(1) It is assumed that the field of the parallel 
section extends a distance s/2 beyond the plates. 

(2) The electron with maximum deflection is 
assumed to travel along an equipotential line close 
to the positive deflection plate. Accordingly the 
field is inclined to the vertical by about the angle 
which the plate subtends with the axis. This means 
that the transverse acceleration results in a reduction 
of axial velocity. The term describing the trans- 
verse acceleration is corrected by the cosine of the 
plate angle, and the axial velocity through the drift 
space is corrected accordingly. 

(3) It is also assumed that the axial motion is 
further retarded at the exit end while the transverse 
motion remains unchanged. Thus an electron just 
passing the edge of the positive plate has energy 
e(V + Vinax/2) but a short distance farther its energy 
is again eV and remains so to the screen. The cor- 
responding field is axially directed since it usually 
terminates on grounded structures just ahead of the 
plates. 


The axial velocity is reduced according’ y. 
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Appendix‘III. Theory 


The general formula used for the design was: 


v(t) rnd) 


(al) 


where 

a@ is the length of the straight section plus 14 of the 
spacing, s, 

6 is the terminal length of the flared section in 
polar coordinates, 

¢c is the initial point of the flared section in polar 
coordinates, 

D is the distance to the screen of the exit end of 
the plates, 

+-V is the anode potential, 

V+U/2, V—U/2 are the plate potentials, 

y is the deflection at the screen, 

F’. is the compression ratio, 

s is the spacing of the plates, and 

a is the plate angle. 











Ficure 3. Flared plate with parallel entrance. 


Formula 1 was derived by integrating and summing 
the equations of motion of the electron in the cor- 
responding regions of the plate and the drift space. 

In the parallel section, for example, we have for 
transverse velocity v, and position y, 


ty 
vy= (el ms) dt, (a2) 
0 
; (a (te / 
y=(el ms) | | dt dt. (a3) 
0 0 
For the flared region these quantities are 
, “t; 
v,=(el me) | (1/r)dt, (a4) 
Jt, 
: ‘. >... 
y= (el ma) | | (1/r)dt dt. (a5) 
Jt,dt 


r c 





We also have, according to our assumption, the 
energy equation determining the axial velocity. 
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mv,?=22eV, (a6) 
and therefore, 

| , = er 
| dt=dzr/y2eV/m, (a7) 
| applying eq (a2) to (a5) to the corresponding regions, 
| integrating by means of (a7) and applying the cor- 
| rections outlined before there results eq (a1). 
| 8. References 
| [1] H. Moss, The electron gun of the cathode ray tube, 
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II. Deflection Defocusing Distortion of Single- 
Bend Deflection Plates With Parallel En- 


| trance Sections 


1. Introduction 


Cathode-ray oscilloscope displays generally suffer 
from a broadening of the trace with increasing 
deflection [1].!. This is known as deflection distortion 
and has been treated in the literature in a general 
way [1], 2, 3]. 

Many standard cathode-ray tubes use single-bend 
plates with parallel entrance sections since these 
plates are both convenient and sensitive. A graphical 

| design method for such plates has been described in 
section 1. We now deal primarily with the calcula- 
tion of the deflection distortion due to plates of this 
type. 


| 
| 


| 2. Distortions Due to the Parallel Entrance Section 
Let V be the mean potential in the plate 
region with respect to the cathode, 
[’ be the deflection potential across the 
plates, 
a be the length of the straight section, 
s be the plate spacing, 
e be the charge of the electron, 


m be the mass of the electron, 
F,, F,, F, be field components, 
0}, Vr, 0, be velocity components, 
FE be a symbol for energy, and 


a 


vb be the width of the beam. 

Since electron energy at any point is given by the 
potential of the point in the field, we see that electrons 
entering the plate at A and B (fig. 1) respectively 


| 1 Figures in brackets indicate the literature references at the end of Section II. 
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have energies differing by 


AE=eUw/s. (1) 

| 
Further, we see from figure 1 that this difference is | 
nearly all due to the difference in horizontal velocity, 
since the vertical field components acting on the 
electrons entering the field are equal while the 
horizontal components are opposite. 

We have during transit: 


field= Fy= U/s (2) 
force=el’/s (3) 
accel=el’/ms (4) 
time=a/v, (5) 
therefore, ' 
vy=eal’/mSV,, (6) 


whence the exit angle for small deflections is given by: 


6~tan 6=eUa/msv, 


(7) 
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FIGURE 1. 


The fringe field in the parallel section. 


But v, differs for the two electrons, and to find the 
distortion we differentiate 6 with respect to the en- 
trance ordinate y. | 
| 
| 





We note that generally 


E=mv?/2 (8) 
and 
E=e(V+yU/s) (9) 
or, combining (8) and (9) 
y= 2e(V+yU/s)/m (10) 





and since, as we have seen, for purposes of energy | 
computations v equals v,, we get for 4 


| 
| 


6=al/2s(V+yU/s) (11) | 
and by differentiation | 
dé al? (12) | 


dy _.28*(V-+U/s)? 


Using w for the separation dy of the two extreme 
electrons and substituting eq (11) in (12) we find: 


(13) 


Thus the distortion varies as the square of the deflec- 
tion angle and linearly with beam width. We also 
note that it varies inversely as the plate length. 


Ad= — 2w6?/a. 


3. Flared Section 


We will assume that the field in the flared section is 
cylindrically symmetrical. 

One further assumption will be made which re- 
quires some discussion. We will again assume that 
the entire energy difference for the extreme electrons 
is due to the difference in v,. This assumption may 
be justified as follows. The mean energy difference 
of the electrons at extreme ends of the cross section 
is of the order of U divided by the ratio of plate 
spacing at some arbitrary point to beam width. 
Since the deflecting field decreases as the spacing 
increases we may assume the effective ratio to be 
about 2 or 3 for plates in which the initial spacing 
is comparable to the beam width. For such plates 
the Y-directed velocity after transit through the 
plates corresponds to a kinetic energy of the order 
of U. Thus it follows that the energy differences 
across the beam are substantial fractions of the 
deflection energy. But it may be observed experi- 
mentally that the distortion is only a small fraction 
of the total deflection. Thus the energy differential 
involved in the vertical motion of the electrons in the 
cross section is negligible. 

Let us express this analytically: 

In figure 2, let 

c=the radius at the entrance, 

b=the radius at the exit, 

a=the angle of flare, and 


g=the final opening. 
Then 
mv? /2=ea(V+yU/ar—e). (14) 
ccmataanes | 
. 




















FIGURE 2. 


The flared section. 





Here is a correction for the fact that some energy is 


contained in the vertical motion of the electron. 
Our assumption now states that e may be larger than 
ylU/axr, but the differential of « with respect to y is 
negligible compared to the differential of (yl’/ax), 
at least in those parts of the field where deflection 
is effective. With these assumptions we proceed as 
before: 


F,=(U)/ar) approx (15) 


(acceleration) y=el’/mar. (16) 


But since the acceleration now depends on x we have 
to integrate to find v,: 


ty 
vy=el’/ma | (1/xr) dt (17) 
Jt. 
and since dt=dzr/v, (18) 
Dh 
vy Uma | 1/(v,r dy. (19) 
c 
For small angles we may sum the increment: 
*D 
6 (el/ma) | i/(vPr)dr. (20) 
Substituting for v2 from (14) we have 
a dx 
: anes - (21) 
2aJ. 2(V+yl/ar—e) 


To find”the defocusing effect we differentiate (21) 
with respect to y, remembering that de/dy=0; hence, 
o6 (OLY " Udr 
Oy 2a). (V+yU/ar—e)*ax? 


(22) 


We may now integrate (21) and (22), disregarding 
small quantities, to get 


6=(U'/2aV) log (b/c), (23) 
and 
IS plinth 
6——d -( ) (24) 
. gy 2a°V*\_ be 


Thus, dé depends on UL, and @ on U’; hence in the in- 
clined section dé again depends on @ and on w=dy, 
buttits‘*dependence on plate shape is now complicated 
as it involves log (b/c). 


4. Total Distortion 


A’ formula may now be established which allows 
direct investigation of the distortion of various plates 
designs. Let D be the distance to the screen; then 
the total broadening of the line will be 

Ay=D(d6,+d6,) 


(25) 


where d6,, /@ are the distortions of the straight and 





flared sections, respectively. 


Substituting in (25) 
from (12) and (24) we 


get, after some rearranging: 


wDL"a (q—s) 
Ay=——3>7 | l+83 - , 26 
id 237)? [ i aga (26) 
Here c and 6 have been eliminated by 
ac=8, (27) 
and 
ab=yq, (27a) 


respectively. For practical designs (26) may be 
further simplified since g is much larger than s, Le., 


(i435) 


Equation (28) is useful in comparing the distor- 
tions of various possible designs. Thus, when using 
the plate designs described in section I of this paper, 
the parameters in eq (28) may be found from the 
electrical characteristics required of the plate, and it 
is possible to take distortions into consideration when 
discussing a particular design. An example of the 
application of formula 28 to three experimental plates 
for a tube is given in table 1. Plate C is about twice 


wDU"a 
2V%s 


? — 


(28) 


TABLE 1 
A B Cc 
W= 030 in. W= 030 in. W=030 in. 
D= 10.50 in. D= 12.25in. D= 11.25in. 
U= 150 v (full U= 125 v (full U= 85 v (full 
scan scan) scan) 
a 220 in a . 210 in. a= . 300 in. 
8 060 in. 8 O50 in. &= . 060 in 
V=1675 v V=1675 v V=1675 v 
a 600 rad a= . 535 rad a= . 500 rad 
& & 
450) . 445 : . 400 
ad ad aa - 
Ay 128 Ay 124 in. Ay= 036 in 


as long as the other plates and may be seen to have 
only about one-fourth as much distortion. For a 
given plate length, increased sensitivity gives 
slightly better distortion characteristics, but gen- 
erally distortions are determined by the desired 
characteristics of the plate. 


5. Post Accelerator Field 


Post acceleration generally reduces a scan pattern 
by the so-called compression factor shown in figure D 
of the plate design curves (see 5, section I). Since 
distortions are essentially a scan pattern they, too, 
are reduced by post acceleration, and eq (28) must 
be modified accordingly: 


DU Ss 
“a ai ( ; ro. ): 


In evaluating the distortions of monoaccelerators 
versus post accelerators it must be remembered that 
the deflection voltage requirements in the latter in- 


Ay (29) 
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crease by the factor Ff. In eq (29) therefore, at 
equal deflection, Ay would be increased by a factor 
of F. since for a given plate and for a fixed deflection, 
U7 would have to be increased as post acceleration is 
applied. 


6. Agreement With Experiment 


Maloff and Eppstein [1] have stated that distor- 
tions due to deflection in well-designed tubes are 
several times less severe than those to be expected 
from elementary theory. It is interesting therefore 
to see how eqs (28) or (29) compare with experiment. 
A typical experiment is shown in figure 3. The 

















Figure 3. Deflection defocusing. 


plate here is similar to plate A of table 1 but a com- 
pression factor of /,=2 applies. The agreement 
between calculated and experimental values are seen 
to be good in this instance. When making observa- 
tions of deflection defocusing some points have to be 
carefully considered. Most modern plates are mar- 
ginal, i.e., they intercept some of the beam toward 
the edge of the scan, and therefore the spot there 
tends to be less distorted than it would be theoreti- 
cally. Another point of experimental significance 
is the mode of focusing. Consider for instance 
figure 4 in which 2 positions are shown for the screen. 














EF ~ FEN B 
DJUSTMENT FOR FOCUS ADJUSTMENT FOR 
R TION BEST CENTER FOCUS 


Figure 4. Two possible focus adjustments relative to screen 


position. 


Position a@ shows an “underfocused’”’ beam while 
position 6 corresponds to an overfocused beam. It 
is clear that deflection distortion, which consists of 
a focusing action with deflection, will be more 
effective in position 6 than in position a. The con- 
clusion is that defocusing depends on the adjustment 
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of the focusing lens and on the nature of beam forma- 
tion about which not much is known. 


7. Corrections for Defocusing 


Qualitative reasoning suggests that no effective 
correction can be found which does not simultane- 
ously result in severe deflection linearity distortion. 
This may be seen as follows: Since A y 1s a negative 
function of y, deflection results in overfocusing. 
Nothing can be done about this in a nearly parallel 
section which is present in any deflection plate, for 
in such a region the assumptions made here are 
rigidly true. 

At the exit end of the plates however, where the 
beam is already deflected across the axis, our assump- 
tions do not hold, and the field may be shaped in such 
a manner as to deflect the electron traveling closer 
to the positive plate more than the electron at the 
other end of the cross section. If this is done it is 
clear that the correction is applied in a region of 
rapidly increasing separation of the entire beam 
from the axis. Thus what is required here is a 
variation of the effective strength of the deflection 
field in a transverse direction which is effective over 
the cross section of the beam. Consider now the 
effect of such a field on the central electron in a beam. 
As the deflection is increased the central electron 
travels in regions of increasingly effective deflection 
fields, and consequently the overall deflection 
becomes nonlinear. 

The foregoing conclusion is particularly true of 
arrangements employing corrective static lenses such 
as positive wires lined up parallel to the exit edge of 
the plates. 

In contrast with the qualitative argument given 
above, detailed mathematical treatment has led 
other authors [3, 4,5], to propose plate shapes for 
which various degrees of correction are claimed. 


8. Conclusions 


Deflection defocusing in single-bend plates with 
parallel entrance sections is approximately propor- 
tional to the square of the deflection, it is linearly 
dependent on beam width, and is roughly inversely 
proportional to plate length and is therefore suscep- 
tible to changes in plate design. 
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The Functional Synthesis of Linear Plots' 
J. P. Vinti and R. F. Dressler ” 


(January 26, 1960) 


In practical engineering or experimental work one often encounters a function F of 
many variables, F (z's, y’s, 2’s), represented only by the families of curves obtained by 
plotting F against each of the z’s on Cartesian graph paper, against each of the y’s on semi- 
log paper, and against each of the z’s on double-log paper. It frequently happens that 
these curves are all approximately straight lines over a limited range of interest. On the 
assumption that they are all true straight lines, the present note shows how to synthesize 
all the graphical representations, for any number of parameters, into the most general 
formula possible, expressing F as the product of a multilinear function of the z’s and the 
exponential of a constant-free multilinear function of the y’s and of the log z’s, the coeffi- 
cients in both multilinear functions being independent of the x’s, y’s, and 2’s. 


1. Introduction 


In an investigation by one of the authors on assemblages of elastic shells, the various 
results for certain components of stress and displacement exhibited approximately linear 
behavior over the limited ranges of interest of the relevant parameters, when plotted on Car- 
tesian graph paper, semi-log paper, or double-log paper. These results, derived from various 
numerical calculations and corroboratory experiments, are functions of many parameters, such as 
the various geometrical ratios defining shell shapes, the shell thickness, Poisson’s ratio, the num- 
ber of coupled shells, and other significant quantities. Rather than to retain these extensive 
results in the form of cumbersome families of graphs, it was desirable and useful to combine 
them in such a way as to obtain a single explicit formula for each dependent quantity, in terms 
of the above-mentioned independent parameters. 

The problem of combining such results is, of course, not peculiar to investigations in 
elasticity but often arises in experimental or engineering work of any nature. One can very 
easily construct simple functions which will exhibit some of these “linear’’ properties, but the 
most general answer to the inverse problem is less obvious, especially when the number of 
parameters is large. Then an unsystematic attempt to effect such a synthesis may prove 
infeasible or incomplete, and furnish no assurance that one has indeed constructed the most 
general function with these properties. Our present note, therefore, formulates the general 
problem for any number of variables in a precise manner and derives its most general solution. 
For the simple case of only three parameters, the solution is illustrated by a few examples. 


2. The General Problem 


In the following developments we shall assume that all the pertinent quantities have been 
grouped into independent dimensionless combinations, so as to take advantage of whatever infor- 
mation is provided by the Buckingham Pi theorem.’ Let us then consider a dimensionless func- 
tion F'of the dimensionless independent variables x;, 22, . . ., Lm3 Yi, Yo, - - +» Yps 215 225 


Mes 
and U;, UW, .. ., U,, such that a straight line results when we plot F against any z on ordinary 
Cartesian graph paper, against any y on semi-log paper, or against any z on double-log paper. 
Let the w’s denote all other dimensionless variables on which F may depend, but for which 
no such “‘linear’’ property exists. It is understood, of course, that in practice such linearity 
may hold only approximately, and only for a certain bounded range of the variable used as 
abscissa, and only for certain bounded ranges of the remaining variables, which appear as 
parameters. With no loss of generality, however, we relax these restrictions, assuming that 
each plot is a straight line for all values of abscissa and parameters. 


1 This research was supported by the United States Air Force, through the Office of Scientific Research of the Research and Development 
Command. 


2 On leave to Institute of Mathematical Sciences, New York University, New York, N.Y. 
3 Pp. W. Bridgman, Dimensional analysis, p. 36 (Yale Univ. Press, New Haven, Conn., 1931). 
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We put 


Seas, Mike=Fyis, -. + M2,—Fe, (n=p+4q); (1) 
then F(a,, . . . 2m, Yi, - - » Yn, U) is a function such that the plot of F against any z or of InF 
against any y isa straight line on Cartesian graph paper. Here uw denotes the set w;, Wo, . . ., Uy. 
Then 

InF=A,yy,.+ B, feel. 2, ...-.) 2) (2) 
where A, and B, are functions of 4; . . . Sm, Yi... Yeoiy Year + + + Yn, and of the w’s. We 


then ask: if InF is linear in each y when the other y’s are all held constant, what is the most 
general form for InF, as a function of the y’s, that will represent such a property? 

To answer this preliminary question, we first recall the definition of a multilinear function 
G, of several variables ¢,, tf, . . ., ¢,, as a funetion which is the sum of a constant and a linear 
combination of all the products of the t’s taken one at a time, two at a time, . . ., s ata time, 
without repetition. For example, if s=3, 


G (ty ,t2,t3) =o t+ Ayty + At, + asts + bytots 4 bot st + bytit, + etitots. (3) 

Any such multilinear function satisfies the differential equations 
0°G/ot?7=—0 (A=1,2.. . »@). (4) 
(When the constant term vanishes, we term G a constant-free multilinear function.) It is 
then easily shown by induction that G@(t,,f, . . ., ¢t,) is the most general function of ¢,,f., .. ., f, 

which is linear in each f. 
2.1. The Synthesis for the y’s Alone 
On applying these considerations to (2), we find that 

InF=9(2,,%2, . . -» Lm, U)+N@=™(y1,4o, . . -> Yu), (5) 


where N“*™ denotes a general constant-free multilinear function of the y’s, with coefficients 
which are, a priori, functions of 7,, . . ., 2m, u. On placing 


OXP 9=J (Zits, . . -» Sa;%), (6) 


we obtain 


u 


F’'=f(2,.20, . . «. Inst) exp N“ 
e - b] 


, (7 
which gives the synthesis of the linearities of In¥' versus the y’s. 
2.2. The Complete Synthesis 
Since F is linear in each z, it follows from the property of multilinear functions that 
Fur P P(e, a, . . « sBa), (8) 


where P™™” is a general multilinear function of the x’s, the coefficients being functions of the 
ysandthewu’s. Then 


Par .« + « Se SLO NOV RPO" (2, 2, . . « » Babs (9) 

where by (4) 
oN /dyi=0 (> ie} (10) 
oP /dxr?=0 fo) a |} (11) 





We next show that f is a multilinear function of the z’s. To do so, put every y equal to 
zero in (9). Since N contains no constant term, it then vanishes, so that (9) becomes 


ST (x1, 22, Fe t > “my u) = M™ (2, 20, ee. » Lm). (12) 


Here M“™ is simply the expression for P™” with each coefficient evaluated at each y=0, so 
that it is a multilinear function of the z’s with coefficients depending only upon the w’s. 


We now show that the coefficients in N®™ are independent of the z’s, as follows. Insert 
(12) into (7), so that 
F= Mm exp N& u) (13) 
and require that (13) satisfy (11). It follows that 
OM oN, OMON oN\’ 
—,+M —~4+2 M =f 14 
Or; ” Ors - Or; ee es ) (14) 


where we have omitted the superscripts for convenience. 


Since \M/ is multilinear in the z’s, 
it follows from (4) that 


0°. /ozr?=0, (15) 

so that 
oN ,, MON oN\?_, 
- Or? a Or; Or; +M (5 ) Ae: (16) 


If we now differentiate (16) twice with respect to y,, we find with use of (10) that 


oN \2 
9 he ~ 
2M (37 >,,) 0, (17) 
whence 
0 (ON f gen 1.2...) mM, 
02° 
OY: ( Or; ) .. | oe * (18) 


Equation (18) means that 0N/dz; can be a function only of the z’s and the w’s. But N and 
ON/oxr; are constant-free multilinear functions of the y’s, with coefficients that are, a priori, 
functions of the z’s and the w’s. These results are compatible if and only if 


ON /dxr;=0 ct re) (19) 
so that each coefficient in the constant-free multilinear function N must be independent of 
the z’s. 

We may thus rewrite (13) as 
P=M™(z,,%, . . . Su) CE NM (yess: . , Ba). (20) 


When we return to the original formulation of the problem in terms of the 2’s, y’s, and 2’s, 
it follows that the most general functional form for F is given by 


F=M“ (a2, . . . »%a) exp N™ (gi Ms, « ~.. 5 Hoe En Se, «.. . , Sad (21) 
where M™ is a general multilinear function of the z’s and N™ a general constant-free multi- 


linear function of the y’s and Inz’s, the coefficients in both being functions only of the w’s. 


3. Some Elementary Examples 


As short illustrations of the general result (21), we append a few cases where F has linear 
plots against only three variables. For each we list the specific form that (21) assumes and 
the slopes and intercepts on the appropriate plots. By comparing the behavior of the slopes 
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and intercepts in the various experimentally given families of curves that define the function F 
with these listed formulas, one can readily determine which coefficients vanish, if any, and thus 
obtain a specific formula for F in any actual case. Here “In’’ denotes a natural logarithm 
and “log’’ a common logarithm. For an z or a y the intercept is taken at zero, while for a z 
it is taken at z=1. For the logarithmic plots the slopes and intercepts are those of log F. 
In the following formulas, it is understood that the constants may be functions of the w’s. 


(a) X, Yi» Yo 
F=(kyr+hkz) exp(ayy; +asy2+ by,y2) 


Cartesian plot versus 2: 
slope S=k, exp(ayy: + a2y2t byrye) 
intercept I=k, exp(ayy,; + d2y2+ by yo) 


Semi-log plot versus y;: 


slope S= 0.4343 (a, +- by») 
intercept I=log(k,x+-kz) +-0.4343asy. 
(b) x,y, 2 


F= (k,a-4 k, \é yu atgtby 


Cartesian plot versus x: 
slope S=k,etzgt:ty 
j ‘ ; by 
intercept T=ke™ 22 


Semi-log plot versus y: 


slope S=0.4343a,+6 log z 
intercept T=log(kye+k;) +a log z 


Double-log plot versus 2: 


slope S=a,+ by 
intercept I=log(kya+ kz) +-0.4343a,y 
(c) Yi» Y2s Z 


Wun rem +4 Yo+b Vy V2 13+ bits + bo Yat CY Yo 


Semi-log plot versus y;: 


slope S=0.4343 (a, + byy2) + (6; +-cy,)log z 
intercept I=log k+-0.4343a2y2.+ (a3+boy2)log z 


Double-log plot versus 2: 


slope S=az+by,4+ boys +eyys 
intercept I=log k+-0.4343 (ayy, + doyo+ bsyry2) 
(d) Ys Zi» Z2 


Fake 2,20 wz,t3t ov] 
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where 


H= 2, ©steninzs — 2, Ostepdinz, 
Semi-log plot versus y: 
slope S=0.4343a,+ 5; log 2,+b. log 2.+2.303¢ log 2, log 2 
intercept I=log k+-az log 2+ a3 log 2.+2.3036; log z, log 22 
Double-log plot versus 2;: 
slope S=a.+by+2.303(b3;+-cy)log 2. 
intercept I=log k+-0.4343a,y+ (a;+b.y)log 2. 
Wasuineton, D.C. (Paper 64C2-31) 
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Radiation Field From a Rectangular Source’ 


J. H. Hubbell, R. L. Bach, and J. C. Lamkin 


(January 12, 1960) 


Many radiation shielding problems involve calculations of the response of an isotropic 
detector to radiation of arbitrary angular distribution from uniform rectangular sources. 
In calculations of this type the family of integrals /s(cos @dS/r?) P;(cosé@) and the integral 
JS s(dS/r?)exp(—yt/cosé) are frequently encountered, where @ is obliquity with respect to an 
axis perpendicular to the plane containing the rectangular radiant surface, S, r is the distance 
from an element of source area, dS, to the detector, uw is the attenuation coefficient, and ¢ is 
the barrier thickness. Solutions of the first type of integral facilitate use of Legendre expan- 
sion representations of radiation directional distributions, and may also have application 
in other radiant surface studies, such as illumination and heat exchange engineering. The 
second integral relates to exponentially attenuated radiation from a plane isotropic rec- 
tangular source separated from the detector by a layer of material of thickness ¢. 


Formulas, 
expansions, and numerical results are presented. 


1. Introduction 


The penetration of radiation into structures has been the subject of theoretical [1,2] ? and 
experimental [3] investigations as part of a general program at the National Bureau of Standards 
and other agencies. The theoretical investigations are usually performed in two stages: 
(a) penetration and (b) geometry. ‘The penetration stage refers to generation of data relating 
to interaction of the radiation with barrier materials [4], and geometry refers here to the applica- 
tion of such data to specific configurations of source, barrier, and detector. This paper treats 
the geometrical stage of problems involving rectangular sources, e.g., ceilings of buildings or 
shelters. We draw to some extent from similar treatments in illumination [5 to 8] and heat- 
exchange [9] engineering. ‘The methods developed in this paper, in turn, have application to 
rectangular sources of light and heat as well as gamma radiation. 

In the geometrical stage of many shielding problems, the following situation, generalized 
in figure 1, arises: one has a radiation detector and a surface S through which radiation emerges. 
In the schematization of reference [1], the flow of radiation through each element of surface, 
dS, is independent of position on the surface, and no scattering of radiation occurs after emer- 
gence from surface S. The intensity and angular distribution, g(cos @), of radiation crossing 
the plane S is assumed known from the penetration stage of the analysis. The quantity 
g(cos 6), usually called the flux [4], is the radiation intensity per solid angle (e.g., photons per 
cm? per second per steradian) traveling in the direction 6 with respect to the z-axis (see fig. 1). 
If the detector has an isotropic response it will register a value 7, proportional to an integral 
of g(cos @) over the total solid angle, 2s, subtended by the surface S from the detector: 


Is= fy d2s(6) g(cos 6). (1) 


where 


d2,(0) =dS + cos 6/7? 


and ris the distance from surface element dS to the detector. 
Introducing the quantity ¥s(@), which is the total range of azimuthal angle, ¢, defined by 


the intersection of a 6-cone with the surface S, as in figure 1, we may write (1) in polar coordi- 
nates as 


ot 


1 Work supported jointly by the Office of Civil and Defense Mobilization, the Navy Bureau of Yards and Docks, and the Defense Atomic 
Support Agency. 


2 Figures in brackets indicate the literature references at the end of this paper. 
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*1 "Vs 4 
I; | d(cos 0) g(cos 6) dg (2a) 
Bees | 


a 0 


"i 
| d(cos 0)g(cos O)Ws(8). (2b) 


“|G 


As in the example in figure 1, ¥s(@) may consist of more than one ‘‘piece”’ of azimuth, Ajgs(6), 
so that in the general case it is a summation over all such pieces: 


~~ 
— 


This quantity is identical with the “angular response function,” ¥(@), of reference [1] except for a 
factor 2z.° 

The surface S may now be specified as a rectangle, as shown in figure 2, with corners at z 
and y coordinates (0, 0), (w, 0), (w, 2), and (0,7). The detector is on the z-axis, hereafter referred 
to as the “corner position,” at z=h. Since the solid angle integration indicated in eq (1) is 
dependent on relative rather than absolute distances we may reduce the number of parameters 
by dividing all dimensions and coordinates by h, i.e., 


a=wlh, b=l/h, a=ax/h, and B=y/h. 


This representation of the relative geometry, rather than representation by horizontal and 
vertical “eccentricities” as in references [1] and [2], is used here in order that the interchange- 
ability of length and width might be used as a partial check of the results. 

Figure 3 shows the use of these scaled coordinates in describing a typical rectangular 
source problem. Radioactive material is uniformly distributed over the roof slab of an under- 
ground shelter, and the problem is to calculate the detector response to radiation transmitted 
by the rectangular ceiling, or rectangular portions of the ceiling, into the enclosure below. 


3 Berger and Lamkin [1] included in ¥(@) a normalizing factor 1/2x, here omitted. Thus the quantity ¥s(@) defined in (3), above, is equivalent 
to their azimuthal range parameter Ae. 
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There are several ways to evaluate the integral (1) for a rectangle, either in polar form 
(2a,b) or using the cartesian parameters a, 8. Berger and Lamkin [1] performed the inte- 
grations numerically, using the polar form, for a few representative angular distributions 
g(cos 6) and a variety of rectangular shapes. Spencer [2] has approximated the function 
¥;(0), which may now be written for a rectangle as ¥(0,a,b), by step functions so as to reduce 
(2b) to a combination of values of the cumulative integral a d(cos 6) g(cos @). 


In this paper, advantage is taken of the fact that g(cos 6) may usually be represented as a 
sum of Legendre polynomials, i.e., 


g(cos 6) -, ( 1+5) gP1(cos 6), (4) 
=0 \ - 
where 
*1 
1. | d(cos 6) g(cos 0) P; (cos @). (5) 
—1 


Thus the integral (2b) may be expressed as a sum of elementary integrals 
x § : 
1(a,b)=) (1+) gp (a,b), 6) 
=() a _ 


where the p,(a,b)’s are coefficients of the expansion 


v(0,a,b) aes ( 1+5)p (a,b) P; (cos @). (7) 
l=0 = 


The p,(a,b)’s, which are purely geometrical and independent of the penetration stage of the 
problem (i.e., source energy, barrier thickness and material, ete.), may be evaluated as 


pa,b)= dQ;(0)P (cos 6) (8) 


in which @ and b completely specify the surface S. Evaluation of the p,(a,6)’s will be outlined 
in section 4.1. 

The sum (6) converges rapidly if either the g, or the p; converge rapidly. Since the g,’s 
are provided more readily from transport theory [4] than g(cos @), it is desirable to use (6) 
wherever applicable. 

In the case of penetrating radiations, as from fallout, convergence of (6) may be inadequate 
as a result of the presence of the unscattered radiation component, which has an angular distri- 


bution of the form 


g°(7, cos 0)= = exp (—71/cos 6)/cos 0. (9) 
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The superscript ° refers to unscattered radiation, 7 is the slab thickness, ut, in mean-free-paths, 
and o is the surface density of radioactive material. Since the g,°(7)’s evaluated from (5) and 
(9) converge too poorly for use in the expansion (6) when 0<71, an alternative expansion 


I°(r,a,b)=7- >> r”q,(a,b) (10) 


n=0 


is useful. The geometry-dependent coefficients q,(a,b), which are obtained by substituting 
(9) into (1), expanding in 7” and integrating term by term, are given in section 5.2. 


2. Notation and Definitions 


Notation relevant to this paper, much of which has already been introduced, is listed for 
convenience as follows: 
detector=device with response proportional to the radiation flux at a given point in 
space, 
S=symbolic reference to a specific surface lying in the z,y plane, 
h=distance from surface to detector in the z-direction, 
w,l=width and length of a rectangular surface, 
z,y=rectangular coordinates parallel to w and J, respectively, with origin at one 
corner of the rectangle, 
corner position=detector position on a line perpendicular to the rectangle at one of its corners, 
i.e., at (0,0,h), 
a=w/h, b=L/h, 
a=xr/h, B=y/h, 
dS=element of area of the surface S, 
6=obliquity with respect to the z axis, 
g=azimuth of the plane through the z axis and the center of dS, measured from 
the z,2 plane, 
¥s(9)=range of g covered by all surface elements dS from which the detector can be 
seen at obliquity 6, 
¥(0,a,b)=ys(6) when S refers to a rectangle and the detector is in a corner position, 
r=distance from dS to the detector, 
dQ;(8)=solid angle subtended from the detector by element dS at obliquity 6 
=dS* cos 6/r*, 
g(cos @)=radiation flux (e.g., number of photons) traveling in direction @ per unit solid 
angle, unit time, and unit area normal to direction @, 
gi=coefficients of a Legendre polynomial expansion of g(cos 8), 
pi=coefficients of a Legendre polynomial expansion of Ws(4), 
J;=response of an isotropic detector to surface S radiating according to g(cosé), 
above, 
I(a,b)=Is when S is a rectangle and the detector is in a corner position, 
I(a,b;a,8)=IIs when the detector lies at (ha,hB,h). 
The remaining quantities refer specifically to applications involving penetrating radiations, 
typified by the “fallout problem” shown in figure 3: 
7=yt—mean-free-path thickness of layers of material between primary source and the 
electron, where yu is the attenuation coefficient in reciprocal units of the linear 
thickness ¢, 
o/4r=surface density of contaminant, in the case of the fallout problem, in units of 
photons emitted per unit time from unit surface area into unit solid angle, 
I(r,a,b) = I*(7,a,b) + T° (ra,b) 
=the detector response to the unscattered and scattered components, respectively, 
of radiation transmitted by a slab barrier having thickness 7 and a rectangular 
exit surface of dimensions a,b. 
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3. Arbitrary Detector Positions 


Although the results in this paper relate directly to the response of a detector in a corner 
position, the response at an arbitrary position (a,8) (shown in fig. 4) is of more general interest. 

A compounding rule found in standard illuminating engineering texts [10] yields the 
detector response, J(a,b;a,8), at such a position by combining corner-position data from four 
related rectangles: 


I (a,b ;a,8) = I(a,8) + I(a—a, 6—B)+I(a, b—B)+I(a—a, 8B), (11) 
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Fiacure 4. Examples using the scaled coordinates of figure 3, of noncorner positions (a) within the rectangular area 
directly under the source and (b) outside this area. 


in which the terms on the right are obtained by replacing the parameters (a,b) of (6) or (10) by 
the quantities indicated. Since both (6) and (10) are odd functions in ab, i.e., 


I(4a,+b)=—T(a,b), (12) 


the rule in (11) applies not only to detector positions 0<a<a, 0<8<6 as in figure 4a, but to 
any position —° <a<+ ©,—o<B<-+ o as shown in the example in figure 4b. 


4. Legendre Expansion of the Response Function 
4.1. Isotropic Detector 
The primary task performed in this work is the evaluation of the Legendre expansion 
coefficients p,(a,b) by performing the integrations indicated in (8). 


Use of scaled cartesian coordinates leads to the transformations 
COs 6=hfr=1)Vo2+B?+1 (13) 
and 
dQs5(0)=dS-cos 6/r?=dadB-(a? + B?+-1)~*”? (14) 
so that we may write (8) explicitly as 


"4 mb ee 
pi(a,b)= | da| dBc? 4-81)". P(e? +B +1). (15) 
0 0 
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Evaluation of this family of double integrals is outlined in appendix A. ‘The first four p,(a,b)’s 
were found to be 





., , ab | 
y U0 , a= te bh a, > a? 6 
pola in Vat 8'+ 1 (16) 
b tan-! tan7! : I 
a ‘ . an 
pila, ) =f b2+-1 \ b? 1] ya" 7 \a° ] ( 7 ) 
b 1 l 
ee Se. . ; 
pr(a,9) 2/a?+6?+1 La ITR | soled 
and 
| i l 
p3(a,b)= sae ; tan” - T - tan~' ) 
8 \b?+-1 yb?4+1 ya'+1 ya°+1, 
5 t 
+2 : ; tan7'- + a tan-* —— 
“| (b?+-1)? y6?+1 (a?+1)? var+l 


ab ] 1 
ta b24 Les teri) |} (19) 


The po(a,b) term is the result given by Spencer [2] for a diffuse (i.e., g(cos 6) constant) rec- 
tangular source and isotropic detector, and p,(a,b) is the result given by illumination [5-8] and 
heat-exchange [9] engineering studies for a diffuse rectangular source and a surface-type detector 
parallel to the source plane. 

The terms become progressively more complicated with increasing /, but the even and odd 
general terms as given by (38) and (42) in appendix A are amenable to programing for an 
automatic computer. Table 1 contains p,(a,b)’s evaluated using an IBM type 704 computer 
for 0</<13 over a grid of a and b values 0.1 <a<b<20. These may now be inserted in (6) to 
give the detector response J(a,b) in a corner position. 

These p,(a,b)’s may be checked by reconstructing ¥(6,a,b) according to (7) and comparing 
with the exact results given in reference [1]. As shown in figure 5, ¥ (6,a,b) is reasonably well 
described for the case a=1, 6=2, neglecting terms in excess of /=7. If, however, higher g, 
terms are available and contain appreciable information, as for a severely peaked angular 
distribution, greater accuracy may be achieved by including additional terms in (6). 


4.2. Cosine Detector 


The g, coefficients may in principle be generated for any azimuthally independent uniform 
plane source by performing the integration indicated in (5) by analytic, numerical, or graphical 
means. Thus the Legendre expansion method might also be fruitful in similar rectangular 
source studies in illumination and heat-exchange engineering in which the g(cos 6) is known. 
In these studies, however, the radiation field is usually described in terms of a surface-element 
detector (“plane detector’) whose response is proportional to the cosine of the angle of in- 
cidence. 

The orientation of such a plane detector may be specified by its direction cosines A, B, 
and ( with respect to the a, 8, and z axes. Legendre expansion coefficients 7,4,8,c(a,6) 
(analogous to the isotropic detector p,(a,b)’s) may each be written as the sum of three com- 
ponents (see ref [8], eq (17)), 


Pi.a,B.c(4,b)=Ap,_.(a,b) + Bpis(a,b) + Cp,,(a,b), (20) 
where Pia (a,b) and pi (a,b) are the Legendre coefficients for a plane detector perpendicular 
to S and to the a- or 8-axis, respectively, and the p,, (a,b)’s are the coefficients for a detector 


parallel to S, i. e., perpendicular to the z-axis. 
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TABLE 1.* Coefficients pi(a, b) of the 
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Legendre expansion (7) 





of the angular response function y(@, a, b) 
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See footnote at end of table, p. 128. 
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TABLE 1.* Coefficients pi(a, 


l b=20.0 


b) of the Legendre expansion (7) of 


¥ (8, a, b)—Continued 











the angular response function 








b=20.0 b=20.0 h=20.0 b=20.0 
a=0.5 a=1.0 a=2.0 a=5.0 a=10.0 a=20.0 
—_—--- —_—— —___—____—— a= —— —_— — —__—_____ -—— ——— — — —-~----— i 
0 | 4.6302 (—1) 1.1047 (0) | 1 3673 (0) 1.4593 (0) 1. 5002 (0) 
1 3. 5122 | 7.0240 (—1) 7.6994 (—1) 7.8112 (—1) 7.8379 (—1) 
2 2. 0031 2. & | 2.0124 9.9215 (—2) 5.5376 (—2) | 3. 5245 (—2) 
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*Note: The figures in parenthesis in tables 1, 2, 3, and 4 indicate the power of ten by which the numbers alongside and below in the same 
column are to be multiplied; e.g., 9.9012(—3) =0.0099012. 
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The dependence on the angle of incidence of radiation on a plane detector for geometries 
described by pia (a,b) and pig (a,b) reduces to a cos 6 and 8 cos 6, respectively, modifying the 
integrand in (15) to give 


0 


Pra(a, 6)= as adaP(1/y a?+ B4 1)/(a? 


J 0 


and 


pis(a, 6) =| daa BdBP,(1/y a? B? 1) /(e?-4 6?+1)?. (22) 


Evaluation of these integrals should be similar to that for (15). In the geometry described 
by the p,, (a,b)’s, the angles of emission and incidence are identical, so that the integrand in 
(15) requires only an additional factor of cos 6. Thus, using the relation 


cos 6 P, (cos @)={1P,_,(cos 6) +(1+1)Pi4,(cos 6) }/(21+-1), 


the p.. (a,b)’s may be obtained from the p, (a,5)’s in table 1 as 


ary 


P,,(a, b)= {lp,_,(a, 6) +(1+1) 1 4.1(@, 6) }/(21+-1). (24) 


5. Unscattered Radiation Component 
In cases involving penetrating radiations, such as gamma rays from the plane source 


shown in figure 3, it is usually desirable to make a separate calculation of the unscattered 


radiation component, /°(7,a,b)._ Two reasons for this separation of components are: 
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(a) In contrast to the scattered radiation component, J*(7,a,6), the unseattered com- 
ponent, J°(7,a,6), is analytic in nature and independent of the barrier material, so that it may 
be calculated once and for all for the relevant parameters 7, a, and b. 


(b) Since the angular distribution, g°(7, cos 6), of unscattered radiation from a plane 
: : , o : : 
isotropic source as expressed in (9) tends toward = (1/cos 6) as r becomes small, representation 
Tv 


as a Legendre expansion in positive powers of cos @ becomes impractical. This difficulty 
usually does not arise for the scattered component. 
5.1. r21: Legendre Expansion 
In order to examine the behavior discussed in (b), above, g,°(7)’s were evaluated by sub- 
stituting (9) into (5). The result 


0 > ~ 
9°(7)=7- {P,Ej41(7)}> (25) 


where /;(7) is the exponential integral [11] 


es) _= 
E(x)= | e~™*u-/du, 


1 


(26) 


and P,{ £;,;(7)} is obtained by replacing (cos @)/ in P, (cos 6) by E;,, (7). Since these g,? (7), 
although of fairly general interest, do not seem to appear explicitly in the literature, a tabulation 
. oe . 
of the quantity — (/+ 3)g,°(7), (.e., (7+4) P, {£j.1(7)}), for 0 <1 < 13 and 0.001 <7 < 100 is 
o 

presented in table 2. 

A check on the values in table 2 is shown in figure 6, in which (47/c) e7 g? (7, cos 8), recon- 
structed according to (4) for 0</<7 and 0.2 <r< 10, are indicated by the dotted curves. 
Comparison may now be made with similar results from the exact expression (9), shown by the 


TABLE 2. 


of unscattered radiation from a plane isotropic source 


Coefficients (4m/o) (1+ 3) g°(r) of a Legendre expansion (4) of the angular distribution (41/0) g°(7, cos 6) 


The factor 4x/o implies unit source strength per steradian per unit of source area, and (/+12) is included to allow cumulative multiplication 
directly with the p; (a,b)’s (see table 1) according to (6). 
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ho 3 0.001 0.01 0.1 | 0.2 | 0.5 1.0 
l } } 
—_————_—____—_—_—_— jee 
0 3.1943 (0) 2.0190 (0) 9.1145 (—1) 6.1135 (—1) | 2.7989 (—1) 1.0970 (—1) 
l 1. 4889 1. 4245 | 1.0838 (0) | 8. 6130 | 4. 8996 2. 2275 
2 —6,. 1143 —3. 2090 —7.1755 (—1) | —2. 0852 | 1. 3130 1. 3712 
3 2. —2. 1124 —1.2757 (0) —8. 3188 | —2. 6901 —2. 6555 (—2) 
4 \ & 3. 3973 3.6314 (—1) —9.1715 (—2) —2. 1534 —9. 3776 
5 2. 2. 5450 1.2524 (0) 6.7606 (—1) | 1. 0298 —3. 3952 
6 8. —3. 4034 —9. 6551 (—2) 2. 4954 | 1. 9507 4. 0965 
7 3.4 —2. 8604 —1. 1624 (0) } —5. 1631 —8. 9640 (—3) 3. 9813 
8 3. 3377 —9. 7963 (—2) —3. 2542 —1.4873 (—1) —9. 6373 (—3) 
9 3. 1040 1.0492 (0) 3. 7667 | —3. 8228 (—2) —3.0111 (—2) 
10 —3. 2372 2.3859 (—1) 3. 5348 | 1.0293 (—1) —4. 9996 (—3) 
11 —3. 2983 —9. 3042 —2.6113 5.7819 (—2) 1, 8731 (—2) 
12 3. 1186 —3. 3868 —3. 5306 | 6.4944 1.0139 
13 3. 4559 8. 1428 3. 3172 —6. 1991 —9. 7479 (—3) 
\or 2.0 5.0 10.0 15.0 20.0 100.0 
:s. 
0 2.4451 (—2) 5.7415 (—4) 2.0785 (—6) 9. 5930 (—9) 4.9178 (—11)_ | 1.8418 (—46) 
1 5. 6301 1.4947 (—3) 5. 7483 2.7162 (—8) 1.4107 (—10) | 5. 4717 
2 5. 1875 1. 8564 8. 1118 | 4. 0293 2. 1490 | 8. 9433 
3 2. 1895 1. 6196 8. 8022 4. 7268 2. 6263 | 1.2159 (—45) 
4 —6.1965 (—3) 1.0195 7. 9587 4. 7804 2.8115 1. 5034 
5 —1.4221 (—2) 3.9434 (—4) 6. 1188 4. 2970 2. 7243 } 1. 7498 
6 —5.7195 (—3) —1. 1402 (—5) 3. 9668 3. 4702 2. 4242 1. 9500 
7 3. 9476 —1. 4339 (—4) 2. 0817 2. 5187 1. 9934 | 2. 1009 
8 5. 1606 —9. 8609 (—5) 9.1367 (—7) 1. 6300 1. 5167 2. 2013 
9 5.1102 (—4) —1. 3120 8.0998 (—8) 9.1882 (—9) 1. 0653 2. 2522 
10 —2. 6394 3. 0530 —1, 5133 (—7) 4. 3555 6.8479 (—11) 2. 2561 
1] —1. 5451 2. 4977 —1. 4134 1. 4721 4.0174 2. 2172 
12 8.7741 3.8280 (—6) —5.0961 (—8) 1.6851 (—10) 2.0205 | 2. 1418 
13 1.2960 (—3) —9. 5729 (—6) 4.5106 (—9) —2. 5223 (—11) 9.8335 (—12) | 2. 0335 
| 
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Figure 6. The angular distribution of unscat- 
tered radiation from a plane isotropic source 
reconstructed according to (4) and (25) (dotted 
lines) as compared with the exact expression 
(9) (solid lines) for 0.2<7r<10. 


(4 r/c) e'g%t, cos 8) 








0 a a 6 8 1.0 
cos 6 


solid curves. From both table 2 and figure 6 it can be seen that the Legendre expansion 
approach, using tables 1 and 2, should give reliable results for unscattered_as well as scattered 
radiation in the region 17315. 


5.2. r=1: Power Series Expansion 


The parameter 7 may in some cases be very small, and for this reason it seemed desirable 
to find a method for evaluating J°(7,a,5), ie., 


I°(7,a,b) =~ (dS/r*) exp (—7/cos 8), (27) 
4 Js 
for small 7. No success was had in obtaining a closed-form result, but q, (a,b) coefficients of 
the 7” power series in (10) were evaluated (see appendix B) which provide adequate convergence 
of (10) from 7=0 up to values of 7 for which the g? (a,b)’s converge rapidly. 
The first five q, (a,b) terms are 


: h —= ] a i+1 a?+1 a?+] i-l 92( y!)? a+] jv 
0 ,t ) “§ — > a ( 1 s t vy ; j " a (. n ) : 
q aie 2 inh~‘a a 23 21+-1\a’?+1 b ~ b j=0 (2) T 1)! a? 6?+1 


j ~. 


(28) 
in which the inner sum over 7 vanishes for 7=0, 
5 a eat b ab 
(a,b) =— < b sinh”! +a sinh! —tan™! ! (29) 
\ b?4 ] \ a*+-] \ a’ t 624 ] 
h 
q.(a,b) =<» (30) 
(a,b), 1 t bva?+b+1 
q3(a,b)- eae > om b(6?+-1) sinh! ms +a(a?+1) sinh7! re. = ; : }, 
oO. o 2-3 \ h2 ] \ a + ] 3 
(31) 
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and 


ab . ab’ 


l 
wad —=7 { CDuaD+5 +73 


\ 


For large a and 6, (a<6), the expression within the braces in ( 


quickly with increasing 7. 


over 7 from 2 


y. 


by noting that the series part of the expression vanishes for large 6, that is, 


lim q(a,6)=; sinh? a, 


(32) 


8) becomes small very 
Indeed, over the entire range of a and b, i.e.,0<a<b<o, summation 
0) to <=6 gives values for q (a,b) departing less than 0.01 percent from their 


respective values when the summation is carried to 7>20. <A partial check of (28) is obtained 


(33) 


in agreement with the result of Andrews [12] for the detector response to an infinite strip source 
in a nonattenuating medium. 


The remaining q, (a,b) terms form two distinct sets, 


whether n is even or odd. The general terms for these 


appendix B. 


sets 


as did the p,; (a,b)’s, depending on 
are given in (57) and (58) of 


Values of J? (7,a,b)4x/o for 0.01 <7< 10 and 0.1 <a<b< 20, as evaluated with the help 
of an IBM 704 computer, are presented in table 4. Numbers preceded by asterisks were 
evaluated using the Legendre expansion approach (15) and (25) for 0</< 13, and the others 
using the 7” series (10) for 0<n< 20. 

Since evaluation of (10) for small specific values of 7 may be faster than interpolation from 
table 4, q, (a,b) is given in table 3 for 0 <n <9 and 0.1 <a<b< 20. 


TABLE 3. Coefficients qy(a,b) of the r® power series expansion (10) of the corner position flux I° (7,a,b) of un- 


scattered radiation from a rectangular plane isotropic source 











n b=0.1 b=0.2 b=0.2 b=0.5 b=0.5 b=0.5 
a=0.1 a=0.1 a=0 a=0.1 a=0.2 a=0.5 

0 } 9.9340 (—3) | 1.9675 (—2) 3. 8971 (—2) 4.6222 (—2) 9.1603 (—2) | 2.1600 (—1) 

| —9. 9669 | —1. 9836 —3. 9481 —4. 8047 —9. 5656 —2. 3215 

2 5. 0000 10.0000 (—3) 2. 0000 2. 5000 5. 0000 1. 2500 

3 —1. 6722 —3. 3609 —6. 7547 (—3) —8. 6819 (—3) —1. 7443 —4, 4960 (—2) 

4 4.1944 (—4) 8.4722 (—4) 1.7111 2. 2639 4. 5694 (—3) 1. 2153 

A —8. 4169 (—5) —1. 7086 —3. 4679 (—4) —4. 7282 (—4) —9. 5873 (—4) —2. 6333 (—3) 

6 1.4075 | 2.8716 (—5) 5. 8574 (—5) 8. 2390 (—5) 1. 6783 4. 7647 (—4) 

7 —2.0174 (—6) | —4. 1368 (—6) —8. 4805 (—6) —1. 2320 —2.5213 (—5) —7.4048 (—5) 

8 2. 5302 (—7) | 5.2148 (—7) 1.0744 1.6141 (—6) 3. 3183 (—6) | 1. 0090 

9 —2.8208 (—8) | —5. 8435 (—8) —1. 2101 (—7) —1. 8819 (—7) —3. 8869 (—7) | —1. 2248 (—6) 

n b=1.0 | b=1.0 b=1.0 b b=2.0 b=2.0 

} a=0.1 } a=0.2 a=0.5 a=1 a=0.1 | a=0.2 

0 7. 8327 (—2) } 1.5540 (—1) 3. 6887 (—1) 6. 3951 (—1) 1.1046 (—1) | 2.1946 (—1) 

1 —8. 8020 } —1. 7535 —4. 2718 —7. 9336 —1. 4421 —2. 8755 

2 5. 0000 10.0000 (—2) 2. 5000 5. 0000 10.0000 (—2) 2. 0000 

3 —1. 9154 —3. 8455 —9. 8624 (—2) —2. 1346 —4. 9338 —9. 8916 (—2) 

1 5. 5694 (—3) 1. 1222 2. 9514 6.9444 (—2) 1. 9472 3. 9111 

) 1.3114 —2.6517 (—3) —7.1478 (—3) —1. 8357 —6. 5195 (—3) —1.3113 

iD 2. 6050 (—4) 5. 2848 (—4) 1. 4593 4.1049 (—3) 1.9117 3. 8496 (—3) 

7 ~4. 4892 (—5) —9. 1365 (—5) —2. 5829 (—4) —7.9826 (—4) —5.0012 (—4) —1. 0081 

8 6. 8492 (—6) 1, 3982 4.0448 (—5) 1. 3771 1.1811 2. 3829 (—4) 

9 ~9. 3939 (—7) —1. 9233 (—6) —5. 6906 (—6) —2.1392 (—5) —2. 5399 (—5) —5. 1290 (—5) 

n b=2.0 b=2.0 b=2.0 b=5.0 b=5.0 b=5.0 
a=0.5 a=1.0 a=2.0 a=0.1 a=0.2 a=0.5 

0 5. 25386 (—1) 9. 3072 (—1) 1.4074 (0) 1. 3708 (—1) 2. 7263 (—1) 6. 5729 (—1) 

l —7. 0456 —1.3285 (0) —2. 2916 — 2.3108 — 4.6120 —1. 1372 (0) 

2 5. 0000 1. 0000 2. 0000 2. 5000 5. 0000 1. 2500 

3 —2. 5139 —5. 3007 (—1) —1. 2410 —2. 3179 —4. 6397 —1. 1666 

4 1. 0069 2. 2222 6.1111 (—1) 1. 9451 3. 8944 9. 8090 (—1) 

5 —3. 4104 (—2) | —7. 8243 (—2) —2. 5253 —1. 4685 — 2. 9404 —7. 4123 

6 1. 0095 2. 3951 9.0494 (—2) (—2) 1. 9849 5. 0081 

7 —2. 6625 (—3) —6. 5116 (—3) —2. 8734 —1, 2004 —3. 0320 

Ss 6. 3366 (—4) 1. 5949 &. 2096 (—3) 6. 5448 (—2) 1. 6550 

Y —1. 3730 —3. 5550 (—4) —2. 1351 —3. 2405 —8. 2045 (—2) 
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TABLE 3. 


Coefficients q,(a,b) of the 7" power series expansion (10) of the corner 
scattered radiation from a rectangular plane isotropic source 










n b=5.4 b=5.0 
a=1.0 a=2.0 

0 1. 1879 (0) 1. 8793 (0) 

1 —2. 1837 — 3. 9333 

2 2. 5000 5. 0000 

3 —2. 3782 —5. 0728 

2.0139 4. 4444 

5 —1. 5272 —3. 4397 

6 1. 0353 2. 3716 

7 —6, 2923 (—1) —1, 4665 

8 3. 4493 8.1909 (—1) 

4 —1.7178 —4. 1614 

n b=10.0 b=10.0 
a= 1.0 a= 2.0 

0 1. 2849 (0) 2.0692 (0) 

l —2. 8671 — 5. 2864 

2 5. 0000 1.0000 (1) 

3 —8. 7046 —1. 7838 

1. 4444 (1) 2. 9722 

5 —2. 1687 —4. 4719 

6 2. 9038 6. 0067 

7 —3. 4758 —7. 2190 

8 3. 7462 7. 8159 

+] — 3. 6647 —7. 6824 

n b=20.0 b=20.0 
a= 0.5 a= 10 

0 7.3091 (—1) . 3345 (0) 

l 1.8253 (0) —3. 5577 

2 5. 0000 .0000 (1) 

3 —1. 6854 (1) —3. 3781 

4 5. 6007 1.1222 (2) 

5 —1. 6803 (2) —3. 3669 

6 4. 4847 8. GRS2 

7 —1. 0690 (3 —2. 1430 (3 

S 2. 2935 4. 5990 

9 —4. 4652 —. 9564 


TABLE 4. 


Unscattered flux (4x/a)1°(7,a,b) of ra 
source calculated using tables 1 and 2 (numbers preceded by « 


b=0.1 6=0.2 
a=0.1 a=0.1 
9.93 (—3) 1.97 (—2) 
9.83 1.95 
9. 74 1.93 
9. 44 1. 87 
&. 99 1. 78 
8.13 1. 61 
*8. 48 *1.65 
6. 02 1.19 
*5. 63 "1.15 
3. 64 7.18 (—3) 
*3. 66 *7. 02 
1. 34 2. 62 
*1.37 *2. 67 
6.60 (—5 1.28 (—4) 
*6. 59 of 
*4. 37 7 *8.26 (—7) 
b=1.0 b=1.0 
a=0.1 a=(0).2 
7.83 (—2 1.55 (—1) 
7.75 1. 54 
7. 66 1. 52 
7. 40 1. 47 
7. 00 1. 39 
6. 26 . 24 
*6, 28 *1, 24 
47 8.85 (—2 
*4.45 *8. 93 
2. 56 5. 06 
*2. 57 *5. 06 
8.49 (—3) 1. 67 
*8.53 MY 
85 (—4 7.55 (—4) 
*3. 29 *6. 35 
*1.67 (—6 *3.16 (—6 
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eT et bt et et et et 


| 
Noh 


—f}, 


* 
to th 


* 
tot 


. 1467 


b=5.0 
a=5.0 


. 7262 (0) 


3832 

2500 (1) 
6595 
8403 
7554 
4695 
0971 
3983 
5541 





6=10.0 
a= 5.0 


(Q)) 
0571 (1) 
5000 


. 0356 


R889 
3854 
9201 
3848 
6756 


7323 


(2) 


b=20.0 


a= 2.0 


1679 (0) 
6641 

0000 (1 
S107 

2611 (2 
7850 
8127 
3264 
2953 


8123 (4 


(3) 


diation, at a corn 


b=0.2 


a=().2 


vO (—% 
86H 
UD 
70 


18 


21 


69 (—] 
65 
60 
29 
93 
¥3 


COON 


rho 


9 2 


5. 61 2 


b=10.0 
a 0.1 


4685 (—1) 
9966 
0000 
6256 


4307 (0) 


b=10.0 
a=10.0 


7961 
6127 
0000 
2892 (2) 
8194 
3227 
R285 


(0) 
(1) 


Continued 


position flux I° (r,a,b) of un- 








b=10.0 b=10.0 
a= 0.2 a= 0.5 
2.9217 (—1) 7.0607 (—1) 

— 5. 9833 —1. 4798 (0) 

10. 0000 2. 5000 

—1.7256 (0) — 4, 3227 

2. 8622 7.1701 

—4, 2957 —1.0762 (1) 

5. 7461 1, 4399 
—6. SO8Y —1.7217 
7. 3929 1. 8536 

—7. 2214 —1. 8112 
b=20.0 b=20.0 
a= 0.1 a= 0.2 

1.5182 (—1) 3.0211 (—1) 

— 3. 6878 —7. 3658 

10. 0000 (0) 

—3. 3683 (0) 

1.1195 (1) (1) 

—3. 3586 

&. 9635 (2) 

—2. 1364 (2) 

4. 5831 Y. 

—8. 9221 —1. 7846 (3) 
b=20.0 b=20.0 
a=10.0 a=20.0 

4.2228 (0) 4. 8801 (0) 
—2. 2546 (1) —3. 3719 (1) 
1.0000 (2) 2. 0000 (2) 
—3. 9741 —1.0231 (3) 
1.3972 (3) 4. 4611 
— 4. 3433 —1. 6801 (4) 
1. 2006 (4) 5. 5605 
—2. 9741 —1. 6422 (5) 
6. 6542 4. 3835 
—1.3550 (5) —1. 0685 (6) 


er position, from a rectangular plane isotropic 


usterisks) and using table 3 (no asterisks) 


6=0.5 
a=().1 


4() ] 
32 
01 


65 


uy 
OO 


62 


99 3 
6Y 








h=0.5 b=0.5 
a=().2 a=0.5 
9.16 (—2) 2.16 (—1) 
9. 07 2.14 
8.97 2.11 
8. 69 2.05 
8. 25 1.94 
7.43 1.74 
*7. 44 *1. 75 
5. 44 1, 26 
*5. 44 a 
3. 23 7.38 (-—2 
*3. 23 *7. 38 
1.14 2. 53 
*1.14 *2. 54 
5.04 (—4) 1.04 (—3) 
*5.02 *1.03 
*2.83 (—6) *5.15 (—6 
b=2.0 b=2.0 
a=().1 a=0.2 
1.10 (—1) 2.19 (—1 
1.09 2.17 
1,08 2.14 
1.04 2. 06 
9.70 (—2) 1.93 
8.53 1. 69 
*2 53 *1. 69 
5. 82 1.15 
5 7O *1. 15 
3.14 6.20 (—2 
*3.14 *6. 21 
9.57 (—3) 1. 88 
*9 60 *1. 88 
*3.37 (—4) *6.54 (—4) 
*1.67 — *3.16 6 


pati att 2 isotropic 
- 2 4. Unscattered flux (4x/c)1°(7,a,b) of radiation, at a corner position, a . nt te eas eee 
i ieee using tables 1 and 2 (numbers preceded by asterisks) and using ta , e: ‘ 
8 e ies ae 



































| uf b=5.0 
pee b=2.0 b=2.0 b=5.0 can a=0.5 
=0.5 a=1.0 | a=2.0 stiles ; 
r a=0.5 | ° ics | ine 3 Sears ae With aaeeeniintiniaiiginiaiial 
F i 3 _ | 37 (— 2.73 (—1) 6.57 (—1) 
9.31 (—1) | 1.41 (0) 1.37 (-1) | $4 6. 46 
° 9.18 1. 38 * 2 64 6. 35 
0.01 ee 1. 36 1. 33 rage 5. OF 
‘ 9.05 on | > 5 6. 03 
0.02 8. 67 | 1.30 1. 26 2 31 5. 55 
0.05 8.07 | 1. 20 1.16 « 
0.1 } = | e 4.72 
= 1.02 9.93 (—2) 1.97 Pept 
a 4.02 m@ | *1.02 | 9.97 | %1.98 4.73 
— | we? 6 38 | 1. 27 3.00 
0.5 2.72 4. 62 . 38 (—1) +6, = *1. 6 *2.99 
"3:73 *4. 62 6.39 oa 
3.27 3.46 (—2) . 
9 2 3.02 } 3. 27 0. 20 #150 
1. 44 2. 35 > me } *3 97 *6. 45 -« 
ate *1.44 *2. 35 *3. 02 
> « ‘ id 6 —2) | an = 
2.0 38 (3) 43° oF 45 | | —*9.69 (—3) | *1.90 | *4. 28 (—2) 
*4 4 . Be ? | “ = = 
- 2 oF *6§. 54 (— 1.35 (—3) 
co *1.32 (—3 *1.75 (—3) *1.80 (—3) *3.37 (—4) 6.54 (—4) 
1.32 (—3) | heen *3.16 (—6) | *5.79 (—6) 
10.0 *5.79 (—6) *6.51 (—6) *6. 52 (—6) | 1.67 (—6) | 3. wie i ESE 
ee oe: eum } os b=10.0 b=10.0 
~ b=5.0 | b=10.0 ey | a=0.5 
yan 4 a=5.0 a=0.1 a=0.2 ih 
Tr a=1. ; Bese Seca se ate he — — pinata —— 
aa oT a 2.92 (—1) 7.06 (—1) 
19 (0) 1.88 (0) 2.73 (0) ar {1 | 2. 86 ¢ | 6. 92 
. ee 1. 84 2. 65 1. 44 2 81 6.77 
ee vis 1. 80 2. 58 | eo 2 65 | 6. 38 
0.02 1.15 169 2. 39 | 1.33 2. 65 a 
0.05 1.08 : 53 2 10 1. 21 2. 41 “eo 
0.1 9.92 (—1) oe | sai | 43 
me on 1. 64 1.02 a Ve #4 8 
0.2 Rg a “165 | *1.02 |  *2.02 4. 84 
oO —_ | | é 
> 49 (— | 1. 27 | 3. 01 
797 (— 8.53 (—1) 6.42 (—2) a + 
od 314 “a” “349 | *6. 39 *1. 23 3.01 
*5.1 7. 27 . 
- 3 3.44 hae "1 5 
1.0 2. 47 x = *2 44 *3. 27 *6.47 (—2) 1. 50 
*2.47 se i *4 9 il 
2.0 *6.51 (—2) *7,57 (—2) *7.69 (—2) *9.69 (—3) *1..90 4.28 (—2) 
sess ings . *6.54 (— *1.35 (—3) 
5. 0 ; * _ | *3 37 (—4) 6.54 (—4) 
*1.75 (—3 *1.80 (—3) 1.80 (—3) | exaiaie 
| . * =. 219 (— 
ne *6.51 (—6) *6.53 (—6) *6.54 (—6) | "1.67 (—6) a ‘- — 
: . _ aon ae a b=20.0 b=20.0 
a b=10.0 b=10.0 “ a=0.2 
b=10.0 b=10.0 =o a=10.0 a=0.1 . 
* a=1.0 a=2.0 a=95.0 . ae ba Pa ie 
: : 3 eeeed aa m6 3.02 (—1 
2.07 (0) 3.15 (0) | 3.80 (0) 1.52 (—1) a (-)) 
0 1.28 (0) “. 04 - 3. 64 1. 48 “9 
: 9 (2 3. 04 2. 88 
0.01 1. 26 197 2.94 3.49 1. 45 a 
oe ci 1.83 2. 67 3.10 a 2. 43 
0.05 1.15 es 2 30 2. 58 1. 22 7 
0.1 1. 04 7 O& ee | 2 03 
- 1. 86 . 02 abe 
0.2 8.60 (—1) Re “173 *1.83 *1. 06 *2. 08 
8. 58 ig | | 
. eo | 8.78 (—1) ; : | on 
0.5 18 os ae (3) p. a (-D *8 77 *6.39 (—2) *1. 27 
*5.18 7.3: . ies 
1.0 *3 9] *3. 44 *3.45 *3. 27 *6. 47 (—2) 
? *2 47 3.2 ‘ s 
| ati . 9 (—% 1.90 
20 *6.51 (—2) "7.56 (—2) *7.69 (—2) | *7. 68 (—2) aaah Bi. 
51 (—2 dee 4 on 
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6. Examples 


Comparison of the results in the above sections 4 and 5 is made in table 5 with similar 
results of Berger and Lamkin [1] for a water barrier of 7=4 and a=b=2, using the compound- 
ing rule in section 3 for noncorner positions. Agreement is within 5 percent. 

There exists little experimental information on flux fields from rectangular gamma-ray 
sources. However, an additional partial check on qo(a,b) was supplied by Auxier et al., [3] 
as shown by the circles, triangles, and squares in figure 7. The distributed source was rep- 
resented by an array of matched (+5%) Co™ sources spaced 2 ft apart over a 24- by 32-ft ree- 
tangle. The detectors were pocket ionization chambers with an average response of 14.4 milli- 
roentgens per hour at 1 ft from a 1-millicurie point Co” source. This response calibration 
was made in a geometry such that the combined effects of backscattering from the ground, 
the buildup factor in air, and the exponential attenuation by air were roughly equivalent to 


TABLE 5. Dose rate, in milliroentgens per hour per milli- 150 
curies per square foot of source strength, at a number of 
positions a=x/h, 8=y/h at distance h from a rectangular 























surface of relative dimensions W/h=2,1/h=2 separated from | = 
a plane source of 1-Mev photons by a layer of water of = 
thickness 7=4 e 
€ 
The numbers in parenthesis* are from table 1 of Berger and Lamkin [1]. = \ 
— - 2 
\Qe 0 0.25 0.50 0.75 1.0 4 
&\ | : 
Boe i|- 
| - 
0 0.211 a 
(. 22) : 
25 275 0. 361 -: 
(. 28) (. 37) = 
— 
. 50 .323 425 0. 503 
(.33) . 44) (. 51 
75 . 350 462 . 546 0. 593 : . 
(. 35) (. 48) (. 55) (. 60) re) 
) a 8 12 16 20 
1.0 . 359 . 473 . 559 . 608 0. 624 FEET FROM CENTER, ALONG DIAGONAL 
(. 37) (. 49) (. 57) (. 62) (. 63) 
Figure 7. Detector response in milliroentgens per hour ath=8, 
* Arithmetic error in the Berger-Lamkin result at a=0.25, 8=0.25 is here 0 and 11 jt above the diagonal of a rectangle of dimensions 
corrected. w=24 ft, l=32 ft covered with one millicurie per foot [2] of 
Co”. The circles, triangles, and squares are experimental 
values of Auxier et al. [3] and the curves are results from (28) 
compounded according to (11). 


these combined effects in the actual measurements. Thus, since the experimental data in 
figure 7 are related to a source strength of 1 millicurie per square foot, the dimensional factor 
o/4mr in (10) may be identified as 14.4 milliroentgens per square foot, and we may multiply 
qo(a,6; a8) by this number for an absolute comparison. The curves in figure 7 represent 
14.4-q(a,6; a,8), and agree with the experimental data to within +10 percent despite neglect 
of attenuation and scattering effects. 


The authors thank M. J. Berger and L. V. Spencer for their interest, encouragement, and 
helpful criticism regarding this effort. 
Appendix A. Evaluation of the p, (a, b) Coefficients 
The Legendre polynomials, P;(cos @), may be written in the form [13] 


Pianist (—1)"(20—2))! . sa 34) 
1(COS 6 = 25 SAT — jy —2j 08 ’ (54, 


where m=//2 if l is even, and m=(/—1)/2 if / is odd. Substituting (34) into (15) and using 
the notation 
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nv n! oe 
( r a SY (35) 


for the 7’th coefficient of an n’th degree binomial expansion, we have for the even p,’s 


te £ ‘ > Qs "a *b 
wen (—1)7/21—7\ (41—2) dadB ries 
~ p> oo ( j )( 2l—7 )I. + op + ] )! I+3 2" (36) 


Integration over @ yields[14] 


bj (—1)'4 (21—j\ (41-—29\ (1) .. *® dg , 
. my x 2i+1 +“ 37 
sia 24 oy ps "(27-44 r j )( 2l—) )( i )a i (6?+ 1)'—/*1(q?+ B+ 1)'*? mr 


For /=0, integration over 8, using eq (234.19a,, reference [14], vields the result in (16). 


For the higher.even terms, use of eq (234.18a) leads to (18) and higher terms which may be 
expressed in the general form 


\a-peme oer oe l+k 
7 i—k—1 92k 
: ab et ” ak JA a) 
P= “UC: 21-4 921 > b2 a a 


- ya’ 1 k=0 >) ae 
(7, ) 2k+0) 
: 1 1 
9] 9p oo ed ree ; 3 
XX (2/ aI T 1) 2 92) (a? aay Ganapati 7) (38) 


The odd p,’s as obtained from (15) and (34) may be expressed as 


‘ee bi-j petites 5 * da dB ne 
eis >? ar . S j Core —j ) ea (a? +6? 1)"-3#? _ 


Integration over a yields (see ref. 14 eqs 15, 15a, b) 


rope BOY GT DC Le) 


— jJ=0 


l > tan7! (a/y6?+1)dB 
x 931 3 | ra (B2+-1) '-4+3/2 : 
+a s Ta ~ ‘x dp eee (40) 
nen 2 ¥ nadia ay I—j—i+1)J po (B+1)'*! (2? +641) 4-1 
—  l-j-4 


The integral containing the arctangent was amenable to integration by parts, and the purely 
algebraic integrand was reduced to terms of standard form by successive application of the 
identity 


] ] l ] ; 
(6?+1) (a2+6?+1) @ ign ae l ) (41) 


leading to the _— form for the odd terms, including (17) and (19), 
L\ (L+R+1)\ (214+-2k4+1 
-(,)¢ k: )( l+k ) 
1 ; b Pee 1 a oe aa b 
x) ap Ve ei @ + vet Oe” Vet 


2, 2 l l q 
i ~ jC) (?+ 64 WACC ESD a 34 Toy 1)*-74 :)} (42) 


I bom 


1)" 
Pritt 22° *(2i-+-1) f=0 Qk 








Appendix B. Evaluation of the q,, (a, b) Coefficients 


The integral (27) may be written in cartesian coordinates as 


*a *b ; 
[°(7,a,b) =7 | da dB exp (—rya?+ B+ 1)/(a?+B?--1). (43) 
The substitution 
a= 768°+1 tan 7 (44) 
vields 
, ] tan-! : 
of i we. B ¥ vet 5 é 
I (7,a,6)=7- , Be 7 Je exp (—7y6’+1 sec n)dn. (45) 


Expanding the exponential as a power series of its argument we have 


a 
tan-! 


it a i . 
0 r .b) ve o : ( | .s 7 : 2] sec n (46) 
I°(7,a,6 i. Fri Jo dot ( \ B°+-1 sec n) dn. 


The first integration is then simply over powers of sec n, which can be found in integral tables 
as standard forms [15], giving 


] i o lf l a -ah-! a 
77.8. 0)== = tan ; — 7rsinh-’ |; 
4 Jo ~6°+-1 y8°+1 Ve?+1 


Tr? 3 21ALl 1 21 1 E a 7 3 
+— a—> | « = a rele sinh7! VB qa E +a(p°-4 ) |-. ‘s Sus (47) 


Integration term by term over 6 then yields a power series in 7 of the desired form (10), in 
which the coefficients q,(a,b) are 


qo(a, )- [. re tan! rz dp, (48) 
‘ - 

qi(a, b): —{. sinh” VB x dp, (49) 

g(a, b) od | as, etc., (50) 


as evaluated using (47). 
No success was had in evaluating qo(a,b) in closed form. However, use of the arctangent 


expansion [16] 
D2m / m 1)? \mt+l 


x tan“'z > ( --;) (51) 


m=0 (2m T 1)! 


yields the series of integrable terms 
*db 


l 2. 2*(m!)? a’ \" 
o(a, b) =- > > ——— ) 1B. 52 
Go oJ), Gen (aaa ol (52) 


Integration of these terms (see ref 15, p. 25) gives the result 


b Ve+1 < l a® \{t! 
qo(a,6)=tan"! --—— 9 = ( - ) 
instil ya?+1 a a 21+ 1 \a*+1 


io. 14 7 & ¥" El oe 7 es). Ye 
+7 Pe or l ey i) f=0 (27-44 Wee i) 9) 


“ j=0 
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in which the double summation vanishes for 7=0. Since convergence of the summations in 
(53) is poor when both a and 6 are large as compared with unity, it is advantageous to rewrite 
the first member of the right-hand side, using the substitutions 


b oo e+ 
tan7! - —=-—tan-! p dl (54) 
Va?+1 2 b 
and 
so. 4 2 \itl 
vya°+1 = 1 a , a a ae 
» sal Ts a7 (= —tanh-! ——— =sinh la. (55) 
a@ f= 2.+1\a@+1 yar+1 


This gives 





_= _, Va?+1 Va?+1 oo 1 a? \¢tt! 
qo(@, 6) “5 sinh-! a—tan ot uae — acy ss 


— i=0 4 
b= 1 a2 \ittl il 274(9!)? a?+1 i - 
ta ey 41 (or 24 (7+ 1)! aE bEPI (56) 


which, after factoring the quantity 


-. 1 oe \fe 
5 
@ j=0 21+1 a’+1, 


from the last two members, may be identified with the result in eq (28). Now, since the inner 
summation over 7 in (28) consists of the leading terms of the expansion of (, a?+1/b) 
<tan7! (Va?+1/b) according to (51), the quantity within the braces in (28) is the remainder after 
i—1 terms of this expansion, causing the outer summation over 7 to converge very quickly for 
all values of a and 6. 

The remaining g, (a,b) coefficients for n even and odd, respectively, may be written in 
the general forms 


1 n—2 7 co a n—1 n—1 q2"—2i-1p2i+l 
Jon (a, 5) “ayy (—D" 2 (—1)"( ; )(2i-+2)arrsa(a, 6) (" ; yen 





i=l / —271—1)(27+1) 
(57 
and 
ee ee a 
i i (2n+1)! 2n+1 sei iiaiall: 
2n! ; a ; b 
tae in | 0(b?+1)* sinh - —-+a(a?+1)" sinh7! ——— 
2°" (n!)*(2n+-1) Ve+1 Vae+1 
se Bo! aN)? OW ae pate . ss 
+abya’?+6?+1 30 Qi+1)! (a?+-b?+-1)*{(a?+-1)"-*-?4 (674 pny |} (58) 
izo (20-+1)! 


from which q,(a,b)’s for 1 < n < 4 are given in (29 to 32). 


References 


[1] M. J. Berger and J. C. Lamkin, Sample calculations of gamma-ray penetration into shelters: contributions 
of sky shine and roof contamination, J. Research NBS 60, 109 (1958) RP2827. 

[2] L. V. Spencer, Structure shielding against fallout radiation from nuclear weapons (to be published as an 
NBS Monograph). 

[3] J. A. Auxier, J. O. Buchanan, C. Eisenhauer, and H. E. Menker, Experimental evaluation of the radia- 
tion protection afforded by residential structures against distributed sources, CEX 58.1, ORNL (Jan. 19, 
1959). 


137 








(7| 
[8] 


[9] 


U. Fano, L. V., Spencer, and M. J. Berger, Penetration and diffusion of X-rays, Encyclopedia of Physics 
(Springer-Verlag, 38 II, 1959). 

R. A. Herman, A treatise on geometrical optics, p. 210 (Cambridge Univ. Press, London, 1900). 

Z. Yamanouti, Geometrical calculation of illumination due to light from luminous sources of simple forms. 
Researches Electrotech. Lab. (Tokyo) No. 148 (1924). 

P. Moon and D. E. Spencer, Light distribution from rectangular sources, J. Franklin Inst. 241, 195 (1946). 

A. I. Mahan and W. F. Malmborg, Radiation characteristics of semicircular, circular and rectangular 
surface sources, J. Optical Soc. Amer. 44, 644 (1954). 

QO. Seiber, Die Wirmeaufnahme der bestrahlten Kesselheizfliche, Arch. Warmwirtsch. 9, 180 (1928). 


[10] See, for example, W. B. Boast, Illuminating engineering, Ist ed., p. 89 (MeGraw-Hill Book Co., New 
York, N.Y., 1942). 

[11] G. Plaezek, The function E,,(x _ run—"du, MT-1, NRC No. 1547, Chalk River, Ont. (1946). 

[12] D. G. Andrews, Gamma ray shielding for engineers, U.K.A.E.A., Report Risley 5225 (Sept. 1954). 


[13] 
[14] 
[15] 


[16] 


) 

R. V. Churchill, Fourier series and boundary value problems, Ist ed., p. 177 (McGraw-Hill Book Co., 
New York, N.Y., 1941). 

W. Grébner and N. Hofreiter, Integraltafel I, Eqn. 234.17b. (Unbestimmte Integrale), (Springer-Verlag, 
1957). 

H. B. Dwight, Tables of Integrals and Other Mathematical Data, 3d ed., p. 97 (MeGraw-Hill Book Co., 
New York, N.Y., 1957). 

Kk. P. Adams and R. L. Hippisley, Smithsonian mathematical formulae and tables of elliptic functions, 
p. 122 (Smithsonian Institution, Washington, D.C., 1922). 


Wasuineton, D. C. (Paper 64C 2-32) 


138 








JOURNAL OF RESEARCH ot the National Bureau of Standards—C. Engineering and Instrumentation 
Vol. 64C, No. 2, April-June 1960 


Microwave Attenuation Measurements With 
Accuracies From 0.0001 to 0.06 Decibel 
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The application of certain power stabilization and measurement techniques to the 
problem of attenuation measurement has yielded a measurement system with a stability 
and resolution of the order of 0.0001 decibel. A practical application for this technique was 
recently provided in the calibration of a rotary vane type of variable microwave attenuator. 

In order to take complete advantage of this increased stability it was necessary to apply 
refined techniques to the evaluation and reduction of mismatch error. 

This proved to be by far the most exacting practical application of the cited techniques 
encountered to date, but the results of this calibration showed excellent agreement with the 
mathematically predicted values used in marking the attenuator dial. 


1. Introduction 


The following is a report on some recent work at the Radio Standards Laboratories of the 
National Bureau of Standards at Boulder, Colo. It is believed that this work represents per- 
haps the most accurate measurements of microwave attenuation yet made. 

In the measurement of microwave attenuation, the stability or resolution of the associated 
measuring system places a limit on the accuracy with which small values of attentuation may 
be measured, because the reduction of other sources of error (such as mismatch) will tend to 
make the total error approach this value. In practice this limit has been typically of the order 
of +0.01 db. 

The application of certain power stabilization! and measurement? techniques to the 
problem of attenuation measurement has yielded a measurement system with a nominal 
hundredfold increase in resolution and stability. A practical application for demonstrating 
the capabilities of this technique was recently provided in the calibration of a rotary vane type 
of variable microwave attenuator which is characterized, in part, by a high degree of resolution 
for small values of attenuation. Thus when a recent application called for this characteristic 
in a variable attenuator, and when a preliminary study of the mismatch and stability character- 
istics of an attenuator of this type indicated that such a step was warranted, it proved desirable 
to experimentally confirm the mathematically predicted values of attenuation which were 
engraved on the attenuator dial. This proved to be by far the most exacting practical appli- 
cation of the cited techniques encountered to date. 

In order to take complete advantage of the improved stability, refined techniques were 
applied to the evaluation and reduction of the mismatch error. The capability of the system 
is indicated by the tabulated results of the attenuation calibration, and the estimate of the 
limits of error is supported by an analytic and experimental treatment. 





! Glenn F, Engen, Amplitude stabilization of a microwave signal source, IRE Trans. on Microwave Theory & Tech. MTT-6, No. 2, 202-206 
(Apr. 1958). 


2 Glenn F. Engen, A self-balancing direct-current bridge for accurate bolometric power measurements, J. Research NBS 59, No. 2, 101-105 
(Aug. 1957) R P2776. 
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2. The Calibration System 


A simplified diagram of the calibration system is shown in figure 1. The attenuator under 
test is placed between an amplitude stabilized microwave signal source! and a bolometer 
mount and power meter. The power meter consists of a self-balancing d-c bolometer bridge 2 
with provisions for measuring and recording the d-c bias power required to maintain the bolom- 
eter at its operating resistance of 200 ohms. The two bolometer mounts M, and M, shown 
in the temperature stabilized water bath are for power measurement and signal source stabiliza- 
tion, respectively. A reasonable amount of care was exercised to obtain good performance 
from each item of equipment, with a resultant system performance as shown in figures 2a and 
2b. Figure 2a shows the stability and repeatability with the attenuator alternately set at 
the 0.00- and 0.01-db positions. It will be noted that the stability and repeatability are better 
than 0.0001 db (10 microbels). A recording of the long term system stability is given in figure 
2b where the maximum variation is of the order of +10 microbels. 
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Figure 1. Simplified diagram of calibration system. 


The absolute accuracy of measurement of small changes in d-c power at this level is 
estimated to be of the order of 0.02 ww. This indicates that further improvement could be 
expected in the results if the system stability could be improved in some manner, perhaps by 
use of a frequency-stabilized signal source. 


3. Theory of Measurement 
One could make an attenuation measurement by measuring the microwave powers P, 


and P, absorbed by the bolometer mount .M, when the attenuator is set first on zero, then to 
some other setting. The relative attenuation A at this setting is 


A=10 logo 5 (1) 
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ee FiGuRE 2a. System response to a 0. 01 db change in attenuation. 


3 HOURS 





FicurE 2b. Long-term system stability. 


In making such a measurement, one assumes that the microwave power P absorbed by the 
bolometer element is proportional *® to the d-c power W withdrawn in order to keep the bolom- 
eter resistance constant. Letting W> represent the d-c power required to bias the bolometer 
at its operating resistance when P equals zero, 


Wo—W, . 
A=10 logy WW, (2) 
where W and W, are the d-c bias powers corresponding to P; and P,. (For the bolometer used, 
W)~15mw and P, ~10mw.) 

The apparatus employed permitted direct measurements of differences in d-c power, a 
procedure permitting greater accuracy and convenience than calculation of differences from 
separate measurements. The changes in d-c power level during an attenuation measurement 
are shown in figure 3. 


The constant of proportionality is determined by the substitution error of the bolometer which is known from other experimental data to 
be independent of power level. 
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FIGURE 3. Changes in d-c power level during an attenuation measurement. 


For attenuations less than approximately 3 db, the power difference W,—W), was meas- 
ured directly. If we let W.,=W.—W), and Wo,=Wy— Wi, eq (2) may be written 


] 9 
W,, (3) 
Wo 


A=10 logo 
] 


For attenuations greater than approximately 3 db, the power difference Wo—W>. or Wo 
was measured directly. Equation (2) becomes 


A=10 logio yy (4) 


4. Propagation of Error in Measuring D-C Power Differences 


It is estimated that the error in measuring d-c power differences is within 0.1 percent 
+0.1 yw. When W,—VW, is measured, it can be shown that the limit of error in determining 
the attenuation is 


l : 
e= 10 logy : (5) 
- l 0.luw W., 
l . + 0.0901 => 
Wal Wot Wa 
Wo 
and when W,—VW,, is measured, 
10 loge . . (6) 
, E10 0.1uw+0.001 Wo 
ii Woe 


The calculated limits of error are shown in figure 4. 


5. Mismatch Errors 


The mismatch error‘ in calibrating a variable attenuator depends upon the reflections from 
the system in which the attenuator is placed and upon the changes in characteristics of the atten- 
uator as its dial is moved from the reference position. The graph of figure 4 shows calculated 


4R. W. Beatty, Mismatch errors in the measurement of ultrahigh-frequency and microwave variable attenuators, J. Research NBS 32, 
No. 1, 7-9 (1954) RP2465. 
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FicuRE 4. Limits of error as a function of attenuation. 


limits of error ® for the attenuator used, based upon measurements of the magnitude of the 
changes in the scattering coefficients S;; and S2. of the attenuator. It was assumed that there 
was no phase change in S,,.. The mismatch error is below 0.0001 db for attenuator settings up 
to 0.1 db, and remains below 0.001 db for higher settings. 


6. Results 


The calibration data taken at 9.3897 kMce is shown in table 1. Three sets of data are 
shown in order to give an idea of the resettability of the attenuator, and more significant figures 
are given than one can normally use when interpolating between the marked dial divisions. It 
should be noted that the calibration is not necessarily representative of this type of attenuator, 
because only one such attenuator was calibrated by this procedure, and this was done at only 
one frequency. It should also be noted that this type of calibration is not available on a routine 
basis, because there are as yet no commercially available attenuators which can be read and are 
repeatable to 0.0001 db. 

The estimated limit of error for the complete range of the attenuator, as determined from 
figure 4, is also shown in table 1. Above 20 db, the calibration was made in two parts: Measure- 
ment of the 20 db step, and measurement of the additional attenuation referred to this step. 
For these values, the quoted limit of error is the sum of the errors in the individual steps. The 
accuracy of setting the attenuator dial on the marks is not as good as the accuracy of the 
measurements. 

It is noted that at the low end, the estimated accuracy of the measurement is also better 
than the repeatability of setting the attenuator on the mark. One concludes that further 
improvements in attenuators need to be made if we are to take full advantage of the most 
accurate measurement techniques. 


5 See appendix. 
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TABLE 1. Results of three sets of measurements and calculated errors in single measurements 





Attenuator Measured attenuations corresponding to dial reading } Calculated limit 
dial Cee a ea a tikes. NY oo Sea a ee ___ |) Maximum devia- of error in single 
reading | | tion from average measurements 
| 1 2 3 | Avg | 
i Pe ne —— ” bc es Se 
db db db db | db db | db 
0.01 0.01077 0.01083 0.01044 | 0.0107 | 0.0003 | 0. 000055 
.02 . 02123 .02145 . 02145 . 0214 } . 0002 } . 000068 
.03 . 03015 | . 03035 . 03038 . 0303 . 0001 | . 000080 
.04 . 04043 . 04088 j . 04074 .0407 | . 0003 | . 000091 
-05 . 05181 | - 05216 . 05226 . 0521 | . 0003 . 00010 
i ] 
. 06 . 06102 . 06067 | . 06086 | .0609 . 0002 . 00011 
07 | . 06992 . 07027 . 06985 . 0700 } . 0003 | . 00013 
.08 - 07991 . 08056 . 08017 . 0802 . 0004 | . 00014 
.09 - 09075 . 09104 - 09078 . 0909 . 0001 | - 00015 
al - 10226 | - 10203 } - 10210 | -1021 | . 0002 | . 00016 
| 1 
12 - 11886 - 11884 - 11951 1191 . 0004 - 00019 
.14 | . 13681 - 13760 - 13819 . 1375 | . 0007 . 00021 
16 . 15737 . 15765 . 15703 . 1573 . 0003 . 00024 
.18 . 17806 - 17792 | . 17888 . 1783 . 0006 . 00026 
2 . 20071 . 20089 | . 20041 . 2007 . 0003 . 00029 
25 - 24702 | . 24709 | . 24724 | . 2471 . 0001 - 00035 
5 . 49799 . 49762 . 49795 4979 . 0003 . 00065 
1 1.0037 1. 0037 1. 0037 1.004 . 0000 . 0013 
2 1. 9954 1.9972 1. 9948 | 1, 996 | . 0010 . 0029 
3 2. 9968 2. 9975 2. 9993 2. 998 0015 0047 
5 4.9841 4.9923 | 4.9927 4.990 006 0048 
10 | 9. 9624 9. 9671 9. 9647 } 9. 965 . 003 . 0053 
15 14. 991 | 14. 988 } 15. 001 | 14. 99 | . 000 . 0063 
20 19. 963 19. 956 19. 945 19. 95 | 01 . 0093 
25 24. 999 25.031 24. 987 25.01 -02 014 
30 30. 080 30. 049 | 30. 074 30. 07 | .02 O15 
40 40. 354 40. 281 | 40. 367 | 40. 33 | .05 .02 
50 52. 338 52.041 | 52. 336 52. 24 . 20 06 


7. Appendix 


The analysis (see footnote 4) of mismatch errors in the calibration of variable attenuators 
yielded an equation for the error in terms of the scattering coefficients of two fourpoles corre- 
sponding to two settings of the attenuator dial, and the reflection coefficients of the system in 
which the attenuator was inserted. The measurement of all these quantities may be tedious 
or difficult, and an approximate method has been developed. One obtains reasonably close 
limits within which the error lies from a fairly simple experimental procedure. 

The complete expression for the mismatch error is 


-—20 log, | LLShPo)(1= Sil) —(Su)*P Pe |, tia) 
| — Sule) —S22Pz)—(So1)’P Pz | 
where scattering coefficients are denoted by S,,,, and Tg, Ty, represent, respectively, the re- 
flection coefficients of the system “looking towards” the generator and load. Primes are 
used to designate a setting of the attenuator other than the zero or reference setting. 
For small reflections, the following expression is derived from (1a): 





e=20 logo 


1—(S),,—Si) Pe—(Sx—S2)T 14 [SiS S1,S22—(S31)?4 ST I zl. (2a) 


If the attenuator VSWR issmall, the products S’,,;S’.. and S,,; S2. may be neglected. Then 
2a) becomes 


e=20 logo! l a (S,,.—Si ) | aie (Soo— S50) ry T [So,’- ° (S31)? 








| wed ° (3a) 


It is convenient to determine the magnitudes of the individual terms but not their phases, 
so that the limit of error, allowing random phase variations in Sj; and S:2, but none in S. is 


€ =~ 20) logy [1 oe Si— Si i “7 Soo— Soo lr, T ( Soy 2_ So “) [Iz ]. (4a) 
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The quantity |Sj—S;,| is determined as follows. With the attenuator under test con- 
nected as shown in figure 5, and set to its zero or reference position, tuner A is adjusted for a 
detector null, and tuner B is then adjusted until the reflection coefficient T,; of the equivalent 
generator at terminal plane 2 vanishes. (This condition may be recognized by means of an 
auxiliary reflectometer.) 






DETECTOR 








TUNER TUNER 
1 8 Lin NON-REFLECTING 
RA +l ete 
| GENERATOR ATT TERMINATION 
Ry ' T, 











+ - Y ° ° , 
Figure 5. Schematic of system for measuring |Sy,—S,|. 


Movement of the attenuator dial to some other setting will then give an observable 
detector output. Using the theory® of a directional coupler with auxiliary tuners, one can 
obtain |S;:—S},| in the following way. 

The amplitude of the output power from the sidearm of the directional coupler in figure 5 


1+ AT, 


has the magnitude sie] T, 


» where k is a constant for a given stable generator oper- 





ating level. 

With the attenuator set first on zero, T,=Su. Adjusting tuner A for |b;|=0 changes K 
so'that KS,——1. Adjusting tuner B for I,,;=0 makes the dependence of |b;| on I’, take the 
simple form: 





| Tr 
1b,|=k|1——4]- 
— | Su 
Suppose the attenuator dial is moved to a new position such that I, now equals Sj,; then 
S; k aS : 
\bs|=k|1— 7 I= =Tsul Si—Sil|+ The factor k/|Si:| is obtained by replacing the attenuator by a 
O11 |; 






waveguide section containing a sliding short. Upon sliding the short, |b;| goes through 
small variations so that one may observe |b3\max and {bs|min- It is easily shown that 


k 1 
[Sy] 7g (lela + [balmta) 


One can assume that |b3|max=/6s|min and employ a fixed short circuit with negligible error if 
the VSWR corresponding to Sj; is less than 1.15. It may be that this error is tolerable for 
higher VSWR’s since it is not important to know |Si—S},| to great accuracy. 

The quantity |S..—S;,| is found in the same way as above with the attenuator turned 
end for end. 
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6 R. W. Beatty and D. M. Kerns, Recently developed microwave impedance standards and methods of measurement, IRE Trans. on Instr., 
I-7, Nos. 3 & 4, 319-321 (Dec. 1958). 
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Effect of Oleophobic Films on Metal Fatigue 
H. E. Frankel, J. A. Bennett, and W. L. Holshouser 


(January 15, 1960) 


The fatigue strengths of a low-alloy steel, a magnesium alloy, and a copper-beryllium 


alloy were increased by coating the specimens with certain polar organic compounds. 


Also 


the dispersion of the results was much less for coated specimens of these materials than for 


clean ones. 


Similar tests showed no effect with titanium or 6061 aluminum alloy, and only 


a slight improvement for 17-7 PH stainless steel. 
The full beneficial effect of the coatings was found only with compounds having a car- 
bon chain of at least twelve, and this effect was not significantly reduced when the coated 


specimens were tested in water. 


Organic solvents (benzene and xvlene) had a deleterious 


effect on the fatigue life of materials that were improved by oleophobic coatings. 
It is suggested that the effect of the coatings is principally due to their ability to present 


a barrier to water and oxygen molecules. 


1. Introduction 


It is well known that environment has an im- 
portant effect on the fatigue strength of metals. 
Lehmann [1],! in 1926, showed that the fatigue 
strength of a metal was greatly influenced by 
various media surrounding the specimen. Gough 
and Sopwith [2, 3] presented evidence which indi- 
cated that even the mild corrosion involved when 
fatigue tests were run in air affected the results. 

Additional work by Ferguson and Bouton [4] 
verified the findings of the earlier investigations. 
In more recent findings, Frankel and Bennett [5] 
showed that an oii film increased the fatigue life of 
quenched and tempered low-alloy steel. 

Not all environments which exclude air are bene- 
ficial, however. Karpenko [6] and Karpenko and 
Karlashov [7] showed that certain surface active 
media such as light oil+2 percent of oleic acid or 
water+2 percent of isopentanol can cause a reduc- 
tion in the fatigue strength of structural steel. A 
recent review by Gilbert [8] summarizes much of 
the published work on corrosion fatigue. 

A study of the literature revealed that, in many 
studies of the effect of environment, the media used 
were complex mixtures, not suitable for mechanism 
studies, and the data were valid only for a specific 
set of conditions involving many variables. In an 
effort to obtain more basic information on the effect 
of surface reactions on the fatigue process, the 
National Aeronautics and Space Administration 
sponsored an investigation at the National Bureau 
of Standards. This paper contains a report on one 
phase of this investigation. 

Zisman and his associates [9] at the Naval Research 
Laboratory have made extensive investigations of 
the phenomena occurring at liquid-solid interfaces. 
As a result of this work, they found that certain 
classes of polar organic compounds form films on 
metal surfaces that are unwetted by, or are oleo- 
phobic to, all but the lowest-boiling hydrocarbons 


a, 


! Figures in brackets indicate the literature references at the end of this paper. 





and other nonpolar liquids. These films are also 
highly hydrophobic, and are known to be effective 
in preventing corrosion. It was thought important 
to determine the effect that they might have on 
fatigue. 


2. Materials and Testing Procedure 


Six alloys were selected for the investigation, as 
follows: 


Material Condition 
SAE 4340 steel Quenched, double 
tempered 


Aluminum alloy 6061- Solution treated and 
T6 aged 

Titanium (75—A) Hot rolled 

17-7 PH stainless Annealed (condition 

A) 

Copper —1.75% Solution annealed 
beryllium 

Magnesium (AZ63-— 
H24) 


Hot rolled 


All of the tests were conducted on rotating-beam 
fatigue testing machines of the R. R. Moore type, 
which apply pure bending moments to the test 
section. The machines were operated at approxi- 
mately 2,000 rpm to minimize overheating effects. 
The specimens were machined with a smooth test 
section in which the longitudinal radius of curvature 
was held constant at 314 in.; the minimum diameter 
was 0.310 in. for magnesium and 0.250 in. for the 
other materials. The straight cylindrical shanks 
were held in collet-type chucks in the bearing boxes 
of the machines. After the specimens were rough 
machined, they were ground to size with a 7-in. 
diameter wheel, the plane of which was parallel to the 
length of the specimen. To avoid overheating the 
surface, very shallow cuts were used and the entire 
operation carefully controlled. The preparation of 
the specimens was completed by polishing in a 
longitudinal direction with 400 Aloxite paper in the 
automatic polishing device previously described [10]. 
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For the tests with liquids in contact with the 
reduced section of the specimen, a short length of 
tight fitting polyethylene tubing was slipped over the 
cylindrical shanks. Prior to testing, the specimens 
were individually cleaned by agitation in_ three 
successive containers of high-purity benzene. The 
pieces of polyethylene tubing were cleaned in a 
similar manner; the specimens and tubing were 
kept in benzene until required for testing. 

A test was run with the specimen under one of 
the following conditions: (1) Cleaned and dry, (2) 
coated with or immersed in an organic substance, or 
(3) immersed in water in either the clean condition 
or after coating with an oleophobic film. Whenever 
a liquid medium was used, an air bubble large enough 
to contact the specimen was deliberately left in the 
tube. The organic substances listed below were 
included in this study: 


Name Formula Ty pe 
n-Hexadecane__ CH;(CH,)14 CH;-.. Nonpolar 
Xylene___. a C,H,4(CHs)>_ - - Do. 
SS: °C ae Do. 
Octyl alcohol CH;(CH,),; CH,OH_ Oleophobic 
Dodecyl aleohol CH ;(CH,),;; OH-_--- Do. 
Dodecylamine CH,;(CH,),, NHo2- - - Do. 
Octadecylamine CH;(CH.,),, NH3- — Do. 


As this investigation was primarily concerned with 
surface reactions caused by fluctuating stress, no 
studies were made of environments more corrosive 
than air and water, nor of the effects of time or 
speed of testing. 

All of the results obtained in tests at constant 
stress amplitude are presented in table 1. The 
water and nonpolar liquids either had no significant 
effect or caused a decrease in fatigue life. The 





Material 
Stress Clean = 
amplitude 
1 2 
1,000 psi 

4340 steel___- 95 181 ‘a 
100 83 111 71 
115 28 

Mg AZ-63 -- 10 500 246 
12. ) 5l 
15 17 

Cu-1.75Be 45 258 
5O 89 108 58 
60 21 

Ti-75A 55 989 
75 13 

Al 6061-T6 30 326 310 
37 49 

17-7 PH 100 114 
115 52 


n-hexadecane 
xylene 

. benzene 

. tapwater 

. octyl] alcohol 


a1, 
9 


1 


oleophobic films, however, caused a marked increase 
in fatigue life for three of the six materials investi- 
gated; namely, the steel, the magnesium alloy, and 
the copper-beryllium alloy. There was no sig- 
nificant difference in the results between specimens 
immersed in these substances or merely coated with 
them, and the results are not differentiated. The 
aluminum and titanium showed no effect of any of 
the environments, while the 17-7 PH _ stainless 
showed only slight effect. 

Coating of the specimens with dodecyl alcohol and 
then testing them in a tapwater environment resulted 
in as much improvement as was found with specimens 
tested in air. It was also found that tapwater had 
very little, if any, influence on the fatigue properties 
of those materials not affected by the oleophobic 
films. Continuous application of xylene or benzene 
to specimens of the first three metals listed in table 
1 resulted in the poorest fatigue behavior. 

Figure 1 shows a comparison of the effect of re- 
agents having different carbon chain lengths. Ap- 
parently a length of about 12 carbon atoms provides 
optimum conditions with little change in fatigue life 
beyond this number. It can also be seen that the 
type of polar group does not appear to be as im- 
portant as the number of carbon atoms. 

The data in table 1 and figure 1 indicated that the 
results for specimens coated with dodecyl alcohol, 
dodecylamine, and octadecylamine were sufficiently 
similar to permit grouping for better statistical 
analysis. This has been done in table 2 for four 
materials. There were not enough specimens of the 
other materials to justify such an analysis. In the 
materials that were affected by the oleophobic films, 
the dispersion of the data, as measured by the co- 
efficient of variation, was much smaller for the coated 
specimens than for clean ones at the same stress 
amplitude. 








Median life, 108 cycles 





Coatings 4 
3 4 5 6 7 8 9 
SO 417 432 
69 76 141 189 176 183 173 
25 40) 37 
430 1, 776 1, 680 1, 543 1, 300 
20 31 187 410 315 293 367 
12 ae 
1,195 1, 450 1, 320 
62 92 151 248 226) 221 226 
21 
806 
14 
317 350 307 
59 63 58 
114 169 
55 . 67 
6. dodecyl alcohol 


7. dodecylamine 
. octadecylamine 
. dodecyl alechol plus water 
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TABLE 2. Evaluated results 





Fatigue test 
x | Coefficient of variation 





















































| 
Material | Stress ampli- Clean | Coated > 
| tude | a 
| | No. of Life,» No. of | No. of Life,* No. of Percent Clean Coated > 
| | specimens | cycles | specimens | cycles | improve- | 
| | ment 
a | os — oon | 
| 1,000 psi | | 
ee = oo. | 11 | 177X103 6 | 433X103 | 144 0. 098 0.015 
100 | 87 6 | 180 107_—‘| . 040 014 
| 115 | 12} 26.6 | 5 | 39 46.8 | . 079 . 030 
Mg AZ-63.......| 10 | 7| 477 7| 1,524 | 242 | 059 | 022 
12.5 | 5 | 48 | 7 | 337 603 | 110 - 038 
Cu-1.75Be.--.--- | 45 | 5| 272 | 5| 1,268 | 265 | 063 033 
| 50 | 5 | 89 7 | 231 159 . 089 . 019 
Al 6061-T6......| 30 | 6 311 3 333 7 . 032 048 
| 37 | 10 | 50 6 54 8 076 072 
| | | | 
« Antilog of the log mean. 
b’ Coated with oleophobic substances having carbon chain lengths of 12 or greater. 
The data for 4340 steel are presented graphically 10 T 
in figure 2, and it will be noted that the vertical Uy as | ear ert 
difference between the S—N curves is nearly constant, | 
as the improvement in fatigue life increases at low i te ee | 1s cpio 
. ~ za } 
stresses. The approximate ranges of +1 standard wis a, | 
° . » . . . w 
deviation from the means are shown in this figure, 3 3 saiell | 
emphasizing the small dispersion of results with “a> T hes, | 5 Geman aba 
= w 
coated specimens. ra 7s 
. . = al : 
An attempt was made to determine the number of — ~ | | 
1: : te. : ° ~ w e 0 o | 
cycles at which cracking started in at least one test 2\3 oe’ 9 > 
“ . ° ° S re} 
of each of the materials in the clean and coated cle 2 4 ee ae Eee! ST 
conditions. This was done by measuring the vibra- a | eae 
tion of the bearing boxes: there is a marked increase * 
n vibrati then cracki te. These measure 7-7 PH 
in vibration when cracking starts. 1ese measure- 
ments showed that approximately 85 to 90 percent Ds & | | 
of the total life of the specimen had been expended 0 4 8 2 6 20 


prior to the initiation of cracks. Thus, it can be 
concluded that, with the specimens used in this in- | _ : “So 
vestigation, the improvement obtained with oleo- | F!6uRe 1. Effect of oleophobic films on median fatigue] life 
: ; . . . . for four materials. 
phobic coatings was primarily a measure of their = 
effect on the initiation of a fatigue crack. g Ate gises apples in the teste were ot follows: Me (AZ-18)—1200 pal 
Zisman and coworkers [11] characterize an oleo- | seel_20n0 ps. (For each, material, ehe upper value a a chain length of 1 
phobic substance as a long chain, aliphatic molecule ; 
having a polar group at one end. A consequence of 
this structure is a strong tendency to form tightly 
packed monomolecular films on various substrates. 
It has been pointed out that these materials offer 
the same degree of protection to specimens immersed 
in water as to those in air. Those specimens, given 
no protection, had generally shorter fatigue lives 
when tested with their surfaces partially in water. 
Thus, it is reasonable to assume that at least some 
of the effect of the oleophobic materials is due to 
their ability to present a barrier to water and oxygen 
molecules. 


CARBON CHAIN LENGTH 
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| } SAE 4340 | 
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STRESS AMPLITUDE, 1000 psi 
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Likewise, the reduction in dispersion of results | 
from coated specimens may be due to the replace- z err : 
, > varis » . "Ty ¢ . ra r] 2 
me nt of the variable laboratory atmosphere with a nunaeen tr evea ne enadeuee 
reproducible oleophobic film of the substance. 
Similar results were obtained by Liu and Corten [12] | Fiaure 2. Number of cycles to fracture as a function of the 
who reported a marked reduction in scatter when they stress amplitude for SAE 4340 steel. 
coated their specimens with vaseline. 
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Standard deviations at + 1¢ are also shown. 
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In contrast to the results with steel, Cu-Be alloy, 
and magnesium, the fatigue strengths of the 6061 
aluminum and the titanium were not affected by 
coatings of the oleophobic substances used. The 
improvement with the 17-7 PH stainless, if any, 
was small. This suggests that the effects observed 
may be associated with the nature of the oxide film 
on the metals; the oxides on the unaffected alloys 
are known to be more adherent and less permeable 
than those of iron, copper, and magnesium. 

On the basis of the limited data in table 1, it 
appears that the organic solvents xylene and benzene 
have a more deleterious effect on the fatigue strength 
than does water. This indicates that there are 
phenomena other than the reactions with water and 
oxygen that are important in the effect of environ- 
ment on fatigue, but the exact nature of these phe- 
nomena cannot be determined without additional 
experiments. 


3. Conclusions 


This study of the effect of certain types of polar 
organic substances on the fatigue life characteristics 
of six different alloys revealed the following: 

1. A general group of organic materials which 
form oleophobic films on me ‘tal surfaces were found 
to be effective in increasing the fatigue strengths of 
SAE 4340 steel, AZ-63 magnesium, and a copper 
1.75 percent beryllium alloy. No effect was ob- 
served with titanium (75—A) and 6061-T6 aluminum 
alloy, and the effect on 17-7 PH stainless steel was 
questionable. 

Conclusions 2 to 5 apply to results on the three 
alloys affected by the coatings. 

2. Polar compounds with carbon chain lengths 
of 12 and 18 produced a greater improvement in 
fatigue life than a compound with a chain length of 


8. As two different polar groups, amines and 
alcohols, showed the same improvement, it seems 


likely that the length of the carbon chain may be 
one of the controlling factors. 

The dispersion of the data from the coated 
specimens was much smaller than that from the 
cl 6 ones tested at the same stress amplitude. 

The mean fatigue life of specimens with oleo- 
sila films was the same whether they were tested 
in water or in air. 

5. With the specimens used in this investigation, 


the fatigue life was governed primarily by the 
number of cycles to initiate a crack. Thus, the 


improvement obtained with oleophobic coatings was 
largely due to increase in the number of cycles to 
the initiation of a fatigue crack. 





6. It is suggested that oleophobiec films inhibit 
the impingement of extraneous molecules, such as 
water or oxygen, upon the surface. This reduces the 
rate of detrimental surface reactions that normally 
occur when specimens are stressed in air. : 





The authors thank Dr. J. Kruger, National 
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Ratio-Recording Spectroradiometer 
Harry K. Hammond III, Warren L. Holford, and Milton L. Kuder 


(January 26, 1960) 


A ratio-recording spectroradiometer has been constructed primarily for determining the 
relative spectral irradiance from fluorescent lamps in the visible spectrum. Radiant flux 
from a test source and a comparison source is transmitted or reflected by separate diffusers. 
The irradiance from each diffuser is sampled alternately of the order of 100 times a second 
by a double prism monochromator with cam-linearized wavelength drive from 360 to 760 mu. 
A 14-stage multiplier phototube with S—-20 response is used with an electronic gate and 
integrator circuits to compare the spectral irradiances from the two sources. When the 
phototube outputs are unequal, a servo unit adjusts apertures in each beam to equalize them. 
The measured parameter is the amount of adjustment required for equalization at each 
wavelength. The apertures are adjusted by a cam coupled to a pen which records the ratio 
on a three-cycle logarithmic strip chart. The instrument requires about 8 minutes to 
record thespectrum at 10myperinch. An expanded scale of 1 my per inch is used for evaluating 
the energy in the spectral lines. The speed of wavelength scan may be made inversely 
proportional to the unbalance signal, if desired, to provide ample time to record spectral 


lines accurately. 


Symmetry of beam treatment is demonstrated by interchanging test and 


comparison sources and recording the inverse ratio. 


1. Introduction 


In order to calibrate fluorescent lamp color stand- 
ards on a sound physical basis, the relative spectral 
irradiance from the lamp should be determined. 
These data, when weighted by the response functions 
of the CIE standard observer, yield chromaticity 
coordinates directly. A spectroradiometer has been 
constructed at NBS for the measurement of relative 
spectral irradiance with particular regard to the 
problem of fluorescent-lamp measurement. 

Spectrometers for dispersing light into its spectral 
components have been commercially available for a 
long time. However, a universal design has not 
been adopted for converting a spectrometer into a 
spectroradiometer because of the varied requirements 
of particular laboratories. A recording spectro- 
radiometer was built by Zworykin in 1939 [1]... Ten 
years later, Studer and Jacobsen reported the design 
of a spectroradiometer which was available com- 
mercially for a time [2]. In 1957, however, no 
spectroradiometer was commercially available; so it 
became necessary to design one to meet our require- 
ments. 


2. Requirements 


The general requirements of a spectroradiometer 
are simply stated. Radiant flux from a source must 
be dispersed into its spectral components and 
evaluated with a suitable detector and read-out 
device to provide relative spectral irradiance data. 
Accurate measurements require a high degree of 


spectral purity, but sufficient energy must be trans- 


—_____. 


Figures in brackets indicate the literature references at the end of this paper. 
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mitted by the dispersing element for precise evalua- 
tion by the detector. A practical compromise be- 
tween resolution and sensitivity is necessary, as very 
high resolution would permit only infinitesimal 
energy to be transmitted to the detector. The in- 
strument response must be appropriate to record 
accurately the peaks of sharp spectral lines. It 
must be able also to average the light from a modu- 
lated source, such as a fluorescent lamp operated on 
alternating current. 

For the purpose of evaluating fluorescent lamp 
chromaticity, a photometric accuracy of 1 percent is 
considered adequate. If the photometric errors are 
truly random, then the chromaticity coordinates 
computed from the spectral radiance data will be 
uncertain by much less than this amount. 


3. Basic Components 


Monochromator. Both prism and grating mono- 
chromators are available for dispersing light into a 
spectrum. Grating monochromators are attractive 
because the dispersion is almost constant over a 
practical range of wavelengths. This permits a 
direct connection to a linear wavelength scale and 
provides a pass band of constant spectral width for 
fixed mechanical slits. Gratings, however, have the 
disadvantages of overlapping orders, higher stray 
light, and possible ghosts. Stray radiant flux in any 
dispersion system is minimized by double dispersion. 
In a good double prism monochromator, the stray 
flux will be several orders of magnitude below the 
radiant flux in the pass band. After some investiga- 
tion of monochromators available at a moderate 
price, a Carl Leiss double monochromator with dense 
flint glass prisms was selected. Dispersion data are 








shown in figure 1. The three bilateral slits have 
straight sides 10-mm long and are set independently. 
Data are taken with constant mechanical slits; 
hence the pass band is a function of wavelength. 
For measurement of fluorescent lamps, the entrance 
and exit slits are frequently set at 0.5 mm with the 
middle slit at 0.6 mm. For these conditions, the 
pass band at 546 my is just under 3 my when evalu- 
ated as the width of the mercury line at half height. 
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Figure 1. Dispersion data for spectroradiometer. 


Carl Leiss double monochromator with dense flint glass prisms. 


The original nonlinear prism drive, consisting of 
sector and worm, has been disengaged. The new 
drive consists of a system of levers including an 18- 
in. arm that follows a specially constructed cam 
with radii ranging from 3 to 7 in. The arm, cam, 
and associated gear train render the wavelength 
scale linear over the range 360 to 760 mu. 

Detector. The low level of radiant flux passing 
through the exit slit of the monochromator requires 
that a highly sensitive detector be employed. The 
detector should also provide good response over the 
entire range of spectrum to be measured. The 
equivalent noise input should be low to obtain a 
high signal-to-noise ratio for weak signals. The 
detector selected was a 14-stage multiplier photo- 
tube, RCA developmental type C7237, now desig- 
nated as type 7265. This tube is an end-on type 
with S—20 response, peaking in sensitivity at approx- 
imately 420 my and maintaining 20 percent of the 
peak response at 700 mu. 

Recording System. Some spectroradiometers have 
been designed to use the output of the detector as 
the measured parameter. Instruments can be de- 
signed, however, that are essentially independent of 
the response characteristics of the phototube and 
the spectral transmittance of the monochromator. 
Considerations involved in this type of system have 
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been published by Mitsuhashi in 1955 [3], and an 
instrument utilizing these principles was reported 
by MacAdam in 1953 [4] and described in detail in 
1958 [5]. Our design was developed independently, 
but is similar in many respects to those cited above. 
Measurement is accomplished by rapidly comparing 
the flux from a standard source with that from a 
test source as the spectrum is scanned and optically 
attenuating the beams until the two have equal 
flux. For this system, the measured parameter igs 
the setting of the beam attenuator. Attenuator 
settings are recorded on a strip chart, giving a plot 
of ratios of relative spectral irradiance from the two 
sources as a function of wavelength. 

Attenuator System. There are two requirements 
that should be considered in the design of an atten- 
uator system. The geometric distribution of flux 
incident on the entrance slit of the monochromator 
should not be appreciably affected by the attenu- 
ator, and the calibration of the system should not be 
a function of wavelength. The latter requirement 
cannot be met by optical wedges or polarizing filters, 
The former requirement cannot be met by a single 
adjustable aperture but for practical purposes is 
satisfied by an array of apertures. To compare pre- 
cisely sources whose spectral irradiance ratios vary 
over a wide range, the aperture arrangement should 
be chosen so that the transmittance will be a rapidly 
changing function of displacement. Our basic at- 
tenuator design consists of an array of 25 square 
openings, io in. on a side, etched in metal foil, and 
arranged as shown in figure 2. Each attenuator 
consists of a pair of arrays, one fixed and one movable 
on an are whose radius is about 4 in. The change in 
area of the openings, and therefore the transmittance, 
is proportional to the second power of the displace- 
ment of the movable array. When arranged so that 
the attenuator in one beam opens as the other closes, 
the transmittance ratio is a fourth-power function of 
the displacement drive, provided the flux through 
the attenuators is relatively uniform. 





FIGURE 2. Allenuator array. 
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4. Entrant Beam Geometry 


The entrant beam geometry is shown in figure 3. 
The rotating sector beam switch, consisting of a 
wheel with four sectors of alternate apertures and 
mirrors, alternately passes light from each source 
into the monochromator. The irradiance from each 
lamp may be sampled at a rate which is adjustable 
from 20 to 120 times a second. Compartments for 
filters are located in each beam for use when desired. 
The optical attenuators are located in each beam and 
operate in opposition, one opening while the other is 
closing. The two light sources are external to the 
instrument case and, therefore, can be easily changed. 








The cone of rays at the entrance slit of the mono- 
chromator required to fill the collimating mirror was 
determined by rear projection of a source placed at 
the exit slit. To prevent the monochromator from 
giving a different weighting to rays from different 
sources, it is necessary that they both subtend the 
same solid angle at the entrance slit, or better, that 
they fill the maximum solid angle of the mono- 
chromator opties uniformly. To achieve this goal at 
some distanee from the slit requires that the sources 
have either relatively large areas of uniform radiance 
or that they uniformly irradiate reflecting or trans- 
mitting diffusers of appropriate size. 
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FIGURE 3. 


5. Photometric Scale 


In order to achieve the desired photometric scale 
on the strip chart, a precision linkage has been in- 
stalled between the attenuators and the servo drive 
mechanism. The two movable elements of the at- 
tenuators are displaced simultaneously by two iden- 
tical eccentric disks on a common shaft. The shaft 
extends to the strip chart recorder where it terminates 
in a cam-following lever. A relatively simple cam, 
driven by the recorder servomotor, converts the 
fourth power function to the desired logarithmic 
function. The location of cam and follower are 
shown in figure 4. The recorder plots the area ratio 
of the attenuators on a logarithmic scale. Cams 


may be interchanged by the removal of one thumb | 


nut to provide different photometric scales. For 
fluorescent-lamp measurements, however, the cam 
consistently used provides a three-cycle log scale on 
the 10-in. width of the strip chart. The intense 
lines of the mercury: spectrum are easily accom- 
modated with this scale. 
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Figure 4. Photometric cam provides 3-cycle log scale. 
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6. Electronics 


The dynode voltage supplied to the phototube is 
continuously adjusted automatically to maintain 
constant average anode current. This feature 
provides constant detector sensitivity as spectral 
irradiance from a source decreases, until a minimum 
photon level is reached and maximum gain is applied. 
This feature also prevents the phototube from giving 
erroneous balance indications due to anode current 
saturation caused by high gain at high flux levels. 
Phototube damage due to excessive gain at such 
levels is also prevented by this automatic voltage 
adjustment. 

Signals from two sources are separately supplied 
to two integrating networks by electronic gates as 
shown in the block diagram of figure 5. The gating 
signals are generated by sets of auxiliary light sources 
and phototubes. One set is used for gating the test 
source signals, another for the signals from the 
reference source. An opaque disk on the beam 
switch motor shaft has precisely positioned apertures 
which produce the gating signals. The gate 
opened after the test or reference flux falls on the 
detector and is closed before the flux is switched off. 
The gated signals are essentially square waves. The 
integrator circuit and subsequent filter are designed 
to provide a smooth direct-current signal to the 
conventional servoamplifier of the recorder. The in- 
tegrating network has a time constant of about 
16 sec. The integrator also serves to average signal 
variations due to modulated radiant flux from an 
a—c operated fluorescent lamp as well as noise pulses 
from the phototube. When the integrated signals 
from the two sources are unequal, the servo drives 
the attenuators until a balance is achieved. 

The magnitude of the unbalance signal is also used 
to control the rate of wavelength scan. One motor 
drives the prisms and the recorder chart and thus 
wavelength correlation is maintained between the 
spectrum and the chart regardless of the scanning 
rate. The speed of the drive motor is modulated 
according to the magnitude of the unbalance signal. 
This modulation is roughly proportional to the 
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reciprocal of the first derivative of the curve on the 
chart. When the curve is flat the rate of scan is 
maximum. As the slope of the curve begins to 
change the rate of scan is automatically reduced to 
prevent inaccurate recording imposed by phase lag 
in the servo response. This feature insures the 
faithful recording of the peaks of sharp spectral 
lines. 

Tachometer generators are used in conjunction 
with both the wavelength drive and photometric 
balance motors. One generator provides motor 
speed information to the wavelength-scan speed- 
control unit. The other generator provides velocity 
damping for the balance motor. ; 

The wavelength drive has been equipped with 
auxiliary gears and electric clutches to provide two 
chart 1 my/in. or 10 may/in., merely by 
energizing the proper clutch. The expanded scale 
of 1 my/in. is most useful for evaluating the flux in 
spectral lines. When a modulated scanning rate is 
not desired, a manual switch provides for selection 
of one of two uniform scanning rates. With a chart 
scale of 10 my/in. the approximate times required 
to sean the visible spectrum at each speed are 8 min 
and 30 min. 


scales, 


7. Meters and Controls 


The front panel of the console is shown in figure 6. | 
The wavelength motor shaft has been extended to the | 
front of the instrument to permit convenient manual — | 
setting toany desired wavelength. Two dials are used, 
one reading hundreds and tens of millimicrons, the 
other reading units and tenths, so that wavelengths 
can be read easily to better than 0.1 my. Calibra- 
tions with spectral-line sources, however, show that 
the original cam has wavelength-scale errors of as 
much as 0.2 mu below 580 my and somewhat larger 
errors at longer wavelengths. 
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FIGURE 
instrument records spectral irradiance 
relative to an incandescent lamp. 


while 
lamp 


Holford records electrical parameters 
of fluorescent 


Photograph was made prior to use of external reflecting diffuser. 
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The channel balance meter in the lower left indi- 
cates the difference voltage between the test and 
reference signals. During operation, the difference 
voltage indications will be very small if the servo 
balance unit is operating properly. 

The meter on the right indicates the dynode supply 
voltage to the phototube. Since the average anode 
current is held constant automatically, a lower flux 
level requires a higher dynode voltage. As the 
dynode voltage is raised to the maximum permissible 
value, the signal-to-noise ratio is lowered slightly. 
Further reduction in flux level results in a noisy 
signal and consequently produces an erratic recorder 
trace. Thus the dynode-voltage meter provides an 
indication of whether there is enough flux to permit 
a reliable ratio determination. 

The meter at the upper left also provides indica- 
tions proportional to the dynode voltage. This 
meter, however, is provided witb zero suppression 
and sensitivity controls so as to provide an indication 
of the location of a spectral line when manually 
manipulating the wavelength knob during wave- 
length calibration. 

Switches for controlling power, speed of scan, and 
wavelength scale are conveniently located on the 
front panel. Pushbuttons to provide fast travel 
of the wavelength drive in either the forward or 
backward direction are located near the wavelength 
dials. Potentiometers for setting electrical charac- 
teristics in the balance circuit and in the speed con- 
trol unit to provide optimum instrument perform- 
ance are located behind the panel below the two 
meters. 


8. Calibration 


During construction of the instrument, particular 
care was given to the design of the components whose 
calibration directly affected the instrumental ac- 
curacy. Mechanical linkages in the wavelength 
drive and photometric system were held to close 
tolerances and materials chosen were those least 
likely to wear or change with time. Several hundred 
hours of instrument operation have revealed no 
changes due to wear. Frequent wavelength calibra- 
tions using narrow slits of 0.15 mm reveal that wave- 
length scale calibration is reproducible to better 
than 0.1 muy. 

The photometric cam is calibrated by utilizing 
ratios of 2 to 1. Two lamps of equal intensity, as 
measured by the spectroradiometer, irradiate the 
diffuser at one entrance port. A reference lamp il- 
luminating a diffuser at the second port is compared 
with each of these lamps in turn to assure their 
equality, and then compared with the sum of the 
two lamps. The resulting irradiance ratio is pre- 
cisely 2. The wavelength drive is not operated 
during the photometric calibration. The three-cycle 
log chart is adequately covered by nine steps of 
ratios of 2to 1. The intermediate calibration points 
are interpolated with data based on micrometric 
measurements of the cam radius. The reproducibility 
of photometric calibrations over the usable portion 
of the cam is about 0.3 percent. 





9. Performance Tests 


Tests have been made to evaluate the performance 
of the spectroradiometer. The investigation of pos- 
sible systematic errors has been aided by the fact 
that the test lamp may be operated at either of the 
two entrance ports with the reference lamp at the 
other. Because photometric accuracy is determined 
by the geometric distribution of the radiant flux at 
the attenuators and the calibration of the photo- 
metric cam, the exchange of entrance ports produces 
reciprocal conditions which provide internal checks 
on instrument performance. A small discrepancy in 
fluorescent-lamp measurements was revealed when 
the reference and test sources were interchanged. 
This caused an investigation to be made of possible 
systematic errors or bias in the instrument. The 
electronic circuits for each channel were interchanged, 
but these circuits were not the cause of the biased 
results. The electronic system was further vindi- 
cated when practically identical spectral-irradiance 
data were recorded from a fluorescent lamp operating 
first on alternating current and then or direct current. 
A part of the small discrepancy was attributed to a 
bias in the servo unit. Equal, but opposite voltages 
to the servoamplifier did not produce equal torques 
at the servomotor because the voltages on the motor 
windings were not exactly 90° out of phase. The 
effect was noticeable only when the signals were 
relatively weak. A trimmer capacitor placed across 
one of the motor windings produced the required 
phase relationship for balanced operation of the servo 
motor. 

Further investigation revealed that interchange of 
reference and test sources introduced no discrepancy 
when the test source geometry is similar to that of 
the calibration source. A fluorescent lamp, however, 
has a geometry quite different from that of the 
incandescent lamp standard. The discrepancy be- 
tween results obtained by interchanging test and 
reference sources was traced to the nonuniform 
irradiation of the diffusing block placed at an angle 
of 45° to the axis of the fluorescent lamp. The non- 
uniform irradiance distribution from the incandescent 
standard was included in the calibration, but the 
entrance port acted as a field stop for the fluorescent 
lamp, thereby producing a different irradiance 
distribution. 

In an attempt to eliminate this difficulty, the 
reflecting diffusers were replaced with mirrors, and 
ground-glass transmitting diffusers were mounted 
parallel to the axis of the fluorescent lamp. Ground 
glass, however, acts like a large number of tiny 
prisms and thus deviates shortwave energy more 
than longwave energy. This differential diffusion 
can be tolerated when comparing sources of the same 
size. When test and standard sources are different in 
size, as is the case with fluorescent and incandescent 
lamps, then different amounts of radiant flux are 
scattered out of the beam at each wavelength and 
erroneous values of relative spectral irradiance are 
obtained. Further consideration of the use of reflect- 
ing diffusers indicated that the problem of non- 
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uniform irradiance of the test diffuser could be solved | 
by relocating it external to the instrument in such a 
position as to obtain uniform irradiance from the 
fluorescent lamp. 


10. Modes of Operation and Applications 


Accurate determination of the ratios of relative 
spectral irradiance from a test and a_ standard 
source by direct comparison requires similar treat- 
ment of each beam. Light from each source must 
traverse equivalent paths to prevent erroneous 
results caused by spectral selectivity of the system. 
No evidence of asymmetry has been discovered, 
but the majority of measurements have been made 
by using a substitution technique. In this tech- 
nique a reference source is calibrated by operating 
the standard source on the test side. This procedure 
has the disadvantage of requiring periodic cali- 


bration of the reference source, but it has four 
advantages: (1) The standard and test beams 


travel the same paths; (2) the life of the standard 
lamp is conserved by operating it only to calibrate 
the reference source; (3) the reference source may 
be a low-wattage, short-life lamp which produces 
less heating of the air in the vicinity of the fluores- 
cent lamp than does the usual 500-w projection 
lamp regularly supplied as a_ color-temperature 
standard; and (4) the spectral irradiance from the 
reference source can be modified by an uncalibrated 
filter. The use of a blue filter reduces the flux from 
the reference source at the red end of the spectrum, 
thereby keeping the flux ratio of fluorescent to refer- 
ence lamp on scale. If a filter were not used, the | 
ratio of flux of incandescent standard at 2854°K 
to reference source would plot as a horizontal line 
throughout the spectrum and the curve for fluores- 
cent lamp ‘“B’’ would be off the bottom of the seale 
beyond 700 mu. See figure 7. 

The spectroradiometer was built primarily for 
the measurement of fluorescent lamps. However, 
the instrument can measure the relative spectral 
irradiance from any source or souree-filter combina- 
tion that will provide a minimum illumination of 
about 1 fe on the diffuser. The transmittance of | 
filters can be obtained by computing the ratio of 
relative spectral irradiance from the source-filter 
combination to that from the source alone. <A 
“100 percent” curve is recorded as the irradiance 
ratio of two sources. If the sources are identical, 
the irradiance ratio will be constant throughout 
the spectrum, and the plot will be a line parallel to 
the wavelength axis. Transmittance data may be 
recorded by inserting the filter in either of the two 
beams. By measuring the filter in each beam and 
averaging the result, any inaccuracies in the photo- 
metric-cam calibration above and below the 100 
percent curve will be averaged. Comparison of 
two filters may be made directly by inserting one | 
in each beam. 

Reflectance data for a specimen illuminated with | 
any desired flux geometry may be obtained by com- | 
puting the ratio of spectral radiance of the specimen | 
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Figure 7. Curves obtained for two daylight fluorescent lamps 


and an incandescent lamp at 2854° K. 

These curves are ratios of the irradiance from the test lamp to that from an 
incandescent refereneée lamp and Corning blue-green filter type 4308. Filter is 
used to limit irradiance-ratio range of lamps to two log-cycles. Daylight fluores- 
cent lamps usually do not differ by as much as shown here. Curve A was ob- 
tained with an old lamp utilizing a discontinued bery!lium phosphor. Curve B 
was obtained with a lamp from current production utilizing a phosphate phos- 
phor. Low values of irradiance from both test and reference sources below 
400 mw and above 700 my cause the trace to be noisy. 


to that of a standard diffuser. Luminescent. or 
fluorescent materials can also be measured if they 
have several inches of uniform luminance and if 
the luminance is of the order of 1 fL. 

The chart drive incorporates several flexible 
features. As previously mentioned, data may be 
recorded with wavelength scales of 10 may/in. or 
1 myu/in. The recorder chart can also be driven 
without advancing the wavelength so as to produce 
a record of the radiant flux of a spectrally homoge- 
neous signal as a function of time. 

Cams have been designed to provide 1 or 3 log 
cycles on the 10-in. strip chart. More precise data, 
however, may be obtained by reading directly two 
dial scales; one, a coarse scale, is on the pen-cable 
drum and the other, a fine scale, is on the servo 
motor shaft. The precision of photometric readings 
is approximately 0.1 percent on the dial seales. 


11. Discussion 


Spectral-irradiance data for a light source can be 
obtained relative to an incandescent lamp of known 
color temperature in relatively short time with good 
precision. Reduction of data from the chart record, 
however, is time consuming and in addition appears 
to be the weakest link in the chain of data precision. 
As indicated above, dial readings of servo motor bal- 
ance positions have been used for calibration and on 
other occasions to provide readings of highest preci- 
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sion. Photometric calibrations can be made with a 
reproducibility of a few tenths of a percent, where- 
as chart readings are uncertain by about half a per- 
cent. Greater precision and increased speed of data 
reduction could be obtained by utilizing an electronic 
encoder to obtain decimal or binary digital readout 
on punched paper tape or other high precision record- 
ing means. Data in this form could then be sup- 
plied without transcription to a high-speed electronic 
computer to obtain chromaticity coordinates. 


The authors express their thanks to Paul E. 
Allison for his contribution of skilled craftsmanship 
in the fabrication of mechanical parts for the instru- 
ment and to Jerome S. Laufer for his assistance 
in the selection of some of the components for the 
instrument. 
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An Intermittent-Action Camera With Absolute 
Time Calibration 
G. Hefley, R. H. Doherty, and E. L. Berger 


(January 12, 1960) 


A detailed description is presented of a film-recording system in which a randomly 


occurring event and its absolute time are recorded simultaneously. 


The system consists of 


a 16-millimeter framing camera capable of intermittent operation at a maximum rate of 
140 frames per second (fps) and a clock capable of reading out time with an absolute accuracy 


of plus or minus 1 millisecond (msec). 


1. Introduction 


The Central Radio Propagation Laboratory of the 
National Bureau of Standards, Boulder Laboratories, 
has developed a recording system which will simul- 
taneously photograph randomly occurring events 
and the precise time of each event (fig. 1). The 
system employs a camera which has a minimum 
consumption of film but can complete a framing 
evele in 7 msec. The timing system is capable of 
resolving time to within a few microseconds. When 
the standard time broadcasts of WWYV are used to 
synchronize the system, the absolute accuracy is 
dependent upon the propagation variables of the 
WWYV signal. 

The high-framing rate of the camera results from 
the use of a unique clutch system between the film 
drive and a continuously revolving flywheel. The 














FIGURE 1. 


Camera and timer. 
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clutch actuation mechanism is similar to a large 
electrodynamic speaker. When a pulse is applied 
to a part corresponding to the voice coil, a cone is 
driven into engagement with the flywheel, thus 
rotating a small drum in contact with the film. The 
rotation of the drum pulls one frame of film past the 
lens opening. The peak power applied to the voice 
coil is approximately 35 kw. Power surges of this 
magnitude may be repeated at intervals of only a 
few milliseconds by utilizing 2 condenser reservoirs 
and 2 thyratrons in a series arrangement. An 
88-microfarad (uf) condenser is discharged through 
the voice coil by one thyratron and after a short 
delay is recharged from a 1,600 yf reservoir by the 
other thyratron. Consequently, the power supply 
required to operate the system needs only to provide 
an average current sufficient to maintain the charge 
on the large condenser reservoir at the maximum 
duty cycle. 

Studies of intermittent, randomly occurring events 
such as the electromagnetic radiation from lightning 
strokes require recording instruments of a special 
tvpe. Because of the excessive film consumption it 
is impractical to take photographic records of oscillo- 
scope waveforms for long periods of time with strip 
cameras. Commercially available framing cameras 
do not attain framing rates great enough for certain 
sferic observations. Frequently, it is necessary to 
combine the data with some form of accurate time 
presentation. For example, (1) sferics observed at 
widely separated stations can be identified on a 
world time basis, and (2) for precise measurement of 
the time interval between sferics. 

A device which appears to satisfy these require- 
ments has been developed at the National Bureau of 
Standards, Boulder Laboratories, Boulder, Colo. 
Although it was designed and built primarily for 
sferic studies, the basic concepts and certain features 
of this instrument may be useful in other research 
and practical applications (fig. 2). 

The timer is composed of conventional decimal- 
counting units arranged to count from 10 ysec to 31 
days with the count controlled by a precision 100 ke 
standard. Each of the decimal-counting units has 
a 1—-2-2—4 binary output, all of which are read out in 
parallel by the same pulse used to trigger the camera. 
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Figure 2. Functional system diagram. 


To prevent any ambiguity in the read-out due to the 
binary output changing state, the time of the read- 
out is controlled by the timer. The read-out is made 
by the first 100-ysec-pulse following the event to be 
recorded. 

Time is displayed for photographic purposes by 
use of small neon bulbs which are illuminated by a 
pulse of 100-ysec duration. Since this pulse is 100 
usec long, smaller increments of time must be re- 
solved by another means such as time markers on an 
oscilloscope trace. 


2. Design Approach 


It is assumed that an electronic trigger can be 
derived from the data to be recorded. The time at 
which the trigger occurs for the purposes of this 
report is considered to be the time of the event. 
This trigger and other triggers derived from it in 
fixed time relationships initiate the oscilloscope 
sweeps for the data presentation and provide the 
commands for the time read-out and the mechanical 
functions of the camera proper. 

Sixteen-millimeter film was chosen as a practical 
compromise between two mutually opposing con- 
siderations. Larger film would have been more 
desirable from the standpoint of picture quality or 
resolution but would have required greater torque to 
accelerate both it and the larger transport mech- 
anism. 

At the time the design was undertaken a film trans- 
port time of about 5 msec was considered adequate 
for the application at hand. In most cases, oscil- 
loscope sweeps of 500 to 1,000 ysec were used for 
data presentation. The objective was, therefore, to 
build a camera that would complete an operational 
cycle in about 6 msec. It was further required that 
the camera be able to cycle at this rate for a number 
of frames. 

In this application the average cycling rate was 
relatively low so there was no unusual requirement 
for a large driving motor. The essence of the 
mechanical problem was to supply the high instan- 
taneous torques necessary to accelerate and deceler- 
ate the film and especially the film transport mech- 
anism in the short period of time required. 

The peak versus average power requirements were 
met by the simple expedient of an appropriate fly- 
wheel mounted directly on the motor shaft. The 





more difficult part of the problem was to transmit 
the available high-peak power to the film in a reliable 
and uniform manner. 

The type of mechanism which would be ideal 
would combine some of the characteristics of both a 
friction and a positive action clutch. That is, a 
clutch which would engage and disengage smoothly 
and quickly with a minimum of external force and 
at the same time be capable of transmitting a very 
high torque. " 

The type of device selected was essentially a cone 
clutch with the engaging surface operating just 
outside the critical angle. 

The command signal available to actuate the 
camera (which was supplied by other equipment) 
was 2 usec in duration with a nominal amplitude of 
100 v. A transducer was required to cause this 
signal to engage the clutch for the period of time 
necessary to advance the film one frame. The 
principle used to accomplish this was to cause a fixed 
quantity of electricity to flow in a coil similar to the 
voice coil in a dynamic speaker upon receipt of each 
command signal. The movement of the voice coil 
was used directly to engage the clutch. 

The use of the usual sprocket wheel to advance the 
film was not considered because of excessive point 
pressures exerted on the film by the large accelerating 
forces. Instead, it was decided to apply the acceler- 
ating force across the entire width of the film by 
means of a precision neoprene-surfaced cylindrical 
roller. 

It was evident from the outset that the inertia of 
the intermittently moving parts would be of major 
importance. The design of the driven parts to 
minimize inertia and yet provide the necessary 
mechanical strength was rather straightforward. 
Managing the film supply, however, in such a manner 
that excessive inertial loads would be prevented was 
a separate problem of at least equal importance. In 
the design requirements it was estimated that the 
camera should be able to take at least 20 frames in 
rapid succession. It seemed totally impractical to 
draw film directly from the supply magazine at such 
high intermittent speeds because of the generally 
complicated mechanism involved and the possibility 
of film breakage. Accordingly, it was decided to 
provide and maintain a free and unsupported loop 
of film between the drive roller and the supply 
magazine. 

One simplifying aspect of the design requirements 
was that no shutter was required. In this case, the 
cathode-ray tubes were simply unblanked and the 
time was read out when the data occurred. 

A crystal-controlled timer which could be syn- 
chronized with standard time signals was chosen to 
maintain a time reference. In order to record the 
time of an event on the film an instantaneous and 
unambiguous parallel read-out scheme was devised. 
For circuit simplicity each increment of time greater 
than 200 usec was displayed on the film in a coded 
binary form. Increments of time smaller than 200 
usec were displayed on the oscilloscope traces. 
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The relative time accuracy of a crystal-controlled | 


timer is dependent on the stability of the crystal; 
whereas, the absolute time accuracy or the time 
comparison between any two geographical points is 


dependent on the accuracy of the synchronizing | 


signal. WWYV (and certain other standard time 
broadcasts) transmit world time accurately, but due 
to the propagation path, instantaneous measurements 
at distant locations can only be resolved to about 
i msec. Through arduous investigation and inter- 
pretation of available ionospheric data WWV time 
can probably be established in retrospect to within 
two or three tenths of a millisecond. 

A low-frequency navigation system known as 
Loran C transmits signals that may be received in 
excess of 1,000 miles with a stability of 1 to 2 usec. 
If these signals were available at all locations and if 
they were properly synchronized with each other and 
with a world time standard, absolute time would be 
available with an accuracy of 1 to 2 usec. 

The timer described in this report could be syn- 
chronized with Loran C signals to provide micro- 
second absolute timing. ; 


3. System Description 
3.1. Camera Mechanics 


The camera that was developed, based on the 
indicated design approach, utilizes methods for film 
storage and transfer that are rather unique and 
require detailed explanations. In addition to the 
film storage and transfer systems, the optical portion 
of the camera is a Dallmeyer, 15 mm wide-angle f1.5 
lens, mounted in a vertical plane for direct exposure. 

The main drive unit (fig. 3) consists of a 3,600- 
rpm induction drive motor, the flywheel, a conical 
clutch and brake, film drive rollers, and the trig- 
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FIGURE 3. Camera mechanical drive. 
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gering mechanism. Mounted directly on the motor 
is a flywheel with a nylon insert used as the driving 
element in the clutch brake system. The driven 
element, in normal position, is spring loaded against 
a stationary, identical nylon insert which acts as 
the braking element. The mechanical triggering 
mechanism €onsists of a large electromagnet and 
a freely suspended coil. In principle, it is entirely 
similar to a large electrodynamic speaker. When 
current is applied to the coil an axial thrust is 
exerted on the driven element, thus engaging it with 
the nylon insert (driving element) in the flywheel. 
The clutch and brake elements are tapered to form 
sections of a cone and the angle of this taper is close 
to the critical angle of the cone of friction. In this 
way, maximum frictional drive and. braking are 
obtained without wedging of the components. 

The kinetic energy stored in the flywheel provides 
the high instantaneous torque necessary to accelerate 
the transport mechanism when the clutch is engaged. 
The maximum duty cycle of the film advance is 
affected by the size of this flywheel since it must 
store sufficient energy so that the motor speed will 
not be significantly reduced during the maximum 
operating rate. 

The mechanism for engaging the clutch is a freely 
suspended coil of insulated copper wire wound on a 
3-in. diameter aluminum shell and centered in the 
air gap of the electromagnet. This coil is connected 
to the shaft of the driven element of the clutch in 
such a way that the shaft may rotate inside the 
coil, but the axial movement of the coil i§ transmitted 
to the shaft. The force exerted by this coil is pro- 
portional to flux density, wire length, art the cur- 
rent. The time during which the force i§-applied 
depends on the length of time the current flows; 
thus, the film advance is controlled by varying the 
period of time during which the current flows through 
the coil. 

The film transport mechanism consists of the film 
drive roller and a tension roller (fig. 4). .The film- 
drive roller is neoprene-coated and has a spline 
which engages with the shaft of the driven element 
of the clutch. The shaft is free to move axially but 
rotational motion is transmitted to the roller. The 
tension roller is spring loaded against the drive roller 
with provisions for release to facilitate film loading. 

The film storage mechanism contains a 200-ft 
supply and a 100-ft takeup magazine with a spring 
loaded light trap that operates when the magazines 
are removed from the camera. The takeup tension 
is applied by operating a small, heavy-duty motor 
(fig. 5) at low voltage in a stalled condition. This 
motor is located directly behind the takeup magazine 
and is conrected to the film spool inside this maga- 
zine. A set of metering rollers (fig. 4) maintain a 
loop of film to keep the shock of acceleration from 
tearing the film and to reduce the mass that the main- 
drive roller must move. These metering rollers are 
geared together in such a manner that, as film is 
pulled into the takeup magazine, an equal amount 
of film is removed from the supply magazine to keep 
a constant loop of film ahead of the drive roller. 
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Figure 4. Camera interior detail. 
a. Shock absorber arm, b. metering rollers, 


c. film loop, d. film drive roller, 
e. film tension roller, and f. optic system. 





Camera exterior detail. 


Figure 5. 


a. Takeup motor 


3.2. Camera Electronics 


The electronic circuitry for the camera is designed 
to deliver high current for short durations and at 
rapid intervals (fig. 6a). A large condenser reservoir 
supplies this high current intermittently and thus 
reduces current which must be drawn from a power 


supply to an average value sufficient to sustain the | 


average pulse rate. The reservoir presently used is a 
1,600-uf condenser which will provide approximately 
ten advances at maximum speed. The power supply 
is capable of delivering 275 ma which will support an 
average rate of 60 fps. The average rate could be 
increased by using a larger power supply, and the 
number of adjacent advances could be increased by 
increasing the 1,600-yf reservoir. The camera drive 
is actuated by discharging an 88-uf condenser with 
a GL-5544 thyratron through the ‘“‘voice”’ coil. The 
88-uf condenser is recharged from the large reservoir 
by use of a second GL-5544. 

The basic scheme is as follows: A main trigger 
pulse is applied to a thyratron trigger circuit. A 
40-usec pulse is generated and fed to the first thrya- 
tron through an isolation transformer. The main 
trigger is also delayed 2 msec and then used in the 
same way to fire the second thyratron. The second 
thyratron recharges the 88-uf condenser from the 
1,600-uf reservoir (fig. 6b). 
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Figure 6. Camera electronic circuitry. 


The complete electronic cycle requires less than 
5 msec, and since the film is moving for 5 msec, 
triggers spaced closer than this would be useless. 
Therefore, a “‘killer’’ circuit was included to eliminate 
triggers spaced closer than 5 msec. 


3.3. Timer Electronics 


A reference oscillator operating at 100 ke and with 
a stability of 1 part in 10° is used to control the timer. 
The 100-ke signal is used to form 10 usec pulses. 
These pulses are fed into a series of decimal-counting 
units so arranged that the count will correspond to 
time from microseconds to days. The resolution 
accuracy of this system is 2 or 3 usec and the time 
between any 2 sferics can be determined within this 
accuracy. The timer is adjusted by utilizing the 
WWY standard time signals. Therefore, the abso- 
lute time accuracy to which the timer is calibrated 
is limited by the propagation variables. This system 
could be synchronized more accurately if more pre- 
cise time signals were available. 
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The timer count from 10 usec to 10 sec is accom- 
plished by using standard decimal-counting units. 
Between each unit other stages were added to provide 
amplification, clipping, differentiation, and to isolate 
the output of one unit from the input to the next 
unit. The use of these stages in addition to cathode 
followers between the units permits longer leads for 
physical separation of the units. The limiting and 
isolation provide quite effective protection against 
spurious signal pickup. 

The count of six necessary for the seconds and 
minutes counters is accomplished by taking advan- 
tage of the fact that the output of the third stage of 
a decimal counter provides a count of six. The 
hours counter is arranged so that it will reset with 
an output pulse at 24 (fig. 7). This is accomplished 
by actuating a relay when the tens-of-hours counter 
reaches the count of two. When this relay is ener- 
gized the normal circuit in the units of bours counter 
is modified in that the output from the second binary 
digit at count of four is not fed to the third binary 
digit, therefore, this counter returns to zero rather 
than counting to five. The output of this second 
binary digit is applied to the second binary digit in 
the tens-of-hours counter, thus resetting it to zero 
and providing an output pulse to the days counters. 


3.4. Timer Calibration 


Dual inputs¥to the minutes, hours, and days 
counters facilitate resetting the counters after non- 
operating periods. One-second pulses can be sup- 
plied to all of these counters through the other input 


No. | No. 2 
MULTIVIBRATOR MULTIVIBRATOR 





for setting purposes. A push-button switch on the 
front of each unit allows the minutes, hours, or days 
counters to be reset within 60, 24, or 40 sec, respec- 
tively. The input to the seconds counter can be 
interrupted by a switch, holding that counter and 
hence all subsequent counters. Stopping the seconds 
counter at 00 and starting it again with the 5-min 
tone from WWYV completes the calibration to within 
1 sec. 

The procedure necessary to calibrate the timer 
more precisely is to trigger an oscilloscope with the 
1-sec pulses from the timer and to display the WWV 
signal on this oscilloscope (fig. 8). If the 1-see pulses 
are systematically slowed down, the WWV seconds 
pulses will appear to move from right to left across 
an oscilloscope with a normal sweep. To slow down 
the timer, two gate generators are also triggered with 
the 1-sec pulses, and these gates are used to remove 
pulses at the input to the l-msec counter. In the 
fast position 10 pulses are removed, slowing the timer 
10 msec per second; in the slow position 1 pulse is 
removed, slowing it 1 msec per second. This allows 
the timer to be set within 1 msec of WWV. If a 
more accurate setting were desired, the pulses could 
be removed at the input to the 100-usec counter, 
and a third gate with a duration between 100 usec 
and 1 msec could be used. 


3.5. Timer Read-Out 


The decimal-counting units have a modified binary 
read-out code. This code is the 1-2-2-4 code con- 
sistent with the binary states associated with the 
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Modification on hours counter. 
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Figure 8. Timer synchronizing circuit. 


feedback loops required to obtain a of 10 
rather than a count of 16. 

The output of each binary digit is fed to a logical 
Snr circuit which operates a neon bulb in the 

‘amera field when an additional gate input is applied. 
This gate input is initiated by the main trigger 
but controlled by the timer in such a way that all 
binary digits are read out in a stable state thus 
precluding ambiguity. Figure 9 shows the block 
diagram and the wave forms associated with the gate 


generator used to produce the gate for all of “the 


count 


“ond” circuits. The read command (fig. 9a) is 
applied to a gate generator. The generated gate 


(fig. 9b) passes the first 200-ysee pulse (fig. 9c). 
This pulse generates a second gate (fig. 9d) which is 
stopped simultaneously by the first gate (fig. 9b). 
Both gates terminate with the next 100-ysec pulse 
(fig. 9e). The 100-usee gate (fig. 9d) thus formed 
is the gate generator output which is one of the 
inputs to the “and” circuit. This gate occurs only 
when all of the binary counters above 200-ysec are 
in a stable state. 

The neon bulbs displaying the binary counts are 
arranged in sets of four for convenience in transposing 
to decimal numbers. The time below 200-yusec is 
resolved on the oscilloscope trace, since 10-and 200- 
usec markers are superimposed on the traces (fig. 
10). If traces of less than 200-yusec duration were 
desired, time could be read out on neon bulbs down 
to 10-ysec by utilizing circuitry to hold the binary 
digits between 10 and 200 usec. 
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FiGgurRE 10. Camera field presentation. 


4. System Performance and Results 


Extensive operation of the camera and timer dur- 
ing a 4-month period demonstrated the system’s 
usefulness and also indicated where modifications 
would be warranted. 

The direction that the film moved was parallel to 
the oscilloscope sweep direction (fig. 11). Sferies 
which occurred at time intervals less than 5 msec 
would overlap, and the traces would be lengthened. 
Rotation of the camera or the scope traces by 90° 
would resolve this problem. Movement of the film 
transversely to the direction of the trace should 
allow sferics separated by only 1 msec to be resolved. 





FicureE 11. 


Photographs of adjacent frames. 


The time measured from the main trigger until the 
camera advance is complete is 7 msec. The film 
motion is about 24 mm or one and one-half 16 mm 
frames. The calculated time necessary to move the 
film 24 mm with the 1-in. drive roller operating at 
3.600 rpm is 5.2 msec. The interval measured be- 
tween the main trigger and the start of the film ad- 
vance is 1.5 msec, so evidently when the film is 
moving it is moving at nearly maximum velocity. 
Figure 11 shows two frames 9 msec apart and 2 others 
with closer spacing (4.5 msec) where the film was 
moving when the second sferic occurred. It can be 
seen that if sferics occur with spacings greater than 
7 msec the camera is an intermittent-action frame 
camera, Whereas, if the sferics are more closely 
spaced the camera becomes a strip camera at a 
speed of 178 ips for short durations. It is obvious 
that camera triggers occurring more often than 5 
msec would be meaningless, and therefore a 5- 
msec blocking circuit is operated following each main 
trigger before the camera is triggered. 

Present limitations on the maximum film velocity 
are the motor speed and the drive-roller diameter. 
Increasing the motor speed creates bearing problems 





due to the heavy flywheel and shock loading. On 
the other hand, increasing the drive-roller diameter 
from 1 in. to 3 or 4 in. appears entirely feasible. 
This change would reduce the film transport time 
from 5.2 msec to 1.7 or 1.3 msec. 

Experiments have been conducted using an alumi- 
num disk and a coil wound with strip copper to de- 
termine the applicability of the induction disk to 
engage the clutch. The coil and the disk are coaxial 
with the film drive. A pulse of current through the 
coil pushes the disk over by repulsion. Initial re- 
sults indicate that this method may be faster, thus 
reducing the 1.5-msec start time, but the disk will 
not keep the clutch engaged for an extended period 
of time as can be done with the voice coil. 'The mech- 
anism associated with the induction disk has fewer 
moving parts and is inherently simpler. This method 
should be excellent for higher film transport speeds. 


Bou.pErR, CoLo. (Paper 64C2-36) 
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Publications of the National Bureau of Standards* 


Selected Abstracts 


A carbon-14 beta-ray standard, benzoic acid-7-C™ in 
toluene, for liquid scintillation counters, W. F. 
Marlow and R. W. Medlock, J. Research NBS 644A, 
No. 2, 143 (1960). 


A carbon-14 beta-ray standard for use in liquid scintillation 
counting has been prepared and standardized. The sample 
consists of benzoic acid-7-C" dissolved in toluene. Samples 
of the solution were oxidized quantitatively in a Paar oxygen 
bomb, and the radioactivity of the carbon dioxide was com- 
pared with the radioactivity of carbon dioxide prepared 
quantitatively from the Bureau’s sodium  carbonate-C' 
standard. 


Some effects of aging on the surface area of portland 
cement paste, C. M. Hunt, L. A. Tomes, and R. L. 
Blaine, J. Research NBS 64A, No. 2, 163 (1960). 


A hardened cement paste cured at room temperature, from 
which part of the evaporable water has been removed by 
vacuum drying, is shown to decrease in surface area with time. 
This change is the opposite of that usually observed during 
hydration and probably represents some colloidal growth 
phenomenon analogous to aging observed in other colloidal 
gels. Both water vapor and nitrogen absorption measure- 
ments show the effects of aging in cement paste. 


Wet or dry paste is shown to undergo less change than paste 
of intermediate evaporable water content, so that if surface 
area after storage is plotted as a function of evaporable water 
content, a curve with a minimum is obtained. With increasing 
storage temperature there is some indication that this mini- 
mum might shift towards lower water content. 


Aging is shown to occur during the initial drying of a cement 
paste, so that even the initial surface area of a cement paste 
depends upon the manner in which the paste has been dried. 


Some solutions for electromagnetic problems in- 
volving spheroidal, spherical, and cylindrical bodies, 
J. R. Wait, J. Research NBS 64B, No. 1, 15 (1960). 


Solutions are presented for the low-frequency electromagnetic 
response of an oscillating magnetic dipole by conducting 
bodies of simpie shape. The quasi-stationary approximation 
is employed throughout which is valid when the relevant 
dimensions of the problem are all small compared to the free- 
space wavelength. This amounts to matching solutions of the 
wave equation within the bodies to solutions of Laplace’s 
equation outside. 


Screw-thread standards for federal services 1957. 
Amends in part H28 (1944) (and in part its 1950 
Supplement). NBS Handb. H28 (1957)—Part II 
(1959) 75 cents. 


Information is presented on the limits of size, basic dimensions, 
thread form, and gaging of tapered and straight pipe threads 
(including Dryseal pipe threads), gas-cylinder-valve outlet 
and inlet threads, hose coupling threads (including fire-hose 


coupling threads), and welding and cutting torch hose 
connections. 


Communication theory aspects of television band- 
width conservation, W. C. Coombs, NBS Tech. Note 
25 (PB151384) (1959) 50 cents. 


New concepts of communication bandwidth utilization and 
conservation are in prospect through difference signal modu- 
lation systems in which only relative changes information is 
transmitted in lieu of absolute amplitudes. By these systems, 
the changes data function is made one of time as well as of 





amplitude, so that better advantage can be taken of redun- 
dancies in the video signal. 


Advantageous conversions of information rate are made 
feasible by reason of a tremendous disparity existing between 
the extremely high information capacity of conventional tele- 
vision systems and the relatively very low perception capa- 
bility of a human observer channel. Information conversions 
more consistent with perception capabilities of the human 
channel are achieved in difference signal modulation without 
usual great expansion of bandwidth that would be required 
to reproduce every possible absolute amplitude of each datum 
in the whole video picture mosaic directly. 


Conversion to binary digital form accrues the further ad- 
vantages of binary systems, including greater immunity to 
noise, greater adaptability to discrete data storage media, 
greater ease of scrambling for security, and greater amenability 
of encoding to optimum parametric controls for most ad- 
vantageous transmission. 


Distribution of mail by destination at the San Fran- 
cisco, Los Angeles, and Baltimore Post Offices, 
N.C. Severo and A. E. Newman, NBS Tech. Note 27 
(PB151386) (1959) $1.50. 


This report presents the results of the application of the 
sampling method presented in a paper by Norman C. Severo 
and Arthur E. Newman entitled “A statistical chain-ratio 
method for estimating relative volumes of mail to given 
destinations,’”’ to appear in NBS Journal of Research 64C, 
No. 1 (Jan.—Mar. 1960). These studies were made at the 
San Francisco, Los Angeles, and Baltimore Post Offices. 


Photographic dosimetry at total exposure levels 
below 20 mr, M. Ehrlich and W. L. McLaughlin, 
NBS Tech. Note 29 (PB151388) (1959) 50 cents. 


Photographic dosimetry at exposure levels below 20 mr is 
difficult because of insufficient sensitivity of commercial 
photographic emulsions, instabilities in the latent photo- 
graphic image, and inaccuracies in sensitometric procedure. 
Some of the endeavors to cover the range below 20 mr more 
adequately, both with film-scintillator combinations and 
with conventional photographie dosimeters, will be briefly 
discussed. 


Resistance diode bridge circuit for temperature 
control, L. H. Bennett and V. M. Johnson, NBS 
Tech. Note 34 (PB151393) (1959) 50 cents. 


The conventional a-c bridge gives irregular performance 
including loss of temperature control when the temperature 
error exceeds a certain critical value. The present note de- 
scribes a simple method of achieving stable temperature 
control over a large range of temperatures. 


Service area of an airborne television station, Ml. T. 
Decker, NBS Tech. Note 35 (PB151394) 75 cents. 


As a step in the evaluation of a proposed airborne television 
network, the service to be expected from an airborne station 
has been calculated for a wide variety of operating conditions. 
The use of basic transmission loss to describe the radio propa- 
gation effects enables the results to be used when equipment 
parameters are changed as system requirements and economic 
considerations dictate. 


Strain gauge calibration device for extreme tem- 
peratures, R. M. McClintock, Rev. Sei. Instr. 30, 
No. 8, 716 (1959). 

A device is described for applying known static strains to 
strain gages under extreme temperature environments. The 
essential part is a constant strength cantilever beam used in 
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Such a way that the mechanical properties of the component 
parts need not be accurately known. The location of the 
strain gages on the beam is not critical since strain is constant 
along its length. Conditions are given for which errors caused 
by thermal expansion can be eliminated. Measurements of 
the strain sensitivity of constantan wire gages are given at 
temperatures between 4° and 300° Kelvin. 


Rapid insertion device for coaxial attenuators, A. Y. 
Rumfelt and R. J. Como, Rev. Sci. Instr. 30, No. 8, 
687 (1959). 


A rapid insertion device for coaxial attenuators is described. 
This device is used as an alignment aid when making inser- 
tion loss measurements of fixed coaxial attenuators at fre- 
quencies up to 6 kMe. It provides a quick connect and dis- 
connect feature, maintains axial alignment of the connectors 
during joining and separation, and reduces cable flexure. 
The improved repeatability and accuracy of insertion loss 
measurements is shown by comparative data with different 
alignment techniques. 


PILOT, the NBS multicomputer system, A. L, 
Leiner, W. A. Notz, J. L. Smith, and A. Weinberger. 
Proc. Eastern Joint Computer Conf. (Philadelphia, 
Pa.), p. 164 (1958). 


At NBS a new large-scale digital system has been designed 
for use on a wide variety of experimental applications ranging 
from automatic search and interpretation of Patent Office 
records to real-time control of commercial aircraft traffic. 
Because of this wide range of intended use, the system had 
to combine in a single installation a variety of characteristics 
not ordinarily associated with one facility, namely (1) a high 
computation rate, (2) highly flexible control capabilities for 
communicating with the outside world, and (3) a wide reper- 
toire of internal processing formats. The over-all system is 
organized around three independently programmed computers 
which intercommunicate in such a way as to permit all three 
to work together concurrently on a common problem, each 
computer working on the piece of the problem for which it is 
individually best suited. The over-all system thus provides 
a working model of an integrated multi-computer network. 


Investigation of creep behavior of structural joints 
under cyclic loads and temperatures, L.. Mordfin, 
N. Halsey, and G. E. Greene, NASA Tech. Note 
D-181, p. 37 (1959). 


EKighty-two structural joint specimens were tested to evaluate 
the effects of cyclic loads and cyclic temperatures on creep 
and rupture. The specimens included riveted joints of 
2024-T3 clad aluminum alloy, and riveted and spot-welded 
joints of 17-7 PH (TH 1050) stainless steel. The results of 
these tests show a wide variance but indicate certain trends 
which permit the estimation of the cyclic creep behavior of 
joints. 

An analysis of the effects of stress concentration on the tensile 
rupture strength of riveted joints is presented in appendix A. 
This, together with data from several other laboratories, 
shows that the effects are small for joints fabricated from 
notch ductile materials and conventional rivets. 


Microwave reflectometer techniques, G. F. Engen 
and R. W. Beatty, JRE Trans. on Microw. Theory 
and Tech., MTT-7, No. 3, 346 (1959). 


A rigorous analysis of the microwave reflectometer is pre- 
sented for what is believed to be the first time. By means 
of this analysis, the correct adjustment of auxiliary tuners is 
described and the errors resulting from incorrect adjustments 
are treated in a quantitative manner. 

It is shown how the reflectometer technique may be further 
simplified while preserving the accuracy of measurement. 
A convenient method of adjusting the auxiliary tuners is 
described, sources of error are discussed, and an example is 
given of the calculation of error limits. 


Cryogenic engineering of hydrogen bubble cham- 
bers, B. W. Birmingham, D. B. Chelton, D. B. Mann, 
and H. P. Hernandez, AS7T’M Bull. No. 240, 34 
(TP164) (1959). 


Low-temperature material problems encountered during the 
development of a liquid-hydrogen bubble chamber are dis- 
cussed. Properties affecting the selection of the bubble 
chamber material, the glass window, and the thermal radiation 
shield are examined. Methods of sealing the glass window 
to the chamber with various gaskets suitable for use at liquid 
hydrogen temperature (—423 F) are presented. In addition, 
a method of continuously maintaining the bubble chamber 
and contents at liquid hydrogen temperature by means of an 
external refrigerator is described. 


On prediction of system behavior, J. R. Rosenblatt, 
Proc. of the New Yrok Univ.-Ind. Conf. on Reliability 
Theory (New York, N.Y.), p. 89 (1958). 


It is customary in theoretical discussions of prediction of 
system performance to assume that the relation between 
variables by which performance is assessed and variables 
describing the parts of the system can be specified by a fune- 
tion. The purpose of this paper is to consider some of the 
problems which arise in choosing the form of such a function 
and in specifying the relevant variables and the nature of 
appropriate experiments, measurements, and data. 

In the context of each of a number of familiar types of math- 
ematical model, the possible consequences of alternative 
choices of time units, subsystem definitions, and independence 
assumptions are discussed. 


Ceramic wafer tubes for modular units, (. P. 
Marsden, sec. on Components and Materials, Flec- 
tronics 32, No. 39, 94 (1959). 


Mechanized assembly and processing means have been 
designed and constructed to show that ceramic wafer tubes 
can be produced with a vield of better than 75 percent and 
with a high degree of reliability and ruggedness in military 
environments. This article discloses the design, methods of 
fabrication and the processing of these tubes. 


Magnified and squared VSWR responses for micro- 
wave reflection coefficient measurements, R. W. 
Beatty, JRE Trans. on Microw. Theory and Tech. 


MTT-7, No. 3, 346 (1959). 


In conventional microwave impedance measuring instru- 
ments, the measured ratio of maximum to minimum detector 
signal level is ideally equal to the voltage standing-wave 
ratio (VSWR) of the termination. In this paper, it is shown 
how radically different types of response are obtainable in 
which the observed ratio may approximately equal the square 
of the VSWR or may be magnified any desired amount. 
Theory is given enabling accurate measurements by interest- 
ing techniques. Accuracies of 0.1 pereent in VSWR to 2.0 
have been achieved using magnified response techniques. 


Method of cooling head-on photomultipliers, G. C. 
Harman, Rev. Sci. Instr. 30, No. 8, 742 (1959). 


A simple method of cooling the photoemissive surface of 
head-on photomultipliers is described. It consists of a small, 
specially constructed Dewar with a window in its top. The 
Dewar is charged with dry ice powder and slipped over the 
head of the multiplier so that it cools only the sensitive 
surface. The rest of the tube remains at room temperature. 
The dark current is reduced by a factor of at least 500. 


Optical measurements on thin films of condensed 
gases at low temperatures, J. Kruger and W. J. 
Ambs, J. Opt. Soe. Am. 44, 1195 (1959). 

A determination of the refractive indices of thin films of 


oxygen, nitrogen, carbon dioxide, water, argon, neon, and 
krypton condensed at 4.2° K was made. These were ob- 
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tained for the 5461 A line of Hg with an ellipsometer. Simi- 
lar measurements were also made on films of these gases 
condensed after passage through a microwave discharge. 
With the exception of N», and Kr measurable differences in 
refractive indices were observed between films condensed 
from discharged and undischarged gases. 


Wavelength definition of the meter, I. C. Gardner, 
Systems of Units, Natl. and Intern. Aspects, Publ. 
Am. Assoc. Advance. Sci., No. 57, p. 58 (1959). 


The different applications of the national prototype meter to 
the calibration of line standards, of end standards, and of 
wavelengths are presented. The advantages of a wavelength 
standard are considered and the progress of international 
conferences for redefining the meter in terms of a wavelength 
standard is set forth. One may expect final and conclusive 
action at the Conference (General) of Weights and Measures 
in 1960. The advantages to be anticipated from the devel- 
opment of an atomic beam source are discussed. 


Digital recording of electrocardiographic data for 


analysis by a _ digital computer, L. aback, 
E. Marden, H. L. Mason, and H. V. Pipberger, 


IRE Trans. Med. Electronics ME-6, 167 (1959). 


A corrected orthogonal 3-lead system has been used to record 
electrocardiograms directly from patients at Veterans 
Hospitals, using three F M channels of magnetic tape. A pilot 
facility has been designed and assembled by NBS to permit a 
medical technician to inspect these on an oscilloscope and 
select a significant cardiac cycle. This is automatically 
sampled at millisecond intervals and the numerical values 
are stored in digital form on magnetic tape acceptable to an 
electronic computer. Upon writing various programs for 
the digital computer, the cardiac researcher will have a 
flexible tool for objective analysis of large quantities of 
biological data by a variety of possible criteria. 


Young’s modulus of various refractory materials as 
a function of temperature, J. B. Wachtman, Jr., 
and D. G. Lam, Jr., J. Am. Ceram. Soc. 42, No. 5, 
254 (1959). 


Young’s modulus as a funetion of temperature was deter- 
mined by a dynamic method for single-crystal sapphire and 
ruby and for polycrystalline aluminum oxide, magnesium 
oxide, thorium oxide, mullite, spinel, stabilized zirconium 
oxide, silicon carbide, and nickel-bonded titanium carbide. 
For the single crystals, Young’s modulus was found to de- 
crease linearly with increasing temperature from 100° C to 
the highest temperature of measurement. For all the poly- 
crystalline materials, except silicon carbide, stabilized zir- 
conium oxide, and spinel, Young’s modulus was found to 
decrease approximately linearly with increasing temperature 
until some temperature range characteristic of the material 
was reached in which Young’s modulus decreased very rapidly 
and in a nonlinear manner with increasing temperature. 
This rapid decrease at high temperature is attributed to 
grain-boundary slip. Stabilized zirconium oxide and spinel 
were found to have the same rapid decrease in Young’s 
modulus at high temperature, but they also had a decidely 
nonlinear temperature dependence at low temperature. 


Thermodynamic properties of helium at low tem- 
peratures and high pressures, D. B. Mann and 
R. B. Stewart, J. Heat Transfer 81, 323 (1959). 


The thermodynamic properties of helium have been com- 
piled and correlated for a temperature range from 3.0 to 
20° K for pressures to 100 atmospheres and for specific 
volumes from 5 to 800 liters per kilogram. The properties 
are presented on both the Temperature-Entropy and the 
Enthalpy-Entropy coordinate systems and include pressure, 
temperature, volume, entropy and enthalpy. 
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Fire research at the National Bureau of Standards, 
A. F. Robertson, Fire Research Abstr. and Rev. I, 
No. 4, 159 (1959). 


A review is presented of work performed at the National 
Bureau of Standards and primarily directed towards the 
unwanted or accidental fire problem. It is emphasized that 
the general field of fire research is so broad that many other 
research activities at the NBS may be considered as perti- 
nent to this subject. Some few of these studies are men- 
tioned; however, emphasis is placed on problems relating to 
ignition, growth, confinement, and extinguishment of fires. 
Pertinent references to NBS work are included. 


Preliminary assessment of the IGY, A. H. Shapley, 
Proc. Natl. Electron. Conf., p. 1 (1958). 


The activities of the International Geophysical Year are 
reviewed, with particular emphasis on studies of the 
ionosphere. 


Factorial experiments in life testing, M. Zelen, 
Technometries 1, No. 3, 269 (1959). 


This paper discusses procedures for analyzing factorial ex- 
periments, where the experiment deals with the life testing 
of components or equipment. These procedures assume an 
underlying general distribution of ‘‘times-to-failure,”’ of 
which the exponential, Weibull, and extreme value distribu- 
tions are special cases. Statistical tests and confidence pro- 
cedures are outlined, and an example illustrating the proce- 
dure for life-test results of glass capacitors is included. 
Small sample approximations, which are adequate for prac- 
tical applications, are given for the proposed procedures. 
This is shown empirically by generating thousands of life- 
test experiments on an electronic computer. An empirical 
sampling investigation is given of the robustness of the pro- 
posed procedures. From the sampling results, it is concluded 
that these techniques are sensitive (non-robust) to departures 
from the original assumptions on the probability distribution 
of failure-times. An investigation is also given of a trans- 
formation which appears to give robust results. These same 
techniques carry over exactly to the situation where one is 
analyzing an array of variance estimates from an underlying 
normal population. 


Other NBS Publications 


Journal of Research, Section A. Physics and 
Chemistry, Volume 64A, No. 2, March-April 1960. 
70 cents. 


A carbon-14 beta-ray standard, benzoic acid-7-C" in toluene, 
for liquid scintillation counters. W. F. Marlow and R. W. 
Medlock. (See above abstracts.) 

A comparison of experimental and theoretical relations 
between Young’s modulus and the flexural and longitudinal 
resonance frequencies of uniform bars. 8. Spinner, T. W. 
Reichard, and W. E. Tefft. 

Determination of copolymer composition 
analysis for carbon and hydrogen. 
Irving Madorsky, and Rolf A. Paulson. 

Some effects of aging on the surface area of Portland cement 
paste. C. M. Hunt, L. A. Tomes, and R. L. Blaine. 
(See above abstracts.) 

Conformations of the pyranoid sugars. I. 
conformers. H.S. Isbell and R. 8. Tipson. 

Tritium-labeled compounds III. Aldoses-/-i. H. 
H. L. Frush, N. B. Holt, and J. D. Moyer. 

Determination of aluminum in precipitation hardening 
stainless steel and high temperature alloys. Lawrence A. 
Machlan, John L. Hague, and Edward J. Meros. 


by combustion 
Lawrence A. Wood, 


Classification of 


S. Isbell, 
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Journal of Research, Section B. Mathematics and 
Mathematical Physics, Vol. 64B, No. 1, January-— 
March 1960. 75 cents. 


Theory of diffraction in microwave interferometry. D. M. 
Kerns and E. 8. Dayhoff. 

Some solutions for electromagnetic problems 
spheroidal, spherical, and cylindrical bodies. 
Wait. (See above abstracts.) 

Kantorovich’s inequality. Morris Newman. 

A symmetric continuous poker model. A. J. Goldman and J. 
J. Stone. 

Moebius function on the lattice of dense subgraphs. R. E. 
Nettleton and M: 8. Green. 

The minimum of a certain linear form. Karl Goldberg. 

Space of k-commutative matrices. Marvin Marcus and N. 
A. Khan. 

Selected bibliography of statistical literature, 1930 to 1957: I. 
Correlation and regression theory. Lola 8. Deming. 

Selected bibliography of statistical literature, 1930 to 1957: IT. 
Time series. Lola S. Deming. 


involving 
James R. 


Journal of Research, Section D. Radio Propaga- 
tion, Volume 64D, No. 3, May-June 1960. 70 cents. 


A theory of radar scattering by the moon. T. B. A. Senior 
and K. M. Siegel. 

A theory of wavelength dependence in ultrahigh frequency 
transhorizon propagation based on meteorological con- 
siderations. Ralph Bolgiano, Jr. 

A preliminary study of radiometeorological effects on beyond- 
horizon propagation. F. Ikegami. 

Trade-wind inversion as a transoceanic duct. M. Katzin, 
H. Pezzner, B. Y.—C. Koo, J. V. Larson, and J. C. Katzin. 

An analysis of propagation measurements made at 418 
megacycles per second well beyond the radio horizon (a 
digest). H. B. Janes, J. C. Stroud, and M. T. Decker. 

On the calculation of the departures of radio wave bending 
from normal. B. R. Bean and E. J. Dutton. 

On the mode theory of very-low-frequency propagation in 
the presence of a transverse magnetic field. D. D. Crombie. 

On the theory of reflection of low- and very-low-radiofre- 
quency waves from the ionosphere. J. Ralph Johler and 
Lillie C. Walters. 

Focusing, defocusing, and refraction in a circularly stratified 
atmosphere. K. Toman. 

Response of a loaded electric dipole in an imperfectly con- 
ducting cylinder of finite length. Charles W. Harrison, 
Jr., and Ronold W. P. King. 

Impedance characteristics of a uniform current loop having 
aspherical core. Saburo Adachi. 


Hydraulic Research in the United States—1959, H. K. Middle- 
ton, Mise. Pub. 227 (1959) $1.25. 

Frequency dependence of VHF ionospheric scattering, J. C. 
Blair, NBS Tech. Note 9 (PB151368) (1959) 75 cents. 

Some applications of statistical sampling methods to outgoing 
letter mail characteristics, N. C. Severo, A. E. Newman, 
S. M. Young, and M. Zelen, NBS Tech. Note 16 (PB151- 
375) (1959) $2.75. 

Radio noise data for the International Geophysical Year 
July 1, 1957—December 31, 1958, W. Q. Crichlow, C. A. 
Samson, R. T. Disney, and M. A. Jenkins, NBS Tech. 
Note 18 (PB151377) (1959) $2.50. 

Variations of gamma cassiopeiae, 8. R. Pottasch, NBS Tech. 
Note 21 (PB151380) (1959) 75 cents. 

A history of vertical-incidence ionosphere sounding at the 
National Bureau of Standards, 8S. C. Gladdern, NBS Tech. 
Note 28 (PB151387) (1959) $2.00. 

Aerodynamic phenomena in stellar atmospheres—A bibliogra- 
phy, NBS Tech. Note 30 (PB151389) (1959) $1.25. 





The following papers appear in Treatise on Anal. Chem. 
pt. 1. Theory and practice, vol. 1 (The Interscience 
Encyclopedia, Ine., New York, N.Y., 1959): 

Accuracy and precision: Evaluation and interpretation of 
analytical data, W. J. Youden, ch. 3. See. A. Analytical 
chemistry: its objectives, functions and limitations, p. 


+. 

Principles and methods of sampling, W. W. Walton and 
J. 1. Hoffman, ch. 4. See. A. Analytical chemistry: 
its objectives, functions and limitations, p. 67. ; 

Atomic weights, E. Wichers, ch. 6. Sec. B. Application of 
chemical principles, p. 161. 

Concept and determination of pH, R. G. Bates, ch. 160. 
Sec. B. Application of chemical principles, p. 361. 

Electrode potentials, R. G. Bates, ch. 9. See. B. Applica- 
tion of chemical principles, p. 319. 


System loss in radio-wave propagation, K. A. Norton, Letter 
to Editor, Proc. IRE 47, No. 9, 1661 (1959). 

The oxide films formed on copper single crystal surfaces in 
pure water. 1. Nature of the films formed at room tem- 
perature, J. Kruger, J. Electrochem. Soc. 106, No. 10, 
847 (1959). 

Variation of the thermodynamic ideal temperature in the 
potystyrene-cyclohexane system, D. MelIntyre, J. H. 
O’ Marsa, and B. C. Konouck, J. Am. Chem. Soc. 81, 3498 
(1959). 

Energy dissipation in standing waves in rectangular basins, 
G. H. Keulegan, J. Fluid Mech. 6, Pt. 1, 35 (1959). 

Trapped energetic radicals, J. L. Franklin and H. P. Broida, 
Ann. Rev. Phys. Chem. 10, 145 (1959). 

Total photoelectric cross sections of copper, molybdenum, 
silver, tantalum, and gold at 662 kev, W. F. Titus, Phys. 
Rev. 115, No. 2, 351 (1959). 

Diode in feedback loop makes stable transistor bias circuit, 
G. F. Montgomery, Electronic Design 7, 54 (1959). 

Strength of synthetic single crystal sapphire and ruby as a 
function of temperature and orientation, J. B. Wachtman 
and L. H. Maxwell, J. Am. Ceram. Soc. 42, No. 9, 432 
(1959). 

Dimensional stability of impregnated sole leather, B. H. 
Fouquet, J. Am. Leather Chemists’ Assoc. LIV, No. 10, 
544 (1959). 

Theoretical problems of reliability measurement and_pre- 
diction, J. R. Rosenblatt, Proce. of the Third Exploratory 
Conf. on Missile Model Design for Reliability Prediction 
at White Sands Missile Range (New Mexico) p. 29 (1959). 

Dimensional changes in systems of fibrous macromolecules: 
polyethylene, L. Mandelkern, D. E. Roberts, A. F. Diorio, 
and A.§. Posner, J. Am. Chem. Soc. 81, 4148 (1959). 

Attenuation of scattered cesium-137 gamma rays, F. §. 
Frantz, Jr., and H. O. Wyckoff, Radiology 73, No. 2, 
263 (1959). 

Observations of the ionosphere over the south geographic 
pole, R. W. Knecht, J. Geophys. Research 64, No. 9, 
1243 (1959). 

Dependence of mechanical strength of brittle polycrystalline 
specimens on porosity and grain size, F. P. Knudsen, J. 
Am. Ceram. Soc. 42, No. 8, 376 (1959). 

Theory of dielectric relaxation in molecular crystals, J. D. 
Hoffman, Arch. sci. (Geneva), Soe. Phys. et Hist. Nat. 
Geneve 12, 36 (1959). 

Denture reliners—direct, hard, self-curing resin, G. M. 
Brauer, E. E. White, Jr., C. L. Burns, and J. B. Woelfel, 
J. Am. Dental Assoc. 59, 270 (1959). 

Photoproton cross sections of carbon, S. Penner and J. E. 
Leiss, Phys. Rev. 114, No. 4, 1101 (1959). 

A porcelain reference electrode conductive to sodium ions 
for use in molten salt systems, R. J. Labrie and V. A. Lamb, 
J. Electrochem. Soc. 106, No. 10, 895 (1959). 

Memorandum on a procedure for obtaining spectral radiant 
intensities of tungsten-filament lamps, 400-700 mu, L. E. 
Barbrow, Letter to Editor, J. Opt. Soc. Am. 49, No. 11, 
1122 (1959). 

Paramagnetic resonance spectra of active species. Blue 
material from hydrazoic acid, W. B. Gager and F. O. Rice, 
J. Chem. Phys. 31, No. 2, 564 (1959). 

Argon deposition on a 4.2° K surface, S. N. Foner, F. A. 
Mauer, and L. H. Bolz, J. Chem. Phys. 31, No. 2, 546 
(1959). 
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Low-temperature tensile properties of copper and four bronzes, 
R. M. McClintock, D. A. Van Gundy, and R. H. Kropschot, 
ASTM Bull. No. 240, 34 (TP177) (1959). 

Kinetic equation for a plasma with unsteady correlations, 
Cc. M. Tchen, Phys. Rev. 114, 394 (1959). 

A nomenclature for conformations of pyranoid sugars and 
derivatives, H. 8S. Isbell and R. 8. Tipson, Science 130, 
No. 3378, 793 (1959). 

On the convergence of the Rayleigh quotient iteration for 
the computation of the characteristic roots and vectors. 
III (General Rayleigh quotient and characteristic roots 
with linear elementary divisors) A. M. Ostrowski, Arch. 
Rat. Mech, Anal. 3, 325 (1959). 

On the convergence of the Rayleigh quotient iteration for 
the computation of the charactristic roots and vectors. 
IV (Generalized Rayleigh quotient for nonlinear elementary 
divisors), A. M. Ostrowski, Arch. Rat. Mech. Anal. 3, 
341 (1959). 

Spectral study of a visible, short-duration afterglow in 
nitrogen, B. E. Beale, Jr., and H. P. Broida, J. Chem. 
Phys. 31, No. 4, 1030 (1959). 

Atomic electron affinities, F. Rohrlich, Nature 183, 244 (1959) 

Cavity resonators for dielectric spectroscopy of compressed 
gases, H. E. Bussey and G. Birnbaum, Rev. Sci. Instr. 30, 
800 (1959). 

A discussion of federal specifications GS-X-—620 and L—F-310 
for dental X-ray apparatus and dental X-ray film, G. C. 
Paffenbarger, A. F. Forziati, and M. P. Kumpula, J. Am. 
Dental Assoc. 59, 472 (1959). 

Dynamic stability of frozen radicals. 
of the model, J. 
722 (1959). 

An additional lunar influence on equatorial Es at Huancayo, 
R. W. Knecht, Research Note J. Atmospheric and Terrest. 
Phys. 14, 348 (1959). 

On artificial geomagnetic and ionospheric storms associated 
with high-altitude explosions, S. Matsusbita, J. Geophys. 
Research 6, 1149 (1959). 

An intermittent-action camera with absolute time calibration, 
R. G. Hefley, R. H. Doherty, and E. L. Berger, IRE 
Wescon Conv. Record 3, Pt. 6, 129 (1959). 

Nonresonant microwave absorption and electric dipole 
moment of NO in the gaseous state, A. A. Maryott and 
S. J. Kryder, J. Chem. Phys. 31, No. 3, 617 (1959). 


II. The formal theory 
L. Jackson, J. Chem. Phys. 31, No. 3. 


The earth and its environment, 8. Chapman, Proce. IRE 
47, No. 2, 137 (1959). 

Bremsstrahlung cross-section formulas and related data, 
H. W. Koch and J. W. Motz, Rev. Modern Phys. 31, 


No. 4, 920 (1959). 

Note on quiet-day vertical cross sections of the ionosphere 
along 75° W geographic meridian, J. W. Wright, Letter 
J. Geophys. Research 64, 1631 (1959). 

Effect of atomic tests on radio noise, C. A. Samson, Letter 
Nature 184, 538 (1959). 

Short-time stability of a quartz-crystal oscillator as measured 
with an ammonia maser, A. H. Morgan and J. A. Barnes, 
Letter Proce. IRE 49, 1782 (1959). 

Electromagnetic radiation from cylindrical structures, J. R. 
Wait (Pergamon Press, New York, N.Y., 1959). 

Interplanetary gas. I: Hydrogen radiation in the night sky, 
J. C. Brandt and J. W. Chamberlain, Astrophys. J. 130, 
670 (1959). 

The nova outburst. III. The ionization of hydrogen gas by 
an exciting star, J. Jefferies and 8S. Pottasch, Ann. Astrophys. 
J. 2%, 318 (1959). 

Exploratory study, by low temperature X-ray diffraction 
techniques, of diborane and the products of a microwave 
discharge in diborane, L. H. Bolz, F. A. Mauer, and H. 8. 
Peiser, J. Chem. Phys. 31, No. 4, 1005 (1959). 

Some clinical applications of research findings in dental 
materials, G. C. Paffenbarger, Obio Dental J. 33, No. 3, 
218 (1959). 

Liquid helium cryostat with an integral super-conducting 
resonator, EK. Maxwell and F. A. Schmidt, Suppl. Bull. 
Inst. Intern. du Froid, Comm. 1, Delft, Holland, Annexe 
1958-1, p. 95 (1958). 

Sporadic E observed on VHF oblique incidence circuits, E. K. 
Smith, In North Atlantic Treaty Organization. Sporadic 
E ionization; Ionospheric Research Meeting, AGARD 





Avionics Panel, Cambridge, England, 
AGAR Dograph 34), p. 129 (1958). 

Low-temperature transport properties of copper and its dilute 
alloys: Pure copper, annealed and cold-drawn, R. L. 
Powell, H. M. Roder, and W. J. Hall, Phys. Rev. 115, 
314 (1959). 

The double bond isomerization of olefins by hydrogen atoms 
at —195°, M. D. Scheer and R. Klein, J. Phys. Chem. 
63, 1517 (1959). 

Paramagnetic resonance in the free hydroxyl radical, H. E. 
Radford, Nuovo Cimento, Serie X 14, 245 (1959). 

Relaxation processes in multistate systems, K. E. 
Phys. of Fluids 2, No. 4, 442 (1959). 


Sept. 1958 


Shuler, 


The following appear in Proc. Intern. Rubber Conf., Nov. 8 

to 13, 1959 (Washington, D.C.): 

Power loss and operating temperature of tires, R. D. 
Stiehler, M. N. Steel, G. G. Richey, J. Mandel, and 
R. H. Hobbs, p. 73. 

An indoor tester for measuring tire tread wear, G. G. 
Richey, J. Mandel, and R. D. Stiehler, p. 104. 

Measurement of the aging of rubber‘ vulecanizates, J. 
Mandel, F. L. Roth, M. N. Steel, and R. D. Stiehler, 
p. 221. 

Standard materials for rubber compounding, F. L. 


Roth 
and R. D. Stiehler, Proce. p. 232. 


Invariant and complete stress functions for general continua, 
C. Truesdell, Arch. for Rational Mech. and Anal. 1, No. 4, 
1 (1959). 

On the convergence of Gauss’ alternating procedure in the 
method of least squares, A. M. Ostrowski, Ann. Mat. Pura 
Appl. (Bologna, Italy) (IV), 48, 229 (1959). 

Ethane carbon-carbon distance obtained from _ infrared 
spectra, H. C. Allen and E. K. Plyler, J. Chem. Phys. 31, 
No. 4, 1062 (1959). 

The anomalous inversion in cristobalite, R. F. Walker, The 
kineties of high temperature processes, 228 pages (J. Wiley 
& Sons, New York, N.Y., 1959). 

Determination of starch in paper. A comparison of the 
TAPPI, enzymatic, and spectrophotometric methods, J. L. 
Harvey, B. W. Forshee, and D. G. Fletcher, Tappi 42, 
No. 11, 878 (1959). 

Branched-chain higher sugars. III. A 4-C-(hydroxymethyl)- 
pentose, R. Schaffer, J. Am. Chem. Soc. 81, 5452 (1959). 
Vibrational spectrum of cyanate ion in various alkali halide 
lattices, A. Maki and J. C. Decius, J. Chem. Phys. 31, 

No. 3, 772 (1959). } 

Glass research at the National Bureau of Standards, C. H. 
Hahner, The Glass Industry, p. 588 (1959). 

Applications of the molecular refractivity in radio metrology, 
B. R. Bean and R. M. Gallet, J. Geophys. Research 64, 
No. 10, 1439 (1959). 

Use of crystal to display high energy X-ray images, J. § 
Pruitt, Non-Destructive Testing XVII, No. 6, 359 (1959). 

The measuring process, J. Mandel, Technometrics 1, No. 3, 
251 (1959). 

Refractive indices and transmittances of several optical 
glasses in the infrared, G. W. Cleek, J. J. Villa, and C. H. 
Hahner, J. Opt. Soc. Am. 49, No. 11, 1090 (1959). 

Absolute photometry of the aurora—II. Molecular nitrogen 
emission in the sunlit atmosphere, M. H. Rees, J. Atmos- 
pheric and Terrest. Phys. 14, 338 (1959). 

Absolute photometry of the aurora—I. The ionized molec- 
ular nitrogen emission and the oxygen green line in the 
dark atmosphere, M. H. Rees, J. Atmospheric and Terrest. 
Phys. 14, 325 (1959). 


*Publications for which a price is indicated (except for NBS 
Technical Notes) are avaliable only from the Superintendent of 
Documents, U.S. Government Printing Office, Washington 25, 
D.C. (foreign postage, one-fourth additional). Technical Notes 
are available only from the Office of Technical Services, U.S. 
Department of Commerce, Washington 25, D.C. (Order by PB 
number.) Reprints from outside journals and the NBS Journal 
of Research may often be obtained directly from the author. 
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